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NOTICE TO AUTHORS 


1, Communications. Papers must be communicated to the Society by a Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the r, and drawing attention to important new information and to the main 
conclusions. summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in length. Authors are requested to submit MSS. in duplicate. These should be 

using double spacing and leaving a margin of not less than one inch on the 
left-hand side. Corrections to the MSS. should be made in the text and not in the 
margin. Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.__Authors are allowed considerable latitude, but they are requeated to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively : 

A. Corlin, Zs, f. Astrophys., 1g, 239, 1938. 

H. Jeffreys, Theory of Probability, and edn., section 5.45, p. 258, Oxford, 1948. 


3. Notation.—Vor technical astronomical terms, authors should conform closely to the 
recommendations of Commission 34 of the International Astronomical Union (Trans. 
1.A.U., Vol. VI, p. 345, 1938). Council has decided to adopt the I.A.U. 3-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, p. 221, 1932). In general 
matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(London; Royal Society, 1951) except where these conflict with 1.A.U. practice, 


4. Diagrams.--These should be designed to appear upright on the page, drawn 
about eoles the size required in print and prepared for direct photographic 
reproduction except for the lettering, which should be inserted in pencil. 
Legends should be given in the manuscript indicating where in the text the 
figure should ar. Blocks are retained by the Society for 10 years; unless the author 
requires ther before the end of this period they are then destroyed. 


5. Tables.-These should be arranged so that they can be printed upright on 
the page. 


6. Proofs.-Costs of alteration exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society's interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts,--When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MSS. 


Reading of Papers at Meetings 
8. When submitti yom authors are requested to indicate whether they will be 
willing and able to ual the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 
g. Postcards giving the programme of cach meeting are issued some days before the 
meeting concerned. Fellows wishing to receive such cards whether for Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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MEETING OF 1955 MARCH 11 
Professor Sir Harold Jeffreys, President, in the Chair 


‘The election by the Council of the following Fellow was duly confirmed ; 
Valdemar Axel Firsoff, M.A., ‘The Retreat, Pilton, Shepton Mallet, Somerset 
(proposed by P. Moore). 
Ninety-three presents were announced as having been received since the last 
Meeting, including : 
A. C. B. Lovell, Meteor Astronomy (presented by the author); and 
P. A. Moore, Guide to the Planets (presented by the author). 


MEETING OF 1955 APRIL 6 
Protessor Sir Harold Jeffreys, President, in the Chair 


‘The election by the Council of the following Fellows was duly confirmed : 

lan Robin Holton Brickett, 69 Dunvegan Street, Sydenham, Johannesburg, 
South Africa (proposed by W. H. van den Bos); 

George Stephen Brosan, Ph.D., B.Sc.(Eng)., A.M.1.E.E., A.M.1.Prod.k., 
The Polytechnic, Regent Street, London, W.1 (proposed by P. Moore); 

Walter George Grimwood, 21 Southlands Road, Bexhill-on-Sea, Sussex 
(proposed by L. S. 'T. Symms); 

Hugh M. Johnson, Ph.D., State University of lowa, lowa City, lowa, 
U.S.A. (proposed by R. H. Garstang); and 

Charles Alexander Shain, B.Sc., Radiophysics Laboratory, C.S.1.R.O., 
University Grounds, City Road, Chippendale, N.S.W., Australia 
(proposed by J. L. Pawsey). 

‘The President announced that the Council had elected the following an 

Associate of the Society : 

Donald Howard Menzel, Professor of Astrophysics and Director of the 
Harvard College Observatory, Cambridge, Massachusetts, U.S.A. 

Kighty-one presents were announced as having been received since the last 

Meeting, including : 

A. Dauvillier, Le magnétisme des corps célestes, Volume | (presented by the 
author); and 

G. J. Whitrow, La Structure de I’ Univers—revised French edition, translated 
by G. de Vaucouleurs (presented by the author). 
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Professor Dirk Brouwer delivered the George Darwin Lecture (see p. 221), 
taking as his subject “ The motions of the outer planets’’. 


After the lecture, the President delivered a short address, and presented the 
Gold Medal of the Society to Professor Brouwer, to whom it had been awarded 
for his outstanding contributions to celestial mechanics. ‘The President said : 


“Most of you were present in February when Dr. Jackson explained the 
Council's reasons for awarding the Gold Medal to Professor Dirk Brouwer. As 
1 was not on the Council at the time, it may be appropriate for me to say some- 
thing about why | should have supported the award had I been on the Council 
when it was made. ‘This is perhaps not altogether superfluous, because Council’s 
decisions have not always been unanimously approved. I refer to 1848, when 
the Council, being apparently unable to choose between Adams and Leverrier, 
awarded testimonials to them and to ten other gentlemen, mostly for reasons 
unconnected with the discovery of Neptune. But they did not include Galle, 
who was the first person to see the planet and know what it was ! 

‘We often hear that dynamical astronomy is a dead subject, which means 
that it is a hard subject; most of the easy problems have been done. Further, a 
solution that has fitted for 50 years generally needs revision in 100 years; more 
terms have to be included and the labour goes up like, roughly, the cube of the 
number of terms. ‘The understanding of the physical constitution of the bodies 
of the Solar System needs all the accuracy we can get. 

| first heard of Professor Brouwer when he collaborated with de Sitter on 
the revision of the fundamental constants, especially those related to the Figure 
of the Earth. He also did one of the fundamental papers in the theory of the 
motion of Hyperion. He will be telling you himself about his work on the outer 
planets. His work on the secular accelerations of the Sun and Moon is of the 
first importance, though perhaps he finds it as difficult as I do to see what it 
means. 

“ Professor Brouwer: | present to you this Medal as an indication of our 
esteem for your past work, and accompany it with best wishes for all you may 
undertake in the future.” 
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THE MOTIONS OF THE OUTER PLANETS 


George Darwin Lecture, delivered by Professor Dirk Brouwer on 
1955 April 6 


It has been my good fortune to have been closely associated with three men 
who received the honour of the award of the Gold Medal of the Royal Astro- 
nomical Society: Willem de Sitter, under whom I studied at Leiden; Ernest 
Brown, whose research assistant I became after I left Leiden; and Frank 
Schlesinger, my predecessor as director of the Yale Observatory. ‘I'wo of these 
men were George Darwin Lecturers. It is with a sense of humility that | appear 
before you this afternoon to contribute to the series that has produced so many 
memorable lectures. 

I am expected to speak about my own work, which is like inviting the entire 
astronomical world to scrutinize it and assay it on its merits. ‘This is not to be 
undertaken lightly, but | am encouraged by a piece of advice that Ernest Brown 
gave me long ago when I had been asked to perform a task that up to that time 
had always been the responsibility of more experienced men in the University. 
He merely said: ‘‘ Just do as you are told’’. I have followed this advice many 
times since in similar situations, and I remembered it again as I prepared for this 
occasion, 

During the second half of the nineteenth century the motions of the principal 
planets were the subject of extensive investigations: in France by Leverrier, 
continued by Gaillot; in the United States by Newcomb and Hill. These 
investigations formed the basis of the tables that are in current use for the 
construction of the planetary ephemerides published in the national ephemerides. 
Since these tables were constructed there has been an impressive accumulation 
of observations, generally of higher accuracy than the observational material 
that was available for determining the constants of the tables. ‘This alone is a 
sufficient reason for a resumption of activity in this field at the present time. 

A second reason is the rapid development in the field of high-speed calculators 
in recent years which has opened up possibilities for the study of problems in 
celestial mechanics that were beyond the dreams of our predecessors of even 
one generation ago. ‘This applies both to the handling of large bodies of obser- 
vational data and to theoretical investigations necessary to improve the planetary 
tables. 

Concerning the terrestrial planets the discussions of the observations of the 
Sun by Morgan and Scott (1) and by Spencer Jones (2), and of Mercury by 
Clemence (3) are examples of the very considerable increase in our knowledge 
of these orbits over that contained in the famous summary published by Newcomb 
(4). ‘This work on the terrestrial planets is being continued at the United States 
Naval Observatory; a discussion of the observations of Venus by Mr R. L. 
Duncombe is approaching completion, while Dr Clemence has in hand an 
entirely new theory of the motion of Mars. ‘The first part of this, containing 
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the perturbations of the first order, has been published (5); the remaining portions 
of the theory have since been completed, and the comparison with observations 
is far advanced, 

A new departure in dealing with the orbits of the outer planets was the 
* Coordinates of the Five Outer Planets’’ by Eckert, Brouwer and Clemence 
(6). ‘These orbits were obtained by numerical integration. ‘They were adjusted 
to represent the existing observations and will provide accurate ephemerides 
for a century in the future, except that for Pluto no new comparison with 
observations was made. Pluto’s orbit, based on Bower’s elements XIX (7), 
may require appreciable corrections to represent the observations of the twenty- 
first century. 

In the numerical integration of these orbits the attractions by the terrestrial 
planets were ignored because including them would have required a much smaller 
interval. A supplementary volume by Clemence (8) provides the necessary 
corrections. 

This, briefly, is the present condition of the improvements of the orbits of 
the principal planets. I intend to review a few of the details of the work on the 
outer planets and indicate in what direction we hope to carry on further research 
on this problem. 

The masses of the outer planets.-Yor the numerical integration we were 
content to adopt the same masses for Jupiter, Saturn, Uranus and Neptune as 
had been used by Newcomb and Hill. ‘This had the advantage of making the 
numerical integration orbits more directly comparable with the planetary tables. 
For Pluto we adopted the mass which had been obtained from a 
solution based on Wylie’s residuals of Neptune. 


In the course of further work on the numerical integrations of these planets 
it is our intention to provide coefficients necessary for applying corrections to the 
masses of the five planets. 

A recent investigation by Hertz (9) indicates that the mass of Saturn should 
be increased by one per cent. From a new discussion of the observations of 
Jupiter which extends Hill's solution to the present time, he finds for the reciprocal 
of the mass 3497°64 + 0°27 (p.e.). ‘This brings the value obtained from planetary 
perturbations into excellent agreement with that resulting from satellite orbits. 

For the mass of Neptune, the most comprehensive determination is that by 
Kichelberger and Newton from the observations of Triton (10). ‘They find a 
mass in good agreement with Newcomb’s value, but there is a large difference 
between the result from the visual and the photographic observations. In fact, 
there are good reasons to suspect systematic errors in both. In the years 
1939-1942 Alden (11) made a series of photographic observations of this satellite 
with the Yale telescope, then in Johannesburg. An objective grating was used 
to produce first-order spectral images of the planet of nearly equal magnitude to 
the image of the satellite. ‘This material is now being used by Miss Jocelyn Gill at 
the Yale Observatory, with as one ojbect a new determination of the mass of the 
planet. For this reason special care is given to the scale value. 

It is encouraging that other new work in this field is planned or in progress. 
Dr Daniel Harris has undertaken a study of satellite orbits with the Yerkes 
refractor, and has done a considerable amount of work on the satellites of Uranus. 
An extensive series of photographic observations of Saturn’s satellites made 
with the Yale telescope in the southern hemisphere was measured at the United 
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States Naval Observatory by Mr F. H. Hollander. ‘The reductions have recently 
been taken up by Dr Morris S. Davis, a former student of mine, now at the 
University of North Carolina. One of the principal objects of this study is the 
determination of Saturn’s mass. For the mass of Jupiter Dr Clemence has 
undertaken the preliminary steps toward the realization of a proposal made by 
Hill (12) in 1873. Hill selected 13 minor planets near the 2: 1 resonance region 
with Jupiter for which the principal Jupiter perturbations in longitude have 
coefficients between 1°-4 and g’, and periods between 60 and 120 years. These 
planets are being observed at the United States Naval Observatory; the old 
observations have been collected and are being reduced to a uniform system. 
Dr Clemence expects that with the aid of electronic calculators it will soon be 
possible to calculate orbits of the requisite precision. 

In the Introduction of the volume of the “ Coordinates of the Five Outer 
Planets ”’ a solution is made for the mass of Pluto from its action on both Uranus 
and Neptune, based upon the residuals obtained from Wylie’s discussion (13). 
Since then Woolard (14) has made a new compilation of the pre-discovery 
observations of Uranus and compared them with the orbit obtained by numerical 
integration. I have recently made a new solution for Uranus in which the pre- 
discovery observations are included. ‘The accidental error of these old observations 
is disappointingly large. ‘The r.m.s. residuals in longitude and latitude are 
sg and 4"1, respectively. The solution depends to a large extent on the 
residuals of Flamsteed’s observations, one in 1712, three in 1715. ‘These are not 
accordant. Rejection of the observations in 1715, as is favoured by Kourganoff 
(15), yields a mass of Pluto comparable with or greater than the Earth's mass. 
It appears to me that the only safe procedure is to give all pre-discovery observa- 
tions equal weight and to refrain from rejecting any. Weights 0-05 and ov1o in 
Wylie’s system were assigned to the longitude and latitude equations, respectively. 
‘The result of the solution, if is the fractional correction to the mass of Pluto 
used in the integration, is 


0°58 + (m.e.). 
A new solution for Neptune from four-year means gives 
—0°06 + 0°14 (m.e.). 


‘This mean error is 60 per cent greater than in the published solution. ‘l'his 
must be ascribed to the use of four-year means rather than the individual normals 
given by Wylie, an indication of lack of independence of the normals. A weighted 
mean of these solutions is: 


p= + 0°13 (m.c.) 


My = + 0°10 (p.e.)}. 

‘The Neptune residuals look very satisfactory ; the longitude solution depends 
entirely on the Lalande residual of 1795, but the modern latitudes cannot be 
satisfied with a much smaller mass. ‘The Uranus residuals are less sensitive to 
the mass of Pluto. ‘The latitudes contribute practically no weight at all to the 
determination, but, regardless of the adopted mass of Pluto, they show a 
systematic trend. 

Essentially the same systematic trend in the residuals in latitude was found by 
Newcomb when he compared the observations with his first theory of Uranus (16). 
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Newcomb represented the residuals by an expression which becomes 


+0°18—0"28T 
if 7 is expressed in centuries from 1g00. He concluded that at the epoch 1840 
Auwers’ equatorial declinations are too small by 0”-35 and that this error is 
diminishing at the rate of o”-28 per century. This conclusion seemed reasonable 


enough in 1870. Wylie obtained the result 
+0°13-0"21T 
from the observations since 1830, compared with Newcomb’s second tables. 
He states that the declination system of the FK3, to which he reduced the obser- 
vations, is not likely to be affected by a systematic error of this size. ‘The obser- 
vations since 1870 compared with the numerical integration orbit give : 
This persistent unexplained deviation in the motion of Uranus destroys much 
of the confidence in the solution of Pluto’s mass from Neptune. For if the 
latitudes of Uranus show the presence of an unknown systematic effect, how can 
we reasonably expect that the latitudes of Neptune are free from a similar unknown 
effect that is absorbed in the solution for the inclination and node and the mass 
of Pluto ? 

If we assume, against better judgement, that the systematic trend in the 
Uranus residuals indicates an error in the FK3 declinations, we must correct 
the Neptune residuals in latitude by subtracting the linear term found from 
Uranus, If this is done, the solution from the latitudes of Neptune gives 
jo +0°39, increasing the mass of Pluto. 

Kuiper (17) objects to a mass of Pluto comparable to the mass of Venus on 
the grounds of Pluto’s apparent faintness and small measured diameter. Until 
recently, judging by the evidence from Neptune alone, it appeared to me that it 
would be difficult to bring this physical estimate into agreement with the gravita- 
tional determination. I have now come to the conclusion that the gravitational 
determination is not entitled to the weight that the formal solution indicates. 
The longitudes of Uranus and Neptune combined give 

—0°20 + 0°24 (m.e.), 


+ 0°20 (p.e.)}. 


Much of the weight depends on the observations of Neptune by Lalande in 1795. 
But Jackson (18) showed in 1930 that a wide range in the elements of Neptune 
is possible without requiring the presence of improbably large systematic errors 
in the observations. For a slowly moving object such as Neptune these systematic 
errors are in part absorbed in the solution and the resulting weight 1s fictitiously 
high. ‘Thus, after a great deal of work has been done on the problem, we must 
for the present leave the question in this state of uncertainty. As Kourganotf 
remarks, the perturbations by Pluto in the longitude of Uranus will soon increase. 
By 1965 the difference in longitude between mp=o and mp= qq will 
amount to 0:3 and the weight of the solution from the Uranus longitudes will 
have increased fourfold compared with the weight of the solution that included 
observations to 1938. Similarly, the perturbations on Neptune must eventually 
increase so that the solution from the longitudes of Neptune will become less 
dependent on the Lalande residual of 1795. 
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Only with the passage of time may we expect a reliable gravitational deter- 
mination of Pluto’s mass and a clue to the cause of the puzzling deviation of 
Uranus in latitude. ‘The latter may well be of an exciting nature. 

Accuracy of the numerical integration orbits,—As a test of general interest 
Clemence and I (19) evaluated from the coordinates and velocity components of 
the numerical integration orbits the three integrals of areas and the energy integral 
for five dates, two of which are very near the beginning and the end of the 
integration, while one coincides with the osculation date. ‘The differences 
from constancy are on the whole in accordance with expectation on the basis 
of accumulation of error of rounding, though in the case of one of the integrals 
of areas the discordance among the individual values is approximately six times 
that expected. We tentatively ascribe this to truncation errors which may arise 
in Jupiter and Saturn. ‘The effect is of little practical interest since it does not 
reach the last printed decimal place in the coordinates. Yet some of the details 
have at least the interest of a curiosity, As stated in the Introduction to the 
volume, the tenth differences of the accelerations of Jupiter and Saturn show 
from time to time systematic trends resembling harmonic functions of relatively 
short period. Further examination showed that they occur in Jupiter near each 
perihelion, in Saturn near each conjunction with Jupiter, recurring in cycles of 
about 60 years. An unforeseen circumstance is that the arbitrarily chosen 
osculation date, 1941 January 6-0 = JD 2430000°5 is very near a Jupiter peri- 
helion and a conjunction of Jupiter and Saturn. It would have been difficult to 
choose a less suitable date for overlapping the forward and backward integrations. 

All the available evidence indicates that for the interval of four centuries for 
which the numerical integrations were carried out, these orbits are of a quality 
that would be extremely difficult to attain with general theories. ‘There remains 
the objection that in the past we have become accustomed to providing planetary 
ephemerides on the basis of general planetary theories. 

It has often been said that the drawback of the numerical integration 
procedure is that the orbit is determined only for the period for which the inte- 
gration has been carried out, and that no information is obtained concerning the 
general properties of the orbit. In this connection I refer to a recent paper by 
Stiefel (20), in which the relation between numerical integrations and general 
planetary theories is considered. Dealing with a system of the Sun and only 
two planets, he shows that it is, at least in principle, possible to derive from the 
result of the numerical integration of the orbits a representation that corresponds 
to that of a general planetary theory in the ordinary form. Perhaps with powerful 
electronic computing facilities a method of practical value could be developed. 
Even then an extension to more than two planets would lead to formidable 
complications. 

Whatever the merits of such an approach may be, it would indeed be 
regrettable if in the future purely numerical methods were to supplant analytical 
methods completely. ‘he numerical methods should rather be looked upon as 
an aid for testing the results of analytical theories and for stimulating further 
theoretical developments. ‘This sentiment was expressed by Danjon (21), and 
I concur wholeheartedly, I shall now turn to the subject of general planetary 
theories. 

Classical planetary theory in rectangular coordinates.—\n the introductory 
section of Hill’s ‘‘ Researches in the Lunar Theory " (22), the author remarks that 
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in elliptic motion the rectangular coordinates can be developed in periodic series 
in the mean anomaly each coefficient of which can be stated as a single Bessel 
function. The law of the series is evident, and the approximation can be carried 
as far as we wish. ‘The longitude and the latitude have much more complicated 
expressions ; the coefficients in the series cannot be expressed by a finite number 
of Bessel functions. Hill then continues: “ If this is true in the elliptic theory, 
how much more likely is a similar thing to be true when the complexity of the 
problem is increased by the consideration of disturbing forces?” 

lor numerical integration Encke had introduced the differential method in 
rectangular coordinates, being attracted by the simplicity and symmetry of the 
equations. He subsequently made an attempt to apply the same principles to 
the problem of general perturbations (23), but in the form developed by him the 
perturbation in the radius vector appears as a fourth quantity to be obtained 
in the course of the integration of the equations. ‘The introduction of two 
superfluous constants of integration is a severe drawback. It is curious to note 
that both Briinnow, in an appendix to Encke’s article, and Hill (24) use Encke’s 
equations as stepping stones toward the derivation of the familiar equations in 
polar coordinates. 

There is an indication of a promising approach to obtaining suitable equations 
in rectangular coordinates in Dziobek’s treatise (25), but Dziobek was satisfied 
with a final result in symbolic form which he described as unsuitable for practical 
tise, 

‘The principle is the following: For convenience choose an elliptic reference 
orbit nearly coinciding with the orbital plane of the perturbed body so that the 
s-coordinate is of the order of the perturbations. Let x», yg designate the 
coordinates in the reference orbit as functions of the time and the chosen fixed 
elements, Let the perturbed coordinates be x, +£, vo +7. ‘The equations for 
€ and » are 


Vy 


dt? + Yon) * X= + AX, 


d*y | 3H Ve OR 


The right-hand members consist of the x and y components of the disturbing 
force computed with the coordinates of the reference orbit with the additions 
AX, AY, which consist of squares and products of € and » but may also contain 
higher-order parts of the disturbing forces. 
If the right-hand members are ignored, the equations have the particular 

solutions 

OVe 

}=1,2,354) 


in which ¢, are any set of four orbital elements that define the motion in the plane 
of an elliptic orbit. ‘Thus the complete solution of the equations without right- 
hand members is 

joe * 


where (, are four arbitrary constants. By the method of the variation of arbitrary 
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constants the same form of solution can be used if C, are obtained by the inte- 
gration of dC, dt from the equations 


y 

“dc, dt 
In a sense this completes the solution, but without further modification the 
method is hopelessly laborious and has a number of practical disadvantages. 
The principal drawback is that the resulting expressions require numerous 
multiplications of series, both in the integrands and outside the integrals. ‘This 
is overcome (26) by using the device due to Hansen of bringing a factor outside 
an integral sign into the integrand by treating the time it contains as constant 
in the integration. ‘This requires this ¢ to be kept distinct from the ¢ already 
present in the integrand, which is accomplished by introducing a different symbol, 
say, 7 for the ¢ so arising. After integration 7 is to be restored tot. Functions 
in the integrand to be expressed in terms of r are written with a bar over them, 

The final expressions for the perturbations £ and » become 

+ |(4,X+A,Y)dt, Ye \(B,X+B,Y)dt, (2) 


in which 


| (ez + ) aear 


a». We 
The functions 


A ( A, (Yo» Xo), Bb, Vo), B, (Vo Vo) (4) 
are written as Poisson brackets because this is the form in which they are expressed 
in terms of partial derivatives of x,, y, with respect to a set of elliptic canonical 
variables. ‘They are of the same nature as A and B which arise in Hansen's 
method in the function W from which the Hansen perturbations ndz and v are 
obtained. It is unnecessary to deal with the perturbations perpendicular to 
the plane of the reference orbit. It is evident that the expressions (2) and (3) 
can at once be generalized for any arbitrary fixed coordinate system. With the 
particular choice of the xy-plane the equation for 2 reduces to a form common 
to that in other types of planetary theory. 

Some similarity with details of Hansen's method is evident; yet the final 
expressions are totally different. A comparison between the two methods 
involves much more than a mere counting of the number of operations required. 
An actual application of the method was made by Davis (27), who also made a 
parallel computation by Hansen's method. A similar comparison was made 
by Clemence (28), who also comments in part in the Celestial Mechanics section 
of the Reports on the progress of astronomy in 1951 (29): “ The calculation of 
the perturbations of higher order than the first is greatly simplified. In addition 
to this important advantage, the developments possess an unusual degree of 
symmetry. ‘This symmetry is obtained at some sacrifice of brevity in the 
numerical development of these factors, which renders Hansen's method the 
readier one if it is proposed to stop with the first approximation, and which to 
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some extent counterbalances the advantages of the new method in the second 
approximation. ‘There is no doubt, however, that in the third and higher approxi- 
mations the new method is preferable to all others. Its use would be indicated 
if new theories of the outer planets were to be constructed and, in most cases 
where it is applicable, if it were desired to obtain a general theory of a minor 
planet comparable in precision with the observations.”’ 

This, | believe, is sufficient justification for the use of the method in rectangular 
coordinates for the development of new theories of the planets Uranus, Neptune, 
and Pluto in which the perturbations of higher order are important. ‘This work 
is Now in progress at the Yale Observatory as a joint undertaking by my colleague, 
Dr van Woerkom, and myself. 

For the development of the odd negative powers of the distance, the method 
given a few years ago by Clemence and myself (30) proved useful. It assumes 
that an approximate development of the reciprocal of the distance is available. 
Denote this by \'. Let the square of the distance be D*, which can be computed 
as a double Fourier series with all the accuracy desired. ‘The square \ * and the 
product D?\"*=1+.6 are then obtained. ‘lhe deviation « from unity shows 
the error in \''. It is then easy to obtain 

= +e? — +...), 

Dr = 1 — be + fe? — +...), 
and from these by multiplication of series all the odd negative powers of D 
required, With adequate automatic computing equipment the method is 
probably not more elaborate than the conventional development. All of the 
operations can be performed in strictly routine fashion, and the operator has to 
learn only one unfamiliar technique, the multiplication of series. ‘The simplicity 
and certainty of the method are its principal advantages. 

For applications to the orbits of the principal planets provisional developments 
of the reciprocals of the distances in most cases are available. ‘This is not so 
for the distances between Pluto and the four major planets. For these, 
provisional developments were obtained by harmonic analysis. Only the 
distance Neptune-Pluto gave serious difficulties, but harmonic analysis from 
So special values gave an expansion that could serve as a first approximation. 

The work now has the highest priority on our computing programme, It is 
a task of which my colleague Dr van Woerkom has taken an important share. 
In fact, all of the machine calculations are performed either by him or under his 
immediate supervision. 

Generalized planetary theory.—Betore the decision was reached to concentrate 
on this problem, a considerable amount of experimentation on a more ambitious 
project had been carried out. ‘The usual form adopted for general planetary 
theories in which powers of the time are permitted in the coefficients of periodic 

terms is imposed on us on account of the difficulty of proceeding otherwise. 
Hill once called this form a “temporary makeshift’. A much more satisfying 
solution would be one like that of the main problem of the lunar theory in which 
linear combinations of the four fundamental arguments, D, /, /’, F appear in the 
arguments of the periodic terms of the Moon’s longitude, parallax and latitude, 
and no powers of ¢ appear in the coefficients. In the satellite systems the mere 
fact that the motions of the pericentres and nodes proceed so much more rapidly 
than in the planetary system compels us to adopt this form. 


pe 
tee 
aa 
fl 2 
AG 
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The literature of the second half of the nineteenth century is rich in investi- 
gations on this form of solution for the planetary problem. Most of these were 
qualitative in character, but Gyldén and Hill were especially concerned with 
devising methods that could be used to obtain a representation of the orbits of 
the principal planets. In the fourth volume of his ‘‘ Collected Mathematical 
Works’, Hill returns to the problem again and again. He devoted particular 
attention to the application of Delaunay’s method to the planetary problem. 
It is perhaps curious that Hill never conceived the modification of Delaunay’s 
method in which the non-critical terms in the disturbing function are treated 
together in one operation. This approach was elaborated in much detail by 
v. Zeipel (31) and also treated by Brown and Shook (32). Unfortunately, the 
procedure requires that all the parameters involved are left indeterminate in the 
formulas. As a consequence, the degree of complexity is forbidding. Thus, 
while the problem has been foremost in the minds of some of the outstanding 
men in the field, we still lack a satisfactory representation of planetary motions 
valid for long intervals of time. ‘This is not to be ascribed to lack of competence 
or perseverance of these investigators but merely to the magnitude of the problem. 

The only completed projects along this line are solutions of the equations of 
the secular variations by Leverrier, Stockwell, Harzer and recently with improved 
values of the planetary masses by van Woerkom and myself (33). ‘The scope of 
this solution is extremely limited compared with that of providing fully general 
theories of planetary motion. Yet it has at least the merit of indicating the form 
of a general analytical perturbation theory in which the time is present only in the 
arguments of sine and cosine terms. 

It seems worth while to explore the possibility of using the solution of the 
secular variations as the basic data for a more general attack on the problem. 
‘This indicates a representation in which the arguments of the trigonometric 
functions appear in a literal form, while the coefficients are to be treated 
numerically. An important decision to be made is whether to aim at providing 
for all the principal planets at once, or begin with a more modest undertaking. 

There is a good reason for concentrating first on the outer planets because, 
though they produce large secular variations in the clements of the terrestrial 
planets, the effects of the terrestrial planets on the outer planets are insignificant, 
It would not be possible at present to extend the solution to the planets Uranus, 
Neptune and Pluto. The orbits of Neptune and Pluto are such that if the 
longitudes of the perihelia and nodes of these planets are permitted to vary without 
restriction, the orbits will at some time intersect. ‘his introduces a singularity 
into the problem that makes it impossible to develop the disturbing function im 
the general form that is required. For this reason, van Woerkom and | were 
obliged to leave Pluto out of the solution of the secular variations. ‘This does 
not necessarily indicate that no solution of the desired form exists. ‘The three 
outermost planets have mean motions so nearly in the ratio 6; 3:2 that it is not 
excluded that their orbits are actually connected by a resonance relationship. 


The completion of the ordinary planetary theory for these planets may help to 


decide whether this conjecture has any merit. At the present time, therefore, 
we have no choice but to begin with the system consisting of the Sun, Jupiter and 
Saturn only. Further, to simplify the problem, we will ignore the inclinations. 
This surely is the simplest case that, if solved, could serve as the nucleus of a 
more ambitious project. 
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Let A be the mean longitude, a the semi-major axis, ¢ the eccentricity and 
w the longitude of the perihelion, # the true longitude, r the radius vector, and put 


raali+p). 


Then 
if ve a(i+p)cos(A+&), y=ali+p)sin(A+e). 

1+f&=(14 p)cosé, == (1 +p)sine, 
then 


x= a(t + €)cosA~ansind, 
y= + €) sin A + an cosa. 

In these expressions, a and A appear explicitly; € and 7 are expressible as 
cosine or sine series the arguments of which are multiples of the mean anomaly, 
/=A-—wm, while the coefficients are power series in the eccentricity. If these 
expressions apply to Jupiter, similar expressions with primed quantities could 
be written for Saturn. 

In the usual form of planetary theory the disturbing function is developed 
by using constants for a, ¢, » and for A a linear function of the time. In a first 
development of a generalized theory it is necessary in the first place to introduce 
for ¢ and w the expressions that are obtained in the solution of the secular 
variations. In addition it 1s possible to incorporate the principal perturbations 
in the mean longitude, transformed to the desired form. 

‘The principal terms in Hill’s solution for the eccentricities and perihelia of 
Jupiter and Saturn (34) are 


sin sin 68 sin, 
cos 9°043749 cos * * cos *’ 
sin sin, 
© cog 0048208 X's | (5) 


x= + + const., 

x = +277°73017 + const., 
/ being expressed in years. ‘The designations of y and y’ have been interchanged 
compared with Hill. Substituted into the expressions for & and », if A, and A,’ 
are linear functions of ¢, the results are: 


Jupiter 
é 7 
Argument cosine series sine series 
Ao 0°04379 0°08 746 
Ao o’ors7i 0°O31 34 
2Aq 2x 0°00095 000045 
2Ag 0000% 0°00034 
2A, 000010 000005 
x x 0°00034 
2x 2x ©'0000 1 


Saturn 


Argure rit cosine series sine series 
Ao \ 0°032451 0°06540 
Ao x 0°04535 0°09645 

2Ay 2x 0°00053 000026 
\ x ©°00157 >*0007% 
2x 0°00122 000061 


0°0008 3 
0°0000}3 


ii 
; 
| 
x 
| 
7 
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The expressions (5) are then used to transform Leverrier’s results for the 
great inequalities in the mean longitudes of Jupiter and Saturn. They are: 


Long-period terms in the mean longitude in radians 


Argument (sine series) 


Jupiter Saturn 
2Ag—5Ag — x 0°00000 0°0000 | 
2Ag—5Ag 2x + 5x" 000005 000013 
2Ag— 5Ag 3X 000062 0°00153 
— 5x + 2x" 0°00167 0°00411 
000011 0'00027 

7 00001 0°00002 


It would have been possible to transform all the larger terms of Leverrier’s theory 
into the required form, but it seemed doubtful that much could be gained in this 
way in the first approximation. ‘The only additional terms actually used are 
the synodic perturbations, which are comparable in magnitude to the quadratic 
terms in the equation of the centre. It should be noted that the principal 
coefficients of the great inequalities in the mean longitudes are intermediate in 
magnitude between the coefficients in the equation of the centre that correspond 
to the first and the second powers in the eccentricities. Hence, in a numerical 
development the inclusion of these terms is bound to accelerate the approximations. 
The principal task in the development of the disturbing forces is the expansion 
of the odd negative powers of the mutual distance. In the particular problem 
before us, the inverse cube of the distance, or rather a““A *, is the basic function, 
a’ being a constant nearly equal to the mean value of Saturn’s semi-major axis. 
The computation of a’ *A* as a cosine series with arguments 


x) x’) + x) 
is readily obtained. ‘This may be written 
a’ 2A? = A,* +O, 


where 
A,? = 1 — 2a cos (Ay — Ag’), 


and 


%=O'545 199. 
rhe principal terms in QO are 

+ 0°13209 cos (Ay’ — x’) 0°08816 C08 (Ay x’) 

+ 002751 COS (Ay — x) + 0°02032 COS (Ag — 2A,’ + x’). 


rhe expansion of a’A~* in powers of O by the binomial expansion 

I 30 15 

2457 
would have led to an exceedingly long process. ‘This was shortened by the 
following device. From 


I 
+ cos p(Ag — Ag’ )} 


we obtain by multiplication of series 
K 1 


Ay 20. 


= 

4 
. 
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Phen 
a’ = 1 + K,). 
Let the principal periodic term in K, be —2a,cos¢,, where «, may be positive 
or negative. Put 


=AP+Q, 
*) 
+ Ky). 


Next, the largest periodic term in K, is selected and the process is continued 
Thus, finally, 


= 


Each of the factors in the product A, *.A,*.... is readily expressed with 
the aid of the well-known Laplacian coefficients. Only the series for Ay * is 
slowly convergent, the remaining factors consist of few terms. In principle, 
this procedure could be continued indefinitely by including in the successive 
\'s smaller and smaller terms. Actually we stopped with A, and the expansion 
of (1 + K,y** was obtained by the binomial theorem. In A, there were three 
terms with coefficients exceeding +0°02. ‘The fifth power of A, was sufficient 
to assure the accuracy required, Inthis manner the development of the derivatives 
of the disturbing function consists of multiplications of harmonic series, all with 
numerical coefficients. Actually exponentials were used by introducing 


u=E+in, s=&—1y. 


Except for the preparatory stages all the calculations were carried out at the 
Watson Scientific Computing Laboratory. I am indebted to Dr W. J. Eckert 
tor making available the facilities of the Laboratory and to Miss R. B. Jones of 
his staff for supervizing the calculations. 

An indication of magnitude of the operations may be obtained from the 
following facts: the basic series containing the equations of the centre and the 
perturbations had been limited to 25 terms for Jupiter, 33 terms for Saturn, 
rounded to five places of decimals. ‘The development of a’*A-* was carried to 
seven places of decimals. ‘The numbers of different arguments of cosine terms 
present in the principal functions are : 


Ag *.4,-*. 248 arguments 
(1+K,)*? 1296 arguments 
2198 arguments. 


Hill’s series for a“°A~* has 183 arguments, but in his development both sines 
and cosines are present. ‘The difference is, of course, that in the ordinary general 
planetary theory the arguments are linear combinations of the two mean anomalies. 
In the more general development in this case, a third argument, x’ — x, 1s present. 
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Thus, in a“A-*, Hill’s terms 

+0°00912144 cos (2/—5l') 0°48984397 cos (/ 

+0°01926179 sin (2/—s5l’) +0°§2282903 sin (/ 
correspond to: 

+ @°0000010 cos (2A, i x) 0°0000188 cos (Ag 
0°0000002 cos (2A, ) + 0°0018606 cos (A, 
0°0004476 cos (2A, x) +-0°0589256 cos (A, 
+ 00062662 cos (2A, “+ y+2y’) 0°7153388 cos (A, 
0°0214232 cos (2A,—: 3x) 0°0294726 cos (A, Xx 
0°0025462 cos (2A, x+4x) 0°0014666 cos (Ay 2x 
-0'0002788 cos (2Ay~- 2x + §x’) 00000084 cos (Ay 3x 
+ 0°0000154 cos (2A,q + 6y’) 
These samples are representative ; for most of Hill’s arguments we find seven or 
eight different arguments. ‘This same multiplication of arguments is found in 
Leverrier’s theory, but Leverrier’s longitudes of the perihelia have the instan- 
taneous values. 

That the ratio 2198 183 1s considerably greater than 7 or 8 must be ascribed 
to our retention of all coefficients of cosine terms numerically equal to or greater 
than 2x10. Many of the very small coefficients are associated with high 
multiples of A, and A,’, and could be safely ignored. As to the number of decimal 
places, Hill begins with eight for the lower multiples, which number is reduced 
for the higher multiples to seven, six and finally to five. There is a further 
advantage on Hill’s side: the principal long-period terms in the mean longitudes 
are obtained in his theory from the development of a’\'', for which he uses 
eight decimal places in most of the coefficients, but ten for the argument 2/— 5/’. 
In the more general theory, use could be made of a’A ' at the cost of some further 
complications in the development. For the present the disturbing forces have 
been obtained with the development of a*A* only. It would therefore be 
necessary to carry further developments to additional decimal places if numerical 
accuracy comparable to Hill’s tables were desired. For the immediate future 
1 will consider our goal reached if it turns out that the procedure converges 
satisfactorily, even if the result cannot compete in the number of significant 
figures with an ordinary planetary theory. 

So far I have dealt only with the development of the disturbing forces. What 
method of integration is to be used ? 

There are tremendous advantages in the use of a method in which the 
coordinates are obtained directly by the integration, especially since it will be 
necessary to substitute the results of the first integration into the developments 
for obtaining a second approximation. By using coordinate systems rotating 
uniformly with the mean mean-motions of the planets, equations in rectangular 
coordinates are obtained which resemble closely those used by Hill for the lunar 
theory. ‘The one term in the differential equations of the lunar theory that, as 
Hill remarks, is a hindrance to their ready integration is the term having the 
factor r * not multiplied by the disturbing mass. In the lunar theory this 
singularity at zero is eliminated with the aid of the Jacobian integral. In the 
planetary problem with two planets there are terms not multiplied by the perturb- 
ing mass with factors r-* and r’~* to contend with. ‘The energy integral contains 
a linear combination of rand r’~'. Obviously there is no possibility of accom- 
plishing the elimination of the singularities at zero, After much experimentation 
the use of these equations was given up. 


5 
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A promising alternative approach is to express the derivatives of the elliptic 
clements in terms of the partial derivatives of the disturbing function with respect 
tothe coordinates. [am now reconciled to the necessity of obtaining perturbations 
in the elements and having to perform the intermediate operations to secure the 
coordinates which must be substituted in the disturbing forces for a new approxi- 
mation. This changing from elements to coordinates adds considerably to the 
number of multiplications of series that is required. ; 

The crucial question is whether the solution of equations for the eccentricities 
and perihelia will give unforeseen difficulties. ‘The indications from the work 
that has so far been done on the problem are favourable. ‘The equations for the 
semi-major axes and mean longitudes are straightforward. ‘The principal 
difficulties are those of the small divisors, present in any form of planetary theory. 
In the Jupiter-Saturn problem the small divisors associated with low multiples 
of the mean longitudes cause a loss in the number of decimal places after integra- 
tion. In another respect these are helpful since they prevent the occurrence of 
other small divisors, except those associated with such very high multiples of 
the mean longitudes that they are of no astronomical significance. Of course, 
they are of importance in connection with questions concerning strictly mathe- 
matical convergence. 

When the developments of the inverse cube of the distance were made, 
about six years ago, the computations taxed the capacity of the calculators then 
available. In the meantime new calculators have been developed that are much 
more efficient for this type of computation. I am now encouraged to proceed 
with further work on this problem. 

| may seem bold to devote so much of this lecture to a project that 1s stll far 


from completed. ‘The only excuse I can offer is that, unless I had limited myself 
to generalities, there seemed to be no other way to explain the possibilities of 


research in gravitational astronomy with modern high-speed calculators. 

I have tried to stress the importance of team work in all these undertakings. 
About eight years ago the Office of Naval Research in Washington decided to 
lend its support to a research project on the motions of the principal planets, 
which was then formulated as a joint undertaking of the United States Naval 
Observatory, the Watson Scientific Computing Laboratory and the Yale 
University Observatory. If we have been able to accomplish some worthwhile 
things, it has been possible only with the magnificent support that we have 
received and with the unselfish cooperation of everyone connected with the 
work. With continued support and continued team work | am hopeful that we 
will succeed in completing the unfinished parts of the project. 
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RR TELESCOPIL IN 1953 AND 1954, AND THE DEVELOPMENT 
OF SLOW NOVAE 


A. D. Thackeray 


(Received 1955 February to) 


Summary 


Recent developments in the emission spectrum of RR ‘Tel include the 
emergence of C tv, [Atv], [K tv}, [Ferv], [Fev]. A special search for the 
hitherto unknown [Fe tv| lines proved successful in the light of Edlén’s 
independent predictions from laboratory data. Fe 11 lines have faded steadily 
relative to [Fe 11] 

The rate of development towards higher stages of ionization is compared 
with that of other slow novae. ‘Time intervals between corresponding stages 
of RR Pic, RR Tel, RT Ser and » Car are found to be approximately in the 
ratios 1:6:17:220. One result of this comparison is to give some grounds 
for supposing that 7 Car brightened from a pre-nova state not very long before 
Halley's observation in 1677 


1. Radcliffe spectra of the interesting slow nova RR Tel, covering the period 
1949-1952, have been described in two previous communications.* It will be 
recalled that the star, which had been known since 1908 to be variable 
(125 to 16™), rose sharply to 7th magnitude in late 1944, remaining fairly 
constant until 1949, since when a gradual decline to about 9th magnitude has 
set in. In 1949 June to August a fairly rapid transition was recorded on low- 
dispersion Radcliffe spectra (Paper 1) from absorption (H, Till) to emission 
(H, Fett). ‘The much more detailed Cassegrain spectra of 1951-52 showed 
an emission spectrum of advanced Eta Carinae type (Paper 11); between 1951 
August and 1952 May there was a marked strengthening of broad bands due to 
[Ort] and [Nett]. 

The present paper is concerned with identifications of lines in spectra of 
RR Tel obtained in 1953 and 1954 with the two-prism Cassegrain spectrograph 
attached to the Radcliffe reflector. In addition, one useful spectrum (120 A/mm) 
extending beyond the Balmer limit to 3300A was secured in 1954 with the 
Newtonian spectrograph, which, equipped with an aluminized plane grating, 
is now freed from all glass absorption, 

Table I lists the observational material secured during the period. Not all 
the plates have been measured, but if the plate was used in the present study 
the last column of the table shows the region covered by the measures. Dr D. 
Koelbloed of the Amsterdam Institute obtained some further material in 1954 
for the purpose of a spectrophotometric study which will give valuable quantitative 
information concerning the physical state of the object. One of his plates 


* A.D. Thackeray, M.N., 110, 45, 1950; 113, 211, 1953 
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(taken for the infra-red region) is included in Table I, and the writes is deeply 
indebted to Dr Koelbloed for the use of this plate. lhe other plates in the table, 
and all the measures, are by the writer. 


Tasie | 
Observations of RR Tel 195% 54 


Exposure Region 


(mun) measured (A) 


Plate Date Emulsion 


1341 March 20 1ro3zaQ) 
€ 1342 March 20 103aF 
€ 1457 May 13 103aF 4630-0793 
1547 June 24 1o3aQ) 10 3869-5042 
> 1§79 2 rojal 270 $930--7005 
1599 ¥y 20 IN 2 6563-8750 
1649 1o3al 
1663 1ozaQ) 
1719 Sept. 1o3al 26-0793 
> March 2! 
19s April 5 1o3aB 29-6371 
> 1955 May 10 103aB 371 
> 2023 June 16 10o3aB 
2023 June 16 
> 2023 June 16 
2022 June 16 
2040 July 3 
2040 July 3 
2041 July 3 103 | 
b 2106 July 26 1ozals 
2165" Sept. 1073 
2196 Oct. ro3zals 
'T 21 Oct. 24 R « 
* Observation by D. Koelbloed 


t Newtonian pectrogran 


2. Since the majority of the 250 lines observed in 1953-54 appear already in 
Table V of Paper II, no useful purpose is served by giving anew a complete 
table of wave-lengths and identifications. ‘lable II therefor: presents a selected 
list of measured lines, which omits, for instance, the large number of tamuliar 
lines due to [Fe u] and Fen. ‘The criteria for inclusion in ‘lable I have heen: 

(i) Unidentified lines. 

(11) Lines not observed prior to 1953 

(it) Lines which have clearly strengthened since 1952 

(iv) A few representative lines, previously observed, to serve as an indication 
of the intensity-scale; e.g. the Balmer lines in the ultra violet, [Ou] 7324, 
tt} 4959, ete. 

‘The arrangement of the table 1s as in Paper I1, an asterisk implying a footnote 
at the end of the table. ‘The identifications are as before based primarily on 
Mrs Moore-Sitterly’s Revised Multiplet Table. \n addition, the writer is deeply 
indebted to Professor B. Edlén for communicating, prior to publication, information 
concerning wave-lengths of permitted lines due to C, N, O in various stages of 
ionization, and also to Dr R. Velasco for similar communication concerning 
Nim. Professor Edlén has also very kindly communicated information 
about [le tv] (see below). 
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Selected lines in RR Tel, 1953 and 1954 


Be Wave- 1953 
x : length (A) Int n Identification 
3343 i 42°9 
3475 78°7 
4490: $3°0, 
3032" 
47° 
gi’ 2 
3697 2 
3703" 2 
2 
14°? 
21° 2 
2 
I 
I 
2 
2 
1 
2 
2 
2 
I 
2 
2 
1 
2 
2 
2 Feu 33°17 
2 Feu 93, 44°51 
2 Nu 
KK Iv ‘0 
2 ?Mn tv 3 
Nu 16 
Nu 6, 41°9 
(blend) 
1 
4 He 
3 


?Fe 
Ow 


2 

3 

3 
3 

3 
3 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 


aF 
2F 96°79 
Vv 3F 20°2 
Vv iF 38-1 
58-07 
12 13°2 
4868-2: 4: 8: Fe 
4900°2 2 Fe tv 
06°2 6 12 Fe tv F 
30°'8 in 2n 430°5 
4955°79 40 50° if 
§033'27 10 
41°55 8 15 Sil 5 
§6°2 3 6 Siu 5 5671 
5060°3 4 Fe iF 603 
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Wave- 1953 
length (A) Int Identification 


716 
1F 848 


5071" 
5054 
5147" 
3F 
5215" 

33° 
53095 
5323 
5412° 
5494 
5515 

35 

44 

52° 


67° 


= 


77° 
§592° 
5666 
5679 
5753" 
5799 
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4490 On the Newtonian spectrogram (120 A/mm) this strong very diffuse feature is barely 
resolved from 4478 and the measured wave-length is very uncertain 

3721'S ‘The contribution of [S 111] to this line is indicated by comparing its intensity with 

those of nexwhbouring Balmer lines 

3891 °65 (Plate 2(a)). Wave-length and intensity uncertain owing to proximity of H 3889. 

3954°4 Identification very uncertain 

4997°3 Strengthened from 1953 to 1954; well resolved from H4(Plate 2(4)). Not 


recorded in Paper II but, like N 111 4640, was probably present, very diffuse, 
blended with the wing of Hd 

4152°46 (Plate 2(h) Unidentified. Observed also in 7 Car. ‘The ascription to [Co 1] 

in Paper II is now ruled out by the subsequent behaviour in RR Tel 

4206°62\ Unidentified. Observed in RR Pic as a single line 4207. Assuming the two 

4209°04f lines to have a common origin the separation may provide a clue to the 
identification 

4571°24 Can searcely be due to Mg | at this stage of ionization 

4368°2 (Plate 2(c¢)). Wave-length and intensity uncertain owing to proximity of Hf. 
Identification with [Fe tv] due to Edlén (private communication). 

4930" Probably due to Fein (1F). [Ot] 4931°8 should be several thousand times 
fainter than the nebulium pair 

§233°9 Fe ti must be a minor contributor. ‘l'his unidentified line appeared in RR Pic. 

§494'1 Mainly unidentified ; this line is too strong for Fe 17 F, and the wave-length 

disayrees 


§679°S ‘The absence of 5177 shows that the contribution of [Fe v1] must be small. 
sgit Observed in N Her 

6312°2 As in Paper II the wave-length disagrees from the prediction for [S 111] but there 
can be little doubt of the correctness of the identification 


6793'2 ‘The lower intensity recorded in 1954 is a result of the insensitivity of the emulsion 


used, In 1953 1034F emulsion recorded this line well, but not the other [K tv} 


line at 6101 


In 1953 this line was regarded as too strong for [A tv] In 1954 it appears 
stronger than the other [A tv] lines ; probably He 11 7177 contributes 
(Plate 2(d)). Just resolved but measured as a single line 


Probably due to C tv as observed in planetary nebulae (see Swings and Jose, 


P.A.S.P., 61, 181, 1949) with [Ni 11] as a minor contributor. 


Representation of elements 


H ‘The Balmer series is represented from Ha to H 18(3691). 
A continuum appears on the Newtonian spectrogram 
extending well beyond the series limit. ‘The Paschen 
series is recorded from 8863 to 8500 A. ‘The stronger lines 
have very wide faint wings—see Hf on Plate 2 (c). 

New lines due to this element appear at 7816, 7281, 8361. 
4086 still appears as a very strong wide band and 5412 is 
also strongly winged. 

4267 is weakly present. 

4050 has been increasing in intensity since 1951. 

5801, 5812, a well resolved pair despite strong wings, 
became quite a prominent feature in 1954. 5801 appeared 
weakly in 1952 but was regarded as unidentified in Paper II. 
7378 appears weakly. 

Weakly present; multiplet 3. 

Sull prominent. 


: 
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HH 
He 
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Ou 
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[Cliv] 
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24! 
Multiplets 1, 2, 17. The 4640 group appeared as a very 
diffuse, unresolved band in 1952 (much stronger than in 
1951); im 1953-54 all Nut lines appeared much sharper, 
with 4634, 4640 and 4097, Hd as well resolved pairs. 

The strong multiplet 1 is the chief feature in the far 
ultra-violet. The very broad band (4048 to 4064) recorded 
as probably due to Niv in 1952 has disappeared. ‘This 
feature probably has the same origin as the broad Nu 
bands. 

‘The line 8446, apparently excited by Ly 8 (Bowen, P.A.S.P., 
59, 196), appears strongly in 1953 and 1954. ‘These 
infra-red plates are much better exposed than those of 
Paper I], and O1 8446 seems to have weakened relative 
to [Ol] 7324, in agreement with the general trend of 
increasing tonization. 

6300, 6363 are still prominent. In addition, the auroral 
line §577 appears to be weakly present. 

The pair at 7324 is prominent, but the 3727 pair cannot 
be detected even on the Newtonian spectrogram. No 
permitted lines are observed. 


Possibly represented by the one line 3757. 


5007, 4959 and 4363 are still very broad bands 
Not found. 
356g, 3967 
( Plate 2 (a)). 
(1954). 
Multiplet 4720, 
is apparently absent but might be expected to appear soon, 
‘The continued 45712, attributed 
doubtfully to Mgt in Paper I, means that some other 
identification should now be sought. 

4481 still appears weakly in 1954 

Represented by multiplets 1, 2, ?4, 5. 

Not found, rather surprisingly in view of the presence of 
Silt. 

Represented by 1F and ?2F. 


present as very strong 


bands, like [Ou], 
3343 appears on the Newtonian spectrogram 
with 


uncertain wave-lengths near 


presence of the lin 


‘The weakness or absence 
of 6717, 6731 closely parallels the weakness or absence of 
[O 11] 3727. 

6312 (3) persists; 3797, 3721 (2F) are both blended with 
H lines, but appear to be present judging by the Balmer 
decrement. 

5517, 5537 (1 F) were recorded in 1953 and 1954. 

A line measured four times in 1954 at §323°9 may be due to 
multiplet 3F (5322-2). 
Represented by 1F, 3F 
is one of the outstanding features of recent developments 
In Paper If 7135 (1) was weakly recorded. 

Represented by 2F. 
unidentified line in 1952. 


8046 (1) was measured once. 
The strengthening of these lines 


4740 was recorded weakly as an 
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Represented by 6101, 6793 (1F) and 4511 (2F). Only 
6793 was found in 1953. 
Not found. 
Doubtfully represented by 6992, 7153 (3/), but the 
wave-lengths are discrepant. 
Doubtfully represented by weak 5783 (1) and 5552 (2/) 
with discrepant wave-lengths. 
Doubtfully represented in 1954 by a very weak line once 
measured at 5071 (4/'). 
Possibly represented by 4529 (2/) which increased in 
intensity from 1953 to 1954. ‘This line is analogous to 
4058 of [Fem]. ‘The next strongest line is blended with 
Mg tt 4481. 
Very doubtfully represented by two members of 2/', but 
1 F should appear too. 
Still represented by very numerous lines. ‘The permitted 
lines have faded considerably. ‘The progressive change in 
ratio between permitted and forbidden lines is well 
illustrated by the lines 4233, 4244 whose intensities for 
each year since 1949 are as follows: 

1949 1950 1951 1952 1953 1954 
4233 4 4 8 6 3 3 
4244 a 2 6 7 8 10 
The absolute intensities are not of course strictly comparable 
from year to year; the intensities for 1949-50, dependent 
only on Newtonian spectra, should certainly be increased 
for comparison with later years. 
A line at 3954 is very doubtfully ascribed to this ton, but 
the more familiar line 4419 does not appear. 
Still represented by multiplets 1 3 F. 
A special search for the multiplets in green and infra-red, 
predicted by Edlén* without precise wave-lengths, was 
conducted in 1953 and 1954. A number of unidentified 
lines, found in 1953, were listed as possibly due to this ion. 
I am much indebted to Professor Edlén for communicating 
privately that the following lines in RR ‘Tel can be definitely 
ascribed to [Fetv] according to the independent analysis 
of his laboratory data: 


Observed wave-length Int. (1954) 
4368-2 
4900°2 8 
4906°2 12 
4 
7189 3n 
722! 3 


* B. kdlén, Paris Colloquium on Novae and White Dwarfs, 1939 
t A. DD. ‘Thackeray, The Observatory, 74, 9°, 1954 
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The lines 4868 and 4906 suffer from blending with HB 
and Fett respectively (Plate 2(c)) The other strong 
line 4900 has been observed in » Car (Paper I1) and in 
other objects. RR ‘Tel, with its sharp lines and continued 
brightness, is by far the best object in the sky for the 
astrophysical study of [Fe tv] lines. 

[Fev] Members of multiplets 1F, 3/ appeared weakly in 1954. 
All lines observed are listed in Table II, (Plate 2(a)). 

[Fevi}, [Fevit] Notfound. Nu 5679 masks [Fe vi] 5678. 

[Cott] The line at 4152 (Plate 2(4)), previously observed in 

7 Car and RR Tel, and tentatively ascribed in Paper II to 
[Cot] on the ground of good agreement in wave-length 
and of relatively high transition probability in the region 
studied, must now be regarded as unidentified in both 
objects, ‘The increase in strength in 1953 and 1954 
points to an ion of higher tonization potential, 
The few lines previously recorded have faded. In 1954, 
4201 alone was once recorded weakly. 
No coincidence with lines in RR‘Tel could be found with 
Velasco's unpublished laboratory lines due to this ion. 
The identification of some lines with [Nim] or [Nity] 
however remains an interesting possibility, 

The strongest unidentified lines are 4152, 4206, 4209, 5033, 5041, §233, 5288. 

4. Summary of developments in 1953-54. ~The chief characteristic of the 
development of RR ‘Tel since spectroscopic observations began in 1949 has been 
a gradual increase in the state of ionization. In 1953-54 this is exemplified well 
by the emergence of Civ, [Fetv], [Fev], [Au], [Atv], [Cliv], [Kiv]; there 
has been no very marked change in individual lines like the emergence of strong 
broad bands of [O11] and [Newt] which occurred between the 1951 and 1952 
scasons. 

A second feature is the fading of Fett permitted relative to forbidden lines 
(see preceding section). In this connection it is worth emphasizing the continued 
presence of Sitt lines but non-appearance of Sitti and Sitv. ‘Che permitted lines 
of higher stages of ionization which have recently appeared are confined to the 
light elements C, N, etc. The absence of Sitt 4128, 4130 suggests dilution of 
radiation. 

Thirdly, there is a striking difference in character between lines of low 
and high stage of 1onization or excitation. Fett lines are sharp, while lines such 
as [Fev], [K iv] 4511 are quite diffuse. (See Plate 2 (a, 6) for comparison of 
Fettand|{Fev]). Among the stronger lines, the series of increasing diffuseness 
is continued by Civ s8o1, [O1n) and [Nei] and the broadest band of all, He u 
4686. ‘Thus diffuseness appears to be correlated with both increasing ionization 
and excitation potential 

5. Rate of development compared with other slow novae It is well known that 
different novae tend to show roughly the same sequence of spectroscopic 
developments, but that the rate of development varies markedly. McLaughlin’*, 
in particular, has derived quantitative relationships between the rates of develop- 
ment and of decline in light and the absorption velocities of nova 

* D. B. MeLaughiin, Ap. J., 8g, 9602, 1937; Qt, 360, 1940; OS, 42%, 1942 
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Sufhcient spectroscopic observations of RR'T'el during the nebular stage 


have now been accumulated to permit a fair comparison with the corresponding 
stages of the well-observed slow nova RR Pic.* 


The observations of spectroscopic 
changes in R'I Ser and Car are very scanty, but nevertheless an attempt to 
compare these other two slow novae with RR ‘Tel is of some interest. 

‘The comparison of stages of spectroscopic development can be supplemented 
by two points on the light-curves. In the cases of RR Pic, RR Tel and RT Ser 


it is known that the epoch of maximum light was preceded by a steep rise some 
months or years earlier ; 


between these dates the light-curves are fairly flat at 
a magnitude or so fainter than the brightest maximum. 


The preliminary rise 
of RR Pic was not actually observed, but must have occurred between 1925 


January 13 and April 13. 7) Car was at its brightest known maximum in 


1843 and is known to have fluctuated between 2nd and 4th magnitude from 
1677 to 1830. 


The main object of the accompanying ‘Table III is to compare the various 


developments during the nebular stages. Consequently the zero-point of the 


tume scale is chosen to coincide with the “ Eta Carinae” stage ({Fe 11] predominant), 


although in » Car itself this stage seems to have persisted without much change 
for at least 40 years. 


Many adopted stages correspond to marginal appearances 
of various states of ionization, especially of Fe. 


‘The dates of marginal appearance 
are of course subject to very great uncertainty, especially when comparing spectra 
taken with widely different spectrographs. 


It must be remembered too that the 
marginal appearance of any emission line will depend upon the strength of any 
continuous spectrum; however, the continuum is usually very weak in these 
objects at the corresponding phases. ‘The marginal appearance of [Fe 1} 
in 7 Car is given as 1951, since this ion is well represented in ‘Table II of Paper II, 
while Gaviola’s observations! from 1944 to 1951°2 do not include lines such as 
4734, 4754 and there is no certain evidence of the ion at all in earlier years. 

As was pointed out in Section 3, the development of RR ‘Tel has involved 
the more rapid fading of Fett than of [Fe] lines. ‘The ratio of permitted to 


forbidden line intensities (e.g. 4233/4244) seems to vary considerably from nova 
to nova at the same stage of ionization. 


No doubt this ratio is an indication of 
a more complex set of physical conditions than the criteria of development used 
in ‘Table IL. 

The last column of ‘Table II] gives an average time-sequence, based on all 
four objects, for the development of the various stages, in arbitrary units. 

In big. 1 the data of ‘lable IIL are plotted with RR Pic used as a basic com- 
parison. It will be seen that there is a roughly linear relationship between the 
rates of development of all four novae, surprisingly close considering the 
uncertainties of the dates. ‘The slopes of the straight lines drawn through the 


points correspond to the following time-scales for corresponding stages of develop- 
ment of the four novae; 


RR Pr RR ‘Tel 


6 


RT Ser Car 
17 220 


‘The above figures are necessarily rough, their probable errors being of the order 
of per cent. 

Spencer Jone lun. ape Obs., | pt 

t H. Shapley, 823, 192 


25 


t k. Gaviola, Ap. 7., 118, 234, 1953 
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It 1s of interest to note that since McLaughlin found that N. Lac 1936 developed 
at 16 times the rate of KR Pic, this rate corresponds to about 3500 times that of 
) Car. However, no close parallelism can be expected between the developments 
of the fastest and the slowest novae; it is the slow novae like RR Tel that seem to 
exhibit the purest examples of sequences of increasing ionization. 


~ 5) 


—— 


1 1 L 
-200 100 0 +100 +200 +300 
RA Pic At Gays) 


big. 4 Epochs of corresponding stages of development of RR Pic (abscissae) and » Car (top), RT Ser 
(middle), RR Tel (bottom) 


The close relationships exhibited in Fig. 1 suggest that extrapolations of the 
lines may be attempted without falling into the realm of pure speculation. ‘Thus, 
by extrapolating the line for RR'Tel we find that the marginal appearance of 
{Fevi] and [Fevit] may be expected at about 1954°8 and 1955°8 respectively. 
No sign of [Fe vt] 5177 could be detected on the last plate of that region (1954°76) 
but its appearance in the near future followed by [Fe vit] may be confidently 
expected. 


3 +100" 

° 

“ 

e y 
~100 10 
A 

0 

a, 

a 

a x 


No. 3, 1955 and the development of slow novac 


247 

Extrapolation of the line for » Car suggests that the [O 111] lines should appear 
about 1980. As was pointed out in Paper Il, however, the lines of oxygen in 
7 Car seem to be surprisingly weak, It will therefore be of special interest to 
see whether the [O11] lines do follow the pattern of development suggested by 
other novae. 

The first two entries in Table II] refer to points on the light-curves and their 
relation to the other entries is probably insecure. Nevertheless it is of interest 
to attempt to estimate from them when a preliminary rise of » Carinae may 
have occurred. A straightforward application of the relative time-scales suggests 
a date between 1670 and 1710, with an uncertainty of about half acentury. Halley's 
observation in 1677, when the star was ranked as 4th magnitude, is usually accepted 
as the first record, and one might therefore suppose that a preliminary rise from 
pre-nova state occurred shortly before. However, Bayer’s U/ranometria (160%) 
apparently recorded the star as 2nd magnitude on second-hand evidence.* As 
pointed out by Innest, the star was probably fainter than 4th magnitude in 
Ptoiemy’s time, for it was not included in the A/magest. Any records of the star 
prior to 1677 would be valuable in this connection and the problem seems to 
deserve the attention of historians of astronomy. 


* Quoted by Agnes Clerke, Problems in Istroph ViICS, 


+ dnn Cape Os) IX (not XI as sometimes quoted) 


Rad: liffe Ohservatory, 
Pretoria : 


1955 /ebruary 2 


ADDENDUM 
Note added in proof I am deeply indebted to Dr A 


Armitage who has supplied the 
following information concerning the history of Eta Carimac 


based on Muller and Hartwig, 
Geschichte u. Lit. des Lichtwechsels, Leypzig, 1918, 1920 

The non-inclusion of Eta Carinae in the oldest Catalogues (Ptolemy, etc.) allows no 
conclusions to be drawn in. view of the fact that bright stars such as 7 Argus (10° farther 
north) are notincluded. ‘The earliest records of the region of Eta, depending on observations 
in India at the end of the sixteenth century, are included in Bayer's and v. Houtmanns’ 
charts. v. Houtmanns’ (1603) does not show the star, while he includes p, q. 5, pe Car, 
whence Eta can hardly have been brighter than 5". Bayer (1603) includes a group of some 
seven stars, designated d, of 4th magnitude; the modern Eta may have belonged to this 
group, but Bayer’s Eta Argus, of 2nd magnitude, is certainly not Eta ¢ 

Lacaille (1757) 


ar 
assigns 2nd magnitude to his Eta Argus, with an accurate position which 
clearly identifies it with the modern Eta Car 

One may conclude that the star must have brightened at some time between 1600 and 
1757, but that the richness of the Milky Way region renders naked-eye records practically 
unusable for estymating the extent of the brightening 


371, 1903 
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EMISSION OBJECTS IN THE SMALL MAGELLANIC CLOUD 
SHOWING THE N,, N, NEBULAR LINES 


E. M. Lindsay 


(Received 1955 March 41) 


Summary 


Out of a total of over 100 objects in the Small Magellanic Cloud discovered 
on ADH spectrum plates showing Ha emission, thirty-three have the N,, N, 
nebular lines in their spectra. These latter have 
yroups 


been divided into two 
The members of Group I, containing thirteen objects, are associated 
with nebulosity and are doubtless early-type stars exciting the surrounding 
nebulosity. Their mean absolute photographic magnitude is — 49+ 0°6 
(p.c.) and mean colour index (blue minus red photographic) — 0°47 + 0°19 
‘The members of Group II, containing seventeen objects, are not associated 


with nebulosity and have a mean absolute photographic magnitude of 


and mean colour index 0°77+ 0°29. ‘Three objects were too 
bright to determine magnitudes on the plates available. It is suggested that 
the members of Group II are planetary nebulae. An additional eleven 
objects showing Ha emission and with the same characteristics of Group II 
are possible planetaries 


A list of emission objects published by Miss Cannon (1) in 1933 contains 
only thirteen such objects as being then known to be present in the Small 
Magellanic Cloud, Of this number, eleven were classified as nebulous and two 
as O-type stars, 

Using a 10°5-in. visually corrected camera and 15° objective prism of the 
Mt Wilson and Palomar Observatories, the Lamont-Hussey Observatory of 
the University of Michigan have completed a survey of the Small Magellanic 
Cloud in Ha light. By this method 120 new emission objects have been discovered 
by Henize and Miller (2). Eighty-eight of these are classified as “stellar” 
and thirty-two as ‘‘nebulous”. As pointed out by the authors, criteria were 
lacking for further classification of the emission objects classed as stellar, so that 
although the nebulous objects would include only gaseous nebulae the stellar 
objects may include small nebulous knots, planetary nebulae and a variety of 
stellar types. 

A catalogue of nebulosities in the Small Cloud has been compiled by Virginia 
Nail, Whitney and Wade (3). ‘The instrument used was the 3 in. f.1°5 Zeiss 
Sonnar camera of the Boyden Station stopped to f 2, with a red glass filter 
(Corning 2403) and Ha interference filter. ‘There were also available for 
examination a go™ red (103a~E) plate, a 60™ blue (103a—O) and a red (103a~-E) 
objective prism plate, all taken with the ADH telescope. ‘The nebulosities were 
classified according to the system of Bok and Wade (4). In addition to their 
types, a Type V was added which indicates objects with stellar images. ‘This 
type contains 41 per cent of the total of 152 catalogued nebulosities. 
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In tests carried out by the writer on the ADH prism (angle 3° 10’; dispersion 
at Hy 240 A mm) attention was concentrated on two regions, the Large and the 
Small Magellanic Clouds. ‘The large amount of material made available by the 
tests enabled a survey to be made of bright-line objects in the Clouds. The 
survey is completed for the Small Cloud; additional plates are needed to complete 
the survey for the Large Cloud, Since it was known that a search for nebulosities 
in these regions was being made at Harvard, no attempt has been made to catalogue 
diffuse nebulosities or stellar but nebulous looking objects. So as to avoid 
overlapping of the work at Harvard and Armagh the survey has been confined, 
with a few exceptions, to stellar-looking objects showing emission lines in their 
spectra. 

A total of 106 emission objects has been catalogued for the Small Cloud, 
The present paper deals only with those which show the N,, N, nebular emission 
lines (AA 5007, 4959). The observational material is as follows 

Spectra. —The most fruitful discovery plates were red plates (103a~E) used 
in conjunction with a Wratten No. 25 A red filter. ‘The filter cuts off sharply 
all light to the red of Hz. Exposures varied from 30" to 2", but nearly all 
the emission-line objects are registered on a 60™ exposure plate of excellent 
quality. In addition to these, red plates (103a—E) without filter and blue plates 
(103a~O and Ila-O) were available with exposures up to go", including widened 
and unwidened spectra. Guiding was carried out by means of a “ gooseneck” 
micrometer eyepiece fitted to the 8-in. guiding telescope attached to the main 
tube of the ADH. Widening was obtained by guiding on the cross-wire of the 
eyepiece as for direct photography, periodically offsetting the cross-wire by a 
known number of turns of the micrometer screw and immediately bringing back 
the guiding star to the cross-wire by the slow-motion guiding mechanism, 

Vagnitudes.—Blue photographic magnitudes have been determined from 
eight ADH Ila—O plates with exposures ranging from 5° to 1"; red magnitudes 
from six ADH 103a-E plates with the 25 A red filter, the exposures ranging from 
1 to go™. ‘These plates were taken on a programme for Dr Shapley and I am 
grateful to him for permission to use them first at Armagh. 1 he magnitudes 
were estimated by means of a scale of stellar images taken in the blue and red 
with the ADH. ‘The scale was calibrated for each plate by means of King’s 
unpublished magnitudes of some of the stars in Sequence | of the Small Cloud (§). 
King’s magnitudes were determined with the photoelectric photometer attached 
to the 60-in. Rockefeller reflector, and are reduced to the International Seale. 


His faintest magnitudes in the blue and red are m,= 15:06 and m,, = 15°64 
respectively, which are approximately the limiting magnitudes on a 30" and 5™ 
exposure in the blue and red respectively with the ADH. ‘The mean difference 


between King’s magnitudes and those given in H.4., 108 for the four faintest 
stars (t, v, wand x) in King’s list is (7.4.—King) = +0™-20. Stars fainter than 
those in H.A., 108 were corrected by —o™-2 and used for the long-exposure 
blue plates. King’s red magnitudes are the only ones available in the Small 
Cloud. It was therefore necessary sometimes to extrapolate the calibration 
curves slightly to reach beyond his limit. In some cases small corrections, not 
exceeding 0-2, were applied to the measured magnitudes for distance from the 
standard stars in accordance with the preliminary estimates for distance correction 
found by the writer (6), 
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Che objects showing nebular lines are listed in Table I, which gives the 
catalogue numbers, right ascensions and declinations (1g00), the positions in 
X and Y coordinates as described in Harvard Annals, 60, No. 4 (1907), and 


spectral lines. 


The right ascensions and declinations are rough measurements 


derived from graph paper with the use of a few standard stars. 


The objects in Table I have been re-arranged in ‘Table II which gives the 


blue and red photographic magnitudes, the colour indices (blue minus red 


photographic), the photographic 


modulus (m 


No R.A Dex 


m 


Xx 


absolute 


Taare | 
Objects in the Small Magellame Cloud 


y 


Emission lines in addition to Ha, N,, Ny 


magnitudes and remarks. 
VW) of the Cloud has been taken to be 18-9 photographic (7). 
correction has been applied to allow for space absorption. 


with nebular lines in thetr spectra 


‘ » 16°79 6552 13260 Hf, Hy 
7% 21°9 9204 g642 Hf 

37" 74 22'3 9522 59156 Hf 

4 o 73 9792 10092 

5 39°9 73 8652 Hf, Hy, H4, He, HZ, A 3727 

6 41°3 3 49°79 10334 7800 Hy 

7 o 41°% 3 56% 10440 8520 Hf, Hy 

8 o 41°% 10464 7428 Hf, Hy, Hd, A 3727 

9 42°9 73 39°'0 10704 8454 
10 44° 7% 11142 7944 Hf, Hy 
45°9 40°5 11490 7242 Hf, A 4363, Hy 
12 460° 74 11640 6180 Hf, Hy 
13 » 47°6 74 31°93 11976 5388 Hf, Hy, Ho 
14 o 73 40°2 12072 Hy 
i o 48 7% 24'0 12072 9486 Hf, Hy, H6, A 3727 
16 71 §7°1 12012 14628 Hy 
‘ 73 13596 10548 Hf, Hy, Hd, A 3727 
18 54°9 73 13806 10158 Hf, Hy, Hd, A 3727 
19 73° 12°9 13848 10098 Hf, Hy, HO, A 3727 
72 00'S 13926 10848 Hf 

} 14106 11535 


1407! 


14610 13980 Hf, Hy, Ho, He, A 3727 

15108 12894 Hf, Hy, Hd, A 3727 

1 06 73 44 16650 094 Hf, A 4363, Hy, Hd, He, HG, A 3727 
26 73 50°0 32 620 Hf, Hy, HO, He, HC, A 3727 
1 73 47°9 17892 7674 Hf, Hy, Ho, A 3727 

73 47°5 17955 7752 Hf, Hy, Ho 
73 18036 7578 Hf, A 4686, Hy, Ho, A 4727 
1 12°0 73 «5! 18054 7578 Hf, A 3727 
7 125 516 18072 7494 HA, Hy, HO, A 3727 
42 1 184 73 464 19740 7650 Hf, Hy, Ho 
33 1 73 40 20535 7896) A 4686, A 4363, Hy, Ho, He, H¢, 

A 3727 


but also in their absolut 


negative; in Group IL positive. 


the numbers in each group 


magnitudes, 


‘Table II includes two types of objects differing chiefly in their colour indices 
In Group I the colour indices are 
‘Table IIL gives the mean colour indices and 
mean absolute magnitudes of each group together with their probable errors and 
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No. 3, 1955 Magellanic Cloud showing the \,, \, nebular lines 


Il 


Magnitudes, red indices and absolute photographic magnitudes of the oljects of Table I 


Remarks 


7 14°02 14°57 ‘ 4°9 He large and diffuse. Region of extended 
nebulosity 
8 14°33 14°61 Hae strong and diffuse. Star in a concentration 
of an oval-shaped nebulosity 
13°2 14°10 ‘ H« large and diffuse. Star involved in nebulosity 
which is connected with the nebulosity of No. 8 
14°67 c 7 He» diffuse and oval. Small cluster probably 
with nebulosity 
15°41 ‘ Circular nebulosity with central star An His 
sphere. Harvard No. 97 
15°41 Very nebulous region 
14°13 4 Very nebulous region 
15°33 ' Nearly symmetrical circular nebulosity with 
central star. An Hi sphere. Hea is double 
On long-exposure red plates there is adjacent 
faint nebulosity. This is No. 188 in the 
Shapley-Wilson Catalogue (8) where it is 
described as ‘“‘ nebulous star, planetary 
Very nebulous region 
Circular nebulosity with central star An Hoes 
sphere. Harvard No. 149 
Very nebulous region 
Very nebulous region 
Very nebulous region 


Ha very faint. No continuous spectrum 
Ha very faint. No continuous spectrum 
Ha very faint No continuous spectrum 
Ha very faint. No continuous spectrum 


Ha moderately strong. | No continuous spectrum 
Ha broad and diffuse, possibly double. No 
continuous spectrum Ov erlapping a fainter 
star. Very nebulous on red plates. Harvard 


No. 16, described as a“ bright nebulous star 
of Type V 

Ha strong No continuous spectrum 
though situated between an Hit sphere and 
extended nebulosity it is apparently not 
associated with any nebulosity 

Ha moderately strong No continuous spectrum 

Ha moderately weak. No continuous spectrum 

Ha moderately strong No continuous spectrum 

Ha moderately strong Star nebulous Faint 
continuous spectrum 

Ha moderately faint No continuous spectrum 

Ha moderately strong. No continuous spectrum 

Nebulous object in neighbourhood of five stars 
There is an impression of faint nebulosity but 
60! exposures on blue and red plates taken 
with the 6o0-in. reflector show that at least part 
of this is due to fainter stars Harvard No. gy 
(X — 13839; Y — 10081) is described as“ nebu 
lous group, three nebulous stars He fairl: 
strong and no continuous spectrum 


Ha moderately strong No continuous spectrum 
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251 
15°92 15°67 0°25 3:0 
2 16°19 0°45 2°7 
4 16°21 15°79 2°7 
14°40 13°93 0°47 4°5 
10 13°90 } 1°36 
14°60 13°07 1°53 4°3 
13 15°66 14°74 3°2 
is 14°01 14°26 4°3 
6 15°16 0°43 
17 14°23 13°27 4°7 
19 14°04 14°01 0°63 4°3 
is6s 14°37 +1°28 


ak 


E. M. Lindsay, Emission objects in the Small 
Tasie II (cont.) 

No. Moy Mor CJ. Remarks 

23 #13°77) Nebulous object with possibly traces of nebulosity 
on long-exposure red plates. Not separated 
from a fainter star. Ha moderately strong 
with no continuous spectrum. Harvard No 
107, described as “‘ nebulous star ’’ of Type V. 
No. 182 in the Shapley-Wilson Catalogue and 
described as “ nebulous star, planetary ?”’ 

32 15°38 0°42 3:1 Hea moderately strong. Possibly faint continu- 

ous spectrum 


Vol. 115 


Notes on Table 11. No. 25 (IC 1644) is too bright to measure on the plates available 
(exposures 30™ blue and g0™ red). It is classified in H.A., 60 as a gaseous nebula. ‘There 
is no extended nebulosity in the immediate neighbourhood and it is difficult to say whether 
or not some adjacent faint stars give the impression of faint nebulosity to the object. 

The emission from No. 26 comes from a condensation of stars and nebulosity. No 
measurements of magnitudes could be made 

No. 33 was off the field of any of the blue plates available. The value m),=11°95 was 
obtained from a single plate, and can only be approximate. It is described in H.C. 224 as 
a gaseous nebula. On a 90™ exposure ADH red plate with red filter it is a stellar object 
in a small cluster of much fainter stars. ‘There is no noticeable nebulosity around it or 
in the neighbourhood. One might suspect a trace of faint nebulosity but, if so, it could as 
readily be interpreted as faint unresolved stars 


Taste Ill 


Mean values of Table LI entries 


(sroup M No Remarks 
I ~o'47 0°19 —4'9+ 06 13 Associated with nebulosity 
Il +0°77 + 0°29 17 Not associated with nebulosity 


It seems clear that in Group I we are dealing with O- or early B-type stars 
which are exciting the surrounding nebulosity and giving the typical emission 
spectra of gaseous nebulae. It is impossible to classify the types since all are 
too faint or are in dense regions. Of the thirteen objects, two (Nos. 21 and 22) 
are in NGC 446 and five (Nos. 27~31) are in the neighbourhood of NGC 456 
and NGC 465, confused regions with overlapping spectra and nebulosity. 
Three are circular nebulosities with a central star—so-called H 1 spheres. Data 
on these latter are given in Table [V and compared with the measurements of 


Taste IV 
11 spheres in the Small Magellanic Cloud 
(A~wArmagh; H-= Harvard) 


Cat. No. Moy M 
A H A H A H A H 
75 15°00 14°90 39 
‘ 79 15°60 —O'lg eve —3°3 
18 97 14°64 14°93 —4°3 
24 115 14°83 15°25 
119 1§°7 15°42 o's + 0°94 —32 —3°S 
28 (149 15°16 14°95 0°34 


Means (omit 79) 14°88 0°54 +0'06 —4°0 —4'0 
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Virginia Nail, Whitney and Wade referred to above. ‘The latter list seven such 
objects. ‘The remaining four are not contained in Table I since they show only 
Ha emission. No measurements have been made on Harvard No. 79 since it 
appears on the Armagh plates as a nebulous object but with no trace of a central 
star. Harvard No. 75 was classified as probably a stellar object in a small cluster 
of faint unresolved stars but possibly nebulosity. 

The Armagh and Harvard blue photographic magnitudes are in fair agreement. 
In four cases, however, the colour indices differ radically. Possibly part of this 
represents the difficulty of estimating magnitudes of stars involved in nebulosity. 
The divergences between the estimates by different observers of the magnitude 
of a comet with a strong coma are well known, and a good example of the difficulty 
is seen in the Armagh estimates of (for instance) the central star of the first object 
in Table IV. From exposure times of 20%, 60%, 3, 10™, 30™ and 60™ the estimated 
magnitudes were 13°55, 14°05, 13°87, 14°18, 14°95 and 15°20 respectively. In the 
first four cases little or no nebulosity was evident and the mean of these was taken. 
On the last two exposures the nebulosity was strong and very strong and it was 
extremely difficult with this subjective method to estimate the magnitude of the 
embedded star. It was found in all cases that when the nebulosity became 
sufficiently dense to be troublesome, the estimates of the magnitudes became 
systematically fainter with increasing density. 

Referring to Table I, it will be seen that three of the objects of ‘Table IV 
(18, 24 and 28) show the N,, N, lines and two also show A3727. It seems unlikely 
that these lines would appear with a central star of positive colour index. It is 
unexpected that the brightest star absolutely (Harvard No. 30; M-= —5:0) 
gives only an Ha sphere with no other lines present. 

In Group II we have objects not associated with nebulosity, generally stellar 
in appearance, with little or no continuous spectrum, showing the typical discrete 
nebular lines and having positive colour indices. It is suggested that these are 
planetary nebulae. ‘The values of the absolute magnitudes are higher than those 
found for planetaries in our own galaxy. Using data on mean parallaxes based 
on proper motions and the rotational terms of galactic rotation, Berman (9) 
and Parenago (10) have arrived at values of —1-5 and +02 respectively for 
the mean absolute magnitudes. ‘The higher value of M,,= —3-7 found here 
may be due in part to the uncertainty in applying corrections for absorption of 
light in our galaxy and in part to a selection effect in the Small Cloud, where only 
the brightest planetaries would be found in the present material. With M= —1°5 
at the distance of the Cloud we would have m= 17-4, which is fainter than any of 
the emission objects so far found. 

The large positive colour indices are in keeping with the results of 
Vorontsov-Velyaminov (11) who has found a visual minus photographic index of 
+1°4 for planetary nebulae from integrated magnitudes. ‘This is explained by 
the composition of spectra of planetaries being mainly discrete emission lines and 
by the strength of the green nebular lines N,, Ny. ‘The strength of Ha will 
also explain the large positive colour indices derived from the difference between 
photographic and red magnitudes. 

A number of objects have been found which have the same characteristics 
as the planetary nebulae above, namely Ha emission, little or no continuous 
spectrum, positive colour indices and comparable absolute magnitudes. ‘These 
are listed in Table V with additional details in Table VI. 
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None of the objects in Tables V and VI is associated with any observable 
nebulosity. ‘Two (Nos. 34 and 36) or possibly three (No. 44) are in or in line 
with small clusters. ‘Taking all together, however, the mean colour indices and 
mean absolute magnitudes are +o80+012 (p.c.) and —33+0°5 (p.e.) 
respectively. ‘The mean colour index is naturally close to the value above, 


Taste \ 


Otnects in the Small Magellame Cloud which are possibly planetary nebulae 
No R.A Xx 


m 


Emission lines 


4 o 41°6 73 43°2 102484 8190 Ha faint, diffuse. No continuous spectrum 

15 44°9 73 11244 %946 He faint. Suspected N,, N,. at 
plate limit. No continuum 

36 © 45°4 74 04°0 11310 7446 He faint, diffuse. No continuum 

17 o 46°2 73 42°8 11568 7584 Ha moderately weak and diffuse. Faint 
continuum 

o 47°% 72 589 11946 10908 Ha weak. Suspected N,, N, at 
plate limit. Possibly faint continuum 

44°97 73 $4'! 12138 7614 Hea faint. Very faint continuum 

40 © 50°6 73 14°8 12096 997% Hae strong. Hf, A 4686, Hy. Generally 


A 4686 is outstanding, but on one plate 
its intensity is equalled by that of H/ 
with Hy also strong. No continuous 
spectrum 


41 © 50°7 73 02°2 12726 10740 Ha weak and diffuse No cont:quum 

42 72 366 12990 Hafairly strong. Suspected Hf, Hy, Ho 
at plate limit. No continuum ? 

43 © $5°9 72 38°4 14094 11772 He faint. No continuum 

44 1 12°9 73. 4t°s 18354 So82 Ha moderately strong and diffuse. No 


continuum ? 


Taste VI 
absolute photographic magnitudes of the objects of Table | 


No Mog Moy C.i Mog Remarks 


Vlagmtudes 


red indices and 


34 16°37 14°59 o-78 3°5 In a small cluster. No nebulosity on go™ red or 
60 blue plate. Harvard No. 15 described as 
‘ nebulosity around four stars '’ of Type I 

1s 16°29) 26 Somewhat nebulous faint star 

% 15°46 15°09 0°37 +4 Inasmall cluster. Some doubt of this object 

14°90 14°34 179 Object nebulous 


of 


Near to but not to be confused with Harvard 
No. 79 which shows no emission lines 

14°73 0°22 

42 1§°72 14°09 1°64 

+3 16°04 


vw 


14 16°23 «15°48 2°7 Five stars close together. "%ifficult to be sure 


which one gives emission 


since the objects were chosen partly for their positive indices; the mean absolute 
magnitude 1s somewhat lower, which is to be expected if they are too faint to show 
emission lines other than Ha. ‘Taking all such objects in Tables I] and V together 
we get C.l.= +078+025 and M,, 35+05. Their distribution with 


respect to absolute magnitude is given in Table VII. The sharp drop in the 
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Distribution of probable (o) and possible (+) planetary nebulae in the Small Magellame Cloud 
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No. 3, 1955 Magellanic Cloud showing the N,, N, nebular lines 


Taste VII 
Distribution of planetary nebulae with respect to absolute magnitude 


Limutsof 50-46 4°5-41 4°0-3°6 34°5 
No 2 5 


5 


numbers at the faintest absolute magnitudes is doubtless due to working at the 
lumit of the plates. Longer-exposure guided spectrum plates with the ADH will 
most probably reveal fainter objects. No. 33 has not been included in Table VII. 
If it were a planetary nebula with the average colour index of + 08 its absolute 
photographic magnitude would be — 60. 

The positions of the objects suspected to be planetaries are plotted on 
Plate 3. No. 33 1s just off the eastern edge of the plate. If it 1s included, the 
longest diameter is about 5° ; the diameter at right angles to this ishalf that. There 
was not sufficient material for a complete search over the whole area of the Cloud 
as given by Shapley (12) from microphotometer tracings. One spectrum plate 
im the direction of the Large Cloud revealed no objects which could be classified 
as planetaries. 

drmagh Observatory, 

Armagh, N. Ireland : 
March 29 
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FRAUNHOFER-LINE PROFILES 
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Summary 


From the experimentally determined cross-sections of the Na atom for 
collisions of the second kind it is shown that the hypothesis of true absorption 
is probably not applicable to the formation of Fraunhofer lines in the solar 
photosphere. Using atomic transition probabilities for radiative excitation 
and ionization, values of the source functions and the ratio of the line to the 
continuous absorption coefficients are derived in terms of kinetic, excitation 
and ionization temperatures. ‘The functional dependence of these quantities 
on the optical depth is such as to permit a numerical integration of the exact 
equation of transfer, and the method of effecting this integration is outlined 


Introduction.—-A_ satisfactory interpretation of Fraunhoter-line profiles, 
involving the minimum of hypothesis, requires a numerical integration of the 
equation of transfer for line formation. Unfortunately the form of the equation 
itself has not yet been settled. ‘Thus in recent years, following Unsdéld’s 
adoption (19) of the heuristic and simplifying assumption of true absorption 
for the formation of lines in the spectrum of 7 Scorpii and the justification he 
offered for its application to the solar photosphere (20), some workers have 
applied the resulting simple equation of transfer to find new values for the 
temperature distribution in the upper photosphere from the limb darkening of 
Fraunhofer lines. K.-H. Béhm (2) has, however, recently shown that, even 
allowing for possible thermal inhomogencities due to granulation and for 
departures from local thermodynamic equilibrium, the thus derived temperature 
distributions fail to account for the limb-darkening both of faint lines and of the 
wings of the Na D-lines. 

Invaluable though Unséld’s hypothesis has proved and is likely to prove in 
the interpretation of stellar-line profiles, it is probably not strictly applicable to 
the solar photosphere. It is the purpose of this paper to propose a method, a 
refinement of one earlier given (13), whereby a numerical integration of the exact 
equation of transfer can be effected. ‘he paper is divided into three sections, 
the first of which examines the hypothesis of true absorption. In the second 
section values of the absorption and emission coefficients in the exact equation of 
transfer are derived, and in the final section these values are shown to lead to a 
simple method of numerical integration of the equation. 

The essential part of the paper is contained in the second section. Upon its 
results, which must be demonstrably correct, depends the method of integration. 
It is suggested, however, that the simple ideas involved may be most readily 
grasped if the paper is read in the order Section 1, Section 2.1 and Section 3; 
reference can then be made back as necessary to the remaining sub-sections of 
Section 2. 
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No. 3, 1955 Interpretation of Fraunhofer-line profiles 


1 Hypothesis of true absorption 

1.1. Test with the Na atom.—\n the strict sense of the words true absorption 
means that the line quantum absorbed is converted by a collision of the second 
kind then and there into kinetic energy of the free electrons and atoms at that 
place in the photosphere. If o,(e) be the cross-sectional diameter for a collision 
of the second kind, then the number of such collisions experienced per see by an 
atom of mass m, is given by 


ann + =) RT) de, (4.1) 
where n, 1s the number per cm* of colliding particles, each of mass my, and « is 
the relative kinetic energy of the particles. From the experiments of Christoph 
(4), see Bates, Fundaminsky, Leech and Massey (1), the values of o,(¢ + «,,) the 
cross-sectional diameter of the Na atom for collisions of the first kind with 
electrons are known, where «,, 1s the excitation energy of the 3*P from the ground 
state 37S. By detailed balancing (6) we find that 


T 


Yo 


)= ae + €y9), (1.2) 


where q, 18 the statistical weight of the lower state, g, of the upper. 

For an optical depth r,(5485 A) = 0°48 in the continuous spectrum, we find 
for a particular model of the solar photosphere (14), the following values of the 
temperature, total gas pressure and partial pressure due to free electrons 

T = 5838 deg. K 

p= 69 1ofdyne cm * at To 5485 A) = 0-48. 

p.=17°7 dyne ecm # 
Not widely dissimilar values are given by other photospheric models for a 
corresponding optical depth. ‘Then from Christoph’s corrected experimental 
values of o,*, converted to a," by equation (1.2) and a numerical integration of 
the expression (1.1), the number of collisions of the second kind* experienced by 
a Na atom in the 3’P state at this depth in the photosphere is 1°52 x 10% sec '. 
The average life of the excited Na atom in so far as it is determined by collisions 
of the second kind is therefore t,<0°66% 10 “sec. But the average life as 
determined by spontaneous emission is 1, A,,, so that the probability that the 
atom does not experience a collision of the second kind in this lifetime is 
exp(—1/4yfo). From ‘Trumpy’s determination of the oscillator strength of the 
Na D-lines (17), A,, = 0°62 « 10" sec |, or the probability that the atom does not 
experience a collision of the second kind in the lifetime 1, A,, sec is 0-976. In 
other words only 2-4 per cent of the excited Na atoms are quenched by a collision 
of the second kind converting the quantum previously absorbed into kinetic energy 
of the free electrons. The percentage is in fact even smaller since the life before 
spontaneous emission is determined by the roughly ten times shorter life produced 
by collisional damping. On this basis it is readily found that only 0°16 per cent 
of the excited Na atoms give up their energy as a result of collisions of the second 
kind. 

Since at this depth in the photosphere the hydrogen atoms are some 3-9 x 104 
times more abundant than the free electrons, it might be expected that the 
excited Na atom would experience many collisions of the second kind with H 
* ‘The number of the first kind is 7-90 = 10* see™', 
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atoms. As Massey and Burhop (8) point out, however, this neglects the fact 
that atoms approach one another relatively slowly ; as a consequence such collisions 
are nearly adiabatic, the excited atom having a chance to adjust itself without a 
quantum transition taking place. In fact collisions with H atoms lead to a 
perturbation of the energy level and to collisional broadening, not to collisions of 
the second kind. ‘lhe probability of a quenching collision will be small (8) if 
0€JhV > 1, where o is the gas-kinetic cross-sectional diameter of the atom and 
the relative velocity from Rosin and Rabi's value of o* = 41-38 x 10° 
for the Na atom (15) and a mean velocity of H atoms of 1-2 x 10°cmsec ', this 
quantity equals 27°2 and the encounters are nearly adiabatic. 

\ more direct comparison may be made with experiment. Norrish and 
Smith (12) have measured o,* for Na atoms in the 3*P state in atmospheres of 
various monatomic rare gases and of many molecules. For an atmosphere of 
He they find, for example, that o,*~o, and for methane =o11 x 10 “cm*, 
In an atmosphere of H atoms, as in the photosphere, the value would presumably 


he nearer that for a He than a methane atmosphere. Assume, however, that it 
. is comparable with the latter and is independent of relative kinetic energy. 
‘Then from equation (1.1) the number of quenching collisions experienced per 
& sec by an excited Na atom will be 

( ) 3°35 10% sec! 

my Ms 

y due to encounters with H atoms at this depth in the photosphere. If, as in the 
; previous discussion, the interval before spontaneous emission is determined by 
Z collisional broadening, we find that only 0-46 per cent of the excited Na atoms 
: experience a collision of the second kind with H atoms 

; 1.2. Runaway quanta. \n so far as the behaviour of the Na atom 1s typical, 


we may conclude that collisions of the second kind play no significant part in 
converting the exciting light quantum then and there into kinetic energy of free 
electrons and H atoms. Nor, so far as | am aware, has it ever been suggested 


that they do. What has been suggested is that the originally absorbed quantum, 
or one of different frequency, is re-emitted and re-absorbed many times until a 
collision of the second kind or a photo-ionization takes place. 


We are interested, however, in finding the source function at a particular 


place. In most observational investigations this ‘ place’’ is somewhat narrowly 
defined, and consists of a region with a superficial area of perhaps five seconds of 
arc square, or y600km square, with an optical depth of o-1 or smaller. The 


quantum, pictured as being handed about from atom to atom until a true absorp- 
tion takes place, has a mean free path, as Unsdld (18) has pointed out in another 
connection, with an optical depth in the line absorption coefficient of the order of 
unity, Assuming that the probability of a collision of the second kind with free 
electrons and H atoms is of the order of ovo1, as indicated in the previous sub- 
section, then after 100 transfers from atom to atom there will from the Poisson 
distribution be a probability of 1 /¢ that a collision of the second kind has taken 
place. In its “drunkard’s walk” of 100 paces, each of optical depth unity, the 
33 quantum will on the average be 10 units of optical depth away from the place 
4 where we are determining the source function before there is a probability 1 e 

that a collision of the second kind occurs. In short, as far as this particular region 
of the photosphere is concerned, the quantum has been scattered, not absorbed. 
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2. The exact equation of transfer and the coefficients of absorption and emission 
2.1. The equation of transfer.—The Fraunhofer lines are usually regarded as 

formed in accordance with the Eddington—Milne equation of transfer 


dl du 
cos (ko | on (Ky 


where « is the fraction of excited atoms experiencing a collision of the second kind. 
Obviously the equation is much simplified when « but this implies truc 
absorption at this particular level which, as we have seen in Section 1, is improb 
able. On the other hand, if « #1, the equation implies that the scattering 1s 
strictly coherent-——a physically impossible assumption which, from Houtgast's 
investigation (7), is not realized in the photospher Moreover the only 
quanta emitted in the frequency of the line are assumed to arise from a population 
of the upper state solely determined by excitation by quanta of the line under 
consideration 

There are therefore good grounds for preferring the exact equation of transfer 
(13) for line formation 

dl, 
(ky + Jo (2.1) 

where «x, and «, are the mass coefficients of absorption in the continuum and the 
line, both for the frequency v, and where ;, and /, are the mass coefficients of 
emission for the frequency v in continuum and line. ‘The latter coefficients are 
defined as the isotropic flux emitted at the frequency » per unit frequency per 
unit solid angle per unit mass. Attention is therefore always focused on the single 
frequency v and no mechanism of absorption or emission ts pre-supposed, 
Introducing the optical depth at the frequency v in the continuum 


Ty } ds, (2.2) 
the exact equation takes the form 


where the source functions, jg «, in the continuum and j,«, in the line, and 
) =*,/Kq are each functions of the optical depth 7). Without some knowledge 
of these functions, however, equation (2.3) 1s numerically almost unintegrable 
Accordingly in the remainder of this section an attempt will be made to evaluate 
these functions in terms of atomic transition coefficients 
2.2. Atomic transition coefficients for continuous absorption and emission. 

Following Milne (g) the atomic transition coefficients for photo-1onization and 
recombination may be readily derived. ‘The number of photo-ionizations from 
the rth state of the neutral atom per unit volume per unit time ts 


nab Avyu, dv, (2.4) 


where #, 1s the number of such atoms per unit volume, u, dv 1s the radiative energy 
density and ,(v) is the atomic probability of photo-ionization. The velocity V 
of the eyected electron is given by 


hv = \mV*~E,, (2.5) 


260 H. H. Plaskett Vol. 115 


where E, is the (negative) energy of the rth state. Correspondingly the number 
of recombinations of V-electrons per unit volume per unit time is 


where n,' is the number of ionized atoms in the ground state per unit volume and 
n, is the number of free electrons with the correct spin per unit volume. yp(V)d¥ 
is the fraction of electrons moving with a velocity between V and V +dV per 
unit solid angle, or 


u(V)dV =(m exp(—mV*2kT)Vi dV. (2.7) 
o{V) and s{V) are the spontaneous and stimulated cross-sectional diameters of 
the ionized atom for capturing an electron moving with velocity V on the rth 
state. ‘Then in thermodynamical equilibrium at temperature 7° we have detailed 
balancing so that dn,/dt=dn,'/dt. Solving for u, and recalling that it is given 
by Planck’s law we find from equations (2.4) and (2.6) 
a (AV) n, fv) di 


From the Boltzmann equation 


“y= 


Sah 
|  lexp (hv/k7) (2.8) 


q, K.—E.veT" ny 
and from Saha’s ionization equation 
n 2(2rm)** b'(T)_, 
(kT) * exp(E,/kT), 


where p,= 2n,k7, since the electron has two possible orientations in an applied 
magnetic field. 4(7) and 6‘(7) are the partition functions for the neutral and 
1onized atom respectively. Recalling from equation (2.5) that dv/dV =mV/h we 
find after a little reduction the following relations between the transition coeffic- 
sents for photo-ionization and recombination 


Since these are strictly atomic constants, they and the relations between 
them remain unchanged no matter what the external conditions in which the atom 
finds itself. ‘The atomic absorption coefficient for continuous radiation is thus 
always given by 

(hy cyb(v). (2.10) 

Consequently if \, be the number of neutral atoms in state ry per unit mass of 
the photosphere, the continuous mass coefficient of absorption due to these 
atoms 1s 


Ko= — Nab Av). (2.11) 


Similarly uf \,’ ws the number of ionized atoms in the ground state per unit mass 
of the photosphere, the mass coefficient of emission for the continuum is from 
equation (2.6) 


hy dV 
— N, Vo,4V) (2.12) 


neglecting the stumulated transitions. 


& 
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2.3. The atomic line transition coefficients.—As Spitzer (16), following 
Weisskopf and Wigner, has pointed out, it is statistically correct when dealing 
with an assemblage of atoms to regard absorption and emission lines as originating 
from transitions between broadened energy levels. ‘The probability of an atom 
finding itself in a sub-level between s, and s, + ds,, measured in frequency units, 
from the centre of the level is given (for atoms in an enclosure in thermodynamical 
equilibrium) by the normalized expression 


ds 


r 

7 (2.13) 
4, 1s the half half-width of the broadened energy level and is given by 5, = 1/4zt,, 
where ¢, is the mean lifetime in state r as determined by radiative and collisional 
damping. ‘The atomic probability of transition between a particular sub-level, 
s,, of the ground state and a level between s, and s, + ds, of state 2 is 
ds, 


By — + (2.14) 


The spontaneous and stimulated downward transition probabilities from a 
particular sub-level, s,, to a sub-level between 5s, and s, + ds, are likewise given by 


ds 
(2.45) 


When the atoms are not in an enclosure in thermodynamic equilibrium, 
these transition probabilities, since they are purely atomic quantities, remain 
unchanged. On the other hand, though the potentially realizable sub-levels 
are still in existence, the actual distribution of atoms among these sub-levels is 
no longer given by expression (2.13). Under these circumstances let the actual 
probability in the ground state be given by the normalized expression W(s,) ds,. 
Then if \,* is the total number of metal atoms in state 1 per unit mass of photo- 
spheric matter, the mass coefficient of absorption for atoms at rest is given by 

In these expressions v,,° is the constant frequency for a transition between the 
centres of the two broadened states, and the integrations are to be carried out 


(2.16) 


2 


subject to the condition that v,, = v,," + Av is a constant. 

he atoms are actually in motion, and for an atom with a velocity component, u, 
in the line of sight the central frequency of the absorbed line is 44° ~ u(v4,4°/c). 
‘The fraction of atoms with velocity components between u and u + du is 

(m®* exp(—m*u*® 2kT,) du, 

where 7, 1s the kinetic temperature and m* the mass of the metal atom. From 
the expression in (2.16) the mass coefficient of absorption for the frequency 
*+ Av will clearly be given by 


hv,," 


Vig 


+ Av) = exp (— m*u* 2hT',) du 


m® 
) ? 


W(s,)ds, 


Substituting 
gop 


(2.17) 
= ; (2.1%) 
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where 6 1s a species of half width for Doppler broadening, and 


a=(s, + Avy db, (2.19) 
equation (2.17) becomes, inverting order of integration, 
+ Av) 12 Bis | W(s,) ds, | (2.20) 


where o 1 a constant for the first integration. 
The first integral in equation (2.20) is well known (10), and provided 6,/b 1s 
small (under photospheric conditions it is of the order of 1-4 x 10 *) we may write 
exp(— y*)dy 


a F(a); 2. 
| (6, bP —20F(a)}, (2.21) 


where exp(—o*)| exp(x*)dx is a tabulated function. Hence from 


equations (2.20), (2.21) and (2.19), the mass coefficient of absorption takes the 
form: 


hv,," [2b s, +Av\?2 


s,+Av fs, + Av 
The integral in equation (2.22) can be written 


+ Av\? + Av\ | 


where, in order to estimate the relative magnitude of the two integrals, we have 
assumed W(s,) to be a delta function. ‘Then from the tabulated values of the 
F(Av b) tunction, a value of b= 477 « 10" sec ' (=0°044A corresponding to a 
Na atom at 6000deg.K) and 4,= 5:08 x 10’ sec ' (=o-00059A) we find for 
\v « 14h that the first term is at least 33°3 times greater than the second, and 
for Av » 3:4) that the second term is at least 81-2 times greater than the first. 
Further for large values of Av b 


#0 that for \v » 64 the error in putting this expression equal to — (Av +)? 1s less 
than 4 per cent. Hence with an error not in excess of this amount we may write 
equation (2.22) for the mass coefficient of absorption 
W(s)ds, 
N° | , 
(s, + Avy 


Av 6b + Av) (2.23) 


Ay 


+ Av\?| 
Av < 14h N° Bas = | exp (* »W(s,)ds,. (2.24) 


I'he expression for the mass coefficient of emission may be derived similarly. 
Where \,* ts the total number of metal atoms in state 2 per unit mass of photo- 
spheric matter and where we neglect stimulated downward transitions we find for 


Av 6b + Av) = Ay (2.25) 


— 


1 

+ Av) 

au 
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2.4. Assumed photospheric conditions.—\in deriving the values of the source 
functions we shall neglect the radiatively stimulated emissions. ‘Thus in thermo- 
dynamical equilibrium the source function, given under these conditions by 
Kirchhoff’s law, will be represented by Wien’s law—a sufficiently accurate 
approximation to Planck’s law for the optical part of the spectrum and photo 
spheric temperatures. Further, we shall suppose that in the photosphere the 
electrons and atoms have a Maxwellian distribution of velocities corresponding 
to a common kinetic temperature 7',, that there is a Boltzmann distribution 
among the stationary states corresponding to an excitation temperature 7’, 
for hydrogen and 7',* for the metal, and that the degree of ionization may be 
represented by Saha’s expression with an ionization temperature 7, for hydrogen 
and 7,* for the metal. 

As is customary we shall assume that the continuous absorption and emission 
arise from photo-ionization of the negative hydrogen ion, and its re-formation by 
the recombination of a free electron with the neutral hydrogen atom in its ground 
state. Free—free transitions, which contribute some to per cent of the continuous 
absorption and emission for A<7000A, are for simplicity neglected, though, 
following Cillié (§), their effect could be readily included, 

If n’ be the total number of neutral hydrogen atoms per unit volume, the 
number, m,’, in the ground state is 

n, =n'q, 6(T,). (2. 
‘The number of negative hydrogen ions per unit volume ts 
n=n'dT.)p., 2h) 
where 
2(27m)** b'(T,) 
HT) KT) 
E, is the (negative) dissociation energy for the negative hydrogen ion. Formally 
the number of negative hydrogen ions in the ground state per unit volume is 


n, = ng, (T,), (2.30) 
though since there is only one state, and it single, g,=(7,)=1. Hence from 
equations (2.27), (2.28) and (2.30) 


6(T,) 


exp(£, k7,). 29) 


Assuming the photosphere to be composed overwheimingly of neutral hydrogen 
in its ground state, the number of hydrogen atoms in the ground state per unit mass 
of photospheric matter is 
N, = my, (2.32) 
where my 1s the mass in gm of the hydrogen atom. Hence from equations (2.31) 
and (2.32) the number of negative hydrogen ions per gm of photospheric matter ts 
1 qy T,) 
V, and N,’ are the quantities to be inserted in equations (2.11) and (2.12) for th: 
mass coefficients of continuous absorption and emission 
Again, if n* be the total number per unit volume of neutral atoms of the metal 
under consideration, the numbers in states 1 and 2 will be 


{| —(E,* — 
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where E,* and E,* are the (neyative) energies of these states. If n*’ is the 


number of once ionized metal atoms in all states of excitation per unit volume we 
have 


(2.35) 
where #*(7;,*) is given by the expression corresponding to equation (2.29). 
If the ratio of abundance of hydrogen to the specific metal be A*~10*, then 
=(n+n')/A*. 
But under photospheric conditions n<n' and n* €n*’, so that 
(2.36) 
from equation (2.27). Hence 


~ myA® $*(T,*)p., (2.37) 
myA* b*(T.*) (7;*)p., (2.38) 


and these are the quantities to be inserted in equations (2.23) to (2.26) for the 
mass coefficients of line absorption and emission. 


2.5. Source function for the continuum.—From equations (2.11), (2.12), 
(2.32) and (2.33) we find 


Jo_ ta,(V) dV 


In this expression “ ) is given by equation (2.7) with the kinetic temperature 
Ty, (7) by equation (2.29), the ratio of atomic constants by equation (2.9), 
and from equation (2.5) dV/dv=h/mV. With these values we find 


jo 2hv® (RT,\ mV* 


In an enclosure containing hydrogen in thermodynamic equilibrium where 
T, = T,= T, equation (2.39) reduces, when account is taken of equation (2.5), to 
hv/kT), (2.40) 
that is to Wien’s law as it should. 


2.6. Source function for the line.—From equations (2.23), (2.25), (2.37) and 
(2.38) we find for Av > 6b 


K(vy_"+Av) By 


W(s,)ds, /3,(° W(s,)ds, 
From Bohr’s frequency condition —E,* =hyv,,°, while the relation of the 


Einstein coefficients, defined for energy density, is 
An n° 8rh( 


Hence Bs 


(2-41) 


ae 
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Similarly for Av < 1-4b, we have from equations (2.24) and (2.26) 


+ Av) 2h(v,.°) 
+ Av) 


In an enclosure in thermodynamic equilibrium containing the metal atoms 
only, 7,* = T and the expressions W(s) ds are given by (2.13). Hence in equation 
(2.41) we may write the numerator of the expression in square brackets : 

« — Av w + 6," ][ 5,7 + 


exp ( —hv,.°/kT,*) 


nm (Av)? +(6,+6,)?’ 
while the denominator similarly has an identical value. Hence in thermodynamic 


equilibrium equation (2.41) reduces, as it should, to Wien’s law. Again, in 
equation (2.42) the numerator of the expression in square brackets becomes 


where the integration is effected as in equation (2.21) and the last step follows 
since Av<1-4b. ‘The denominator of the expression in square brackets has 
similarly an identical value, so that in thermodynamical equilibrium equation 
(2.42) also reduces to Wien’s law. 

A further test for the accuracy of our expressions for the mass coefficients 
of absorption and emission in the continuum and line spectrum may readily be 
applied. Suppose we had a sample of our photospheric matter with hydrogen 
and the fraction 1/A* of the metal in an enclosure in thermodynamical equi- 
librium at some temperature comparable with that of the photosphere. Then in 
the frequency of the line the mass coefficient of emission will be jg + j(v4,° + Av), 
and the mass coefficient of absorption xy +,(v,.°+ Av). Inserting the values of 
these coefficients as given above and remembering that in thermodynamical 
equilibrium 7, = 7,= 7,* = T,;=7,* = T, we find after a little reduction that 

+ Av) _ 
ko + Av) 
reducing to Wien’s law as it should. 

2.7. The ratio of absorption coefficients.—For Av 2 6b we find from equations 
(2.11), (2.23), (2-33) and (2.37), since v= v4,", 


me LAT) 


d(T,) is given by equation (2.29) and ¢*(7;,*) by 


I 2(2m)** b*'(T,*) on 


exp(—hv/kT), 
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where — E,* is the tonization potential of the neutral metal. Inserting these 
values we have 


Vor Av < 1°46 we find from equation (2.24) an expression identical with the nght- 
hand side of equation (2.43), except that the integral is replaced by 


+ Av\ 2) 
exp 4 » W(s,)ds,. (2.44) 
\ 7 ) 


3. Numerical integration of exact equation of transfer 


Multiplying through by the integrating factor 
R= Ryexp ¢ sec | (1 +)dt> 


the exact equation of transfer (2.3) may be written 


1 (0, 6) = | is) exp —secé | +) dt > sec dry. (3.1) 


Ke 


Vhis is an integral equation from which to find 2 I and 7 a8 functions of optical 
Ko 


depth 7, (in the frequency of the adjacent continuous spectrum) from the observed 
limb darkening in the line. Its numerical integration may be effected from the 
knowledge we have gained in the previous section of the dependence of these 
three functions on kinetic (7',), ionization (7, and 7,,*) and excitation (7, and T,,*) 
temperatures. In order to put these functional relations in a more tractable 
form, three simplifying assumptions may be made. 

The first assumption is that 7, = 7,= 7, where 7, is the radiation temperature. 
The justification for this assumption has been given elsewhere (14). If it ts 
legitimate we have from equations (2.39) and (2.5) 


2hy 
exp (—hv/kT,). (3.2) 


‘The resulting value of 7’, as a function of r, can be found in one of the usual ways 
from the numerical integration of equation (3.1) for the case 7 =o, that is from 
the darkening to the limb, J,(0, #) in the adjacent continuous spectrum. 

The second and more dubious assumption is that 7;,* (for the metal) = 7. 
his assumption is questionable first because the frequency of the ionizing 
radiation, at least for ionization from the ground state of the metal, is far removed 
from the frequency v for which 7\{v) is determined as a function of 7, from 
equation (3.2). Secondly, since the radiation field is crossed by the already 
formed absorption lines of the metal, the excitation temperature, 7.,*, must be less 
than 7, so that 7;* must be less than 7,. The effect on the degree of 
ionization may however be small (14). The consequence of this simplifying 
issumption is that equation (2.43) takes the form 


where for the negative hydrogen ion q, = 1 and for the neutral atom b'(7T,) = q,' 
since all the neutral hydrogen atoms are effectively in the ground state. 


Equation 


| 
me 
; 
; 
j 
ge 
tk 
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(3-3) applies to the case where Av >66; for Av <1-4b the equation remains the 
same except that the integral is replaced by that in (2.44). 

Such integrals containing the unknown functions W(s,) and W(s,) involve 
us in a certain difficulty. Ideally we should like to determine these functions of 
s and 7», but the Fraunhofer lines by themselves do not appear to yield sufficient 
information. Applying the mean value theorem 
W(s,)ds, 3 I 


I 


the last step following since W(s,) is a sistigtitiend function. In this expression 
s,’ is a particular value of s,. 

Our third assumption is that since W(s,) is a peaked distribution centred on 
s,=0 we may put s,’=o0. In thermodynamical equilibrium W(s,) is given by 
the expression (2.13), which is very peaked indeed, and in the photosphere will 
differ from this only in that not all the sub-levels will be filled ; it may however be 
expected still to be peaked. With this assumption equation (3.3) for » takes the form 


2 B, E an 


Av <1-4h = €Xp 4 (=) exp {(E, — E,*)/kT,}. 
(3-5) 


In these expressions the quantity in square brackets is independent both of 
r, and of Av. Apart from the unknown abundance ratio of the metal, which we 
wish to find, it involves only known atomic constants. In front of the square 
brackets stands the factor which depends primarily on Av and after stand the 
factors which depend upon 7’,, that is upon the optical depth 7). Of these 
latter factors the dominant one is the exponential function which is completely 
known, — E, being the dissociation potential of the negative hydrogen ion, — E,* 
the ionization potential of the metal while 7’, is known as a function of r, from 
equation (3.2). Of the remaining factors dependent upon optical depth, 6 may 
be calculated from equation (2.18) and 4, from the theory of collisional broaden- 
ing, since both the gas pressure, p, and 7, = 7, are known (14) as functions of 
optical depth, 7). We may therefore write equations (3.4) and (3.5) 


7 


C 
Av 2 6b h(t), (3.6) 


Avsigh Cyfto)exp[— (3-7) 
where f,(7,), fo(t») and g(7,) are known functions of r,, and C, and C, are constants 
to be determined. 

Finally and by virtue of our third assumption we may write equations (2.41) 
and (2.42) in the form 
4 Av) 2h(v,,")? 


h 0 3 
Av < = exp (3.9) 


Av >6b exp(—hv,,° kT,*), (3.8) 


Following a suggestion originally made by ten Bruggencate (3), the numerical 
solution of the integral equation (3.1) may now be made. By a numerical solution 
of this equation for =o we find j,/x9, therefore from equation (3.2) 7,, as a 
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function of r,. Then im a specific frequency in the line, namely v4. = 43° + Av, 
for the case Av 266, arbitrarily choose some value C,’ of the constant in 
equation (3.6). From this equation we then know 7 as an explicit function 
of 7), and we may therefore solve equation (3.1) to find the only unknown 
function of optical depth, that is j’(v,.° + Av,)/«’(vy,°+Av,). At some second 
frequency Av, the corresponding values of 7 as a function of optical depth are 
from equation (3.6) given by 


(7) = (=) (79): 


Now solve equation (3.1) for j’(v4_ + Av,),«'(v4y° + Avg). Then if we have chosen 
our constant C,’ correctly, since from equation (3.8) j(v,2° + Av)/«(v4," + Av) is 
independent of Av, the two thus found values of this function should be identical 
at each optical depth r,. If they are not another value C,” is selected and the 
whole calculation repeated, if necessary again and again, until the equality 
+ Avy) + Avg) 
Avy) + Avg) 
is satisfied at each optical depth. 

This procedure may be contrasted with that employed earlier (12) in the 
numerical integration of the Eddington—Milne equation. In that investigation 
7 was an unknown function of r,, and it was only as a result of very protracted trials 
that its functional form was finally derived. In the present method, on the other 
hand, only a single constant, C,, has to be found and a strictly iterative procedure 
may be followed. A further disadvantage of the earlier method, quite apart 
from the fact that it was applied to the inexact Eddington—Milne equation, was 
that the test for the correct functional form of » could only be applied after a 
tedious calculation of 


dw 
| 9) 


at each optical depth—a calculation which took longer than the preceding 
solution of the integral equation. In the present method, on the other hand, the 
test is immediate and involves no further calculation after the solution of the 
integral equation for the frequencies Av, and Av,. As far as can be judged without 
an actual trial the present method should therefore prove practicable and fairly 
rapid, 

From the solution of the integral equation two results follow. First from the 
value of the constant C,, which contains as unknown only the quantity A* as 
is seen from equations (3.6) and (3.4), we obtain a direct observational determin- 
ation of the abundance ratio of the metal to hydrogen. Secondly from the 
determination of j(v,.° + Av)/«(vy4°+Av) as a function of 7, we find from 
equation (3.8), and a calculated value of the ratio 5,/4,, the value of the excitation 
temperature, that is the ratio of populations in states 2 and 1 of the metal atom, 
as a function of 7). If in addition the integral equation can be solved for the 
case Av < 14) we obtain from the constant C, in equations (3.7) and (3.5) an 
independent determination of the abundance ratio, a direct determination of the 
excitation temperature from equation (3.9) and, combined with the already 
determined equation (3.8), an observed value of the ratio 5,/5,. The solution 
for the case Av <1-4b follows with obvious slight modifications the method 
already outlined for Av 266, but depends for its success upon an accurately 
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determined line core, an observationally difficult problem, and upon a reasonable 
change in surface brightness within Av <1-4b. If the observational difficulties 
can be overcome it is probably the better method for faint lines, and the more 
directly informative. 

The solution of the integral equation has been treated for the case of resonance 
lines, but it is not in principle limited to them. The only requirement for the 
application to any other line is that the population in the lower state of the line 
has been previously found relative to the ground state of the atom by the prior 
application of the method to an appropriate line. 


In conclusion I wish to express my thanks to Mr R. J. Bray for many helpful 
discussions on the problem of line formation. 


University Observatory, 
Oxford 
1955 March 16. 
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Summary 

Abt has discussed the possible influence of hyperfine structure on solar line 
profiles. Allowing for hyperfine structure, and with the aid of detailed assump- 
tions concerning physical conditions on the Sun, he reproduced the Utrecht 
profile of the weak manganese line A 5516°8. ‘This paper describes the inter- 
ferometric investigation of the moderately strong manganese lines AA 6021, 
6016, 6013, and the iron line A6027. The method of circular channels 
was used, but with the line under examination at the centre of the hetero- 
chromatic fringe system. A solar absorption coefficient is deduced from 
the iron line profile, and this is used to compute the profiles of the manganese 
lines allowing for hyperfine structure. The only major assumption necessary 
concerning solar conditions is that the lines are all formed in the same level 
of the solar atmosphere. The computed and observed profiles agree within 
the limits of experimental error, and are substantially different from the 
manganese line profiles computed on the assumption of no hyperfine structure. 


Introduction.—Hyperfine structure in atomic spectra is attributed to the 
interaction of the orbital electrons of the atom with nuclear magnetic and electric 
moments. ‘lhe largest effects are observed when the nucleus possesses a large 
magnetic dipole moment. Abt (1) has discussed the possible influence of 
hyperfine structure on certain solar lines of manganese, vanadium and some other 
elements. He found some correlation between the half-widths of the solar 
profiles and the total hyperfine splitting of the line as observed in the laboratory. 
Abt also found evidence that hyperfine structure influenced the curve of growth. 
More particularly, he computed the profile of the line Mn 5516-8 allowing for 
hyperfine structure (known from laboratory experiments) and assuming an 
excitation temperature, a kinetic temperature, a turbulence velocity and a damping 
constant for the solar atmosphere. He used a Milne~Eddington model of the 
solar atmosphere with pure scattering. After allowing for the instrumental 
pattern, he found it possible to reproduce the Utrecht profile of the line. A very 
large number of assumptions about the solar atmosphere are involved in this 
computation, and further, the line A 5516°8 is rather weak (Rowland intensity 0). 

A more convincing approach seems to be to compare profiles broadened by 
hyperfine structure with iron line profiles of similar intensity. No hyperfine 
structure has ever been detected in the iron spectrum. 

This paper describes the investigation of the profiles of three moderately 
strong manganese lines and an iron line with a view to establishing the influence 
of hyperfine structure on the manganese lines. Manganese possesses only one 
stable isotope whose nucleus has a spin of § and a magnetic dipole moment of 
3°47 nuclear magnetons, so that large hyperfine structures (~o-1 A) may be 
expected, 
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1. Experimental method ' 

1.1. The choice of lines.—The manganese lines chosen should possess large 
hyperfine structures which are known from laboratory observations or which may 
be reliably calculated. They should be reasonably intense, to allow an accurate 
enough determination of the profiles, and for the same reason they should be in 
a clear region of the spectrum. This latter consideration eliminates most of the 
violet manganese lines. It seems unlikely that the observed doubling at 
Mn 4030°8, 4033°1, 4034°5 is due to hyperfine structure, for the observed separa- 
tion is about twice the total hyperfine splitting measured in the laboratory. The 
most satisfactory lines seemed to be the triplet 


34°48(7S)4p* Pie. 1a 34°48(7S) 55° AA 6013°5, 6016-6, 6021°8, 
of Rowland intensity 6. 


The iron lines should be of similar intensity to the manganese lines and in 
the same spectral region. ‘They should be of pressure class a or 6 so that, 
presumably, they do not suffer from excessive collision broadening. ‘The only 
suitable iron line was Fe 6027'1, Rowland intensity 4. This line had the additional 
advantage that the excitation potential of the lower term is about the same as that 
of the manganese lines. 

1.2. The observational method.The Oxford solar spectroscope was used in 
conjunction with a Fabry-Perot etalon to achieve the necessary resolving power, 
‘The method of circular channels due to Treanor (§) was adopted, but with a new 
relative position of the channels and the line under observation. ‘The axis of the 
etalon was at such an angle to the axis of the light beam emerging from the 
spectroscope that the line being examined fell at the centre of the heterochromatic 
fringe system. ‘This centre was made bright by small adjustments to the angle of 
the etalon. For photometric studies this method avoids the difficulty of photo- 
metering along a narrow circular channel. Further, the absorption line is long 
enough to allow a long photometer slit, which helps to reduce the effect of plate 
grain. ‘This arrangement of the method of circular channels for line profile 
studies was suggested by Dr M. G. Adam. ‘Treanor used a field lens in the focal 
plan of the spectroscope. This was not employed, as the field of view under 
examination was always quite small, and the removal of the field lens eliminated 
an appreciable quantity of scattered light. ‘The focal lengths of the second 
collimator and camera lenses were 15 cm and 47 cm respectively. ‘The etalon 
plates were mounted with their wedge angles parallel to the spectroscope slit, 
and a small enough slit height was used so that the ghost images were completely 
clear of the central part of the primary fringe system. ‘The focal plane of the 
camera had an aperture 5 cm long by 08 cm high, so that four solar exposures 
could be taken on a single 2 in. x 4 in. plate, as well as the intensity marks. 
Plate 4 shows the appearance of the photographs. 

1°5 mm and 2-5 mm spacers were used and the reflectivity of the silver films 
on the etalon plates was about 93 per cent. ‘This gave theoretical resolving 
powers of 1°8 x 10° and 3-0 x 10° for the two spacers respectively, while photo- 
graphs of the cadmium line 46438 indicated actual resolving powers of about 
12x 10° and 20x 10°. The difference between theoretical and actual resolving 
powers is a result of plate imperfections and possibly the limitation of the number 
of reflected beams by the etalon aperture. 
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For the 1-5 mm spacer the spectroscope slit width was o-25 mm and for the 
2°5 mm spacer, o'1omm. ‘These were the largest slit widths allowable for the 
order-to-order scattered light to be kept sufficiently low. Ilford Rapid Process 
Panchromatic plates were used, and the exposure times were 2-5 min for the 
1°5 mm spacer and 8 min forthe 2-5 mm spacer. All exposures were of the centre 
of the solar disk. 

The plates were calibrated by taking a photograph with the etalon removed 
and a calibrated step wedge mounted immediately in front of the spectroscope 
slit. A similar photograph without the step wedge was taken to make possible 
a correction for the non-uniform illumination of the slit. ‘The exposure times for 
all photographs on the same plate were the same, a condition achieved by the use 
of filters and the adjustment of the spectroscope slit width for the calibration 
photographs. 

The plates were stored for 24 h before processing, to allow the initial fading 
of the latent image to take place. They were developed in Parkhurst developer. 

1.3. Photometry.—-Microphotometer (4b) tracings were made across the 
centre of the heterochromatic fringe system at right angles to the absorption line 
with both low (6°89 x ) and high (48-2 x ) magnifications. ‘The purpose of the 
former was to locate the precise position of the absorption line with respect to the 
centre of the heterochromatic system so that the dispersion could be calculated 
(see below). The width of the image of the photometer slit on the plate was 
0-057 mm and the height o-4 mm. 

The characteristic curve of each plate was determined by photometering along 
each step of the step-wedge photograph parallel to the dispersion and through 
absorption lines. By making use of points in the lines a large number of points 
on the characteristic curve were obtained. ‘The correction for non-uniform 
illumination of the spectroscope slit was derived from a tracing of the photograph 
without the step wedge, at right angles to the dispersion. By means of the 


characteristic curves thus obtained the photometer tracings of the lines were 
reduced to observed profiles. 


2. Treatment of the observed profiles 
2.1. Determination of the dispersion.—At the centre of the heterochromatic 
fringe system the dispersion in a monochromatic fringe is equal to the dispersion 
in the spectrum obtained in the focal plane of the camera lens without the etalon. 
The dispersion in a monochromatic interference fringe is given by the expression 
dx = Ax/f*, 
where f mm is the focal length of the camera lens, and x mm the distance from 


the centre of the system. Now the distance between heterochromatic and 
monochromatic centres is 
p=frw/r 


(see Treanor (5)), where w A/mm is the dispersion without the etalon. ‘Thus 
when x=p, w=dA/dx. If, however, x= p+ 2, where z is the distance between 
the heterochromatic centre and the line under consideration, the dispersion in 
the line is 
dA/dx = (p+ 2z)A/f* = w(1 + 2/p). 

p was calculated to be 81-8 tracing squares with the large magnification, while < 
varied from o to 7 tracing squares. w was measured on tracings of a spectrum 
without etalon and found to be 0-096 A/tracing square on the high-magnification 
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tracings, 0°673 A/square on low-magnification tracings, and o-58 A/mm on the 
plate. 2 was measured on the low-magnification tracings and dA/dx calculated 
for both low and high magnifications. The observational profiles were thus 
determined with the correct dispersion. 

2.2. The scattered light correction.—The general scattered light in the 
spectroscope is known to be about 1:5 per cent with a slit height of 19 mm. With 
the present slit height of 3 mm the general scattered light was therefore negligible. 

In the method of circular channels there is always a certain amount of light 
scattered from one heterochromatic order to another because the slit image of the 
spectroscope extends over more than one order of interference. ‘This occurs 
even with quite narrow slits, since the spectroscope is by no means a perfect 
image-forming device, and has a rather wide apparatus function. In order to 
estimate the proportion of scattered light, the intensity half-way between hetero- 
chromatic orders was measured. ‘There are two contributions to this intensity. 
First, even if the spectroscope were a perfect image-former, there would be an 
inter-order intensity due to the interferometer alone. Secondly, there is the 
light that is due to the imperfection of the spectroscope. It is only the latter with 
which we are concerned, for the former is not scattered from order to order. The 
purely interferometric inter-order intensity was calculated, assuming the Airy 
formula for the intensity distribution, but remembering the finite slit width, 
and subtracted from the total to leave the true scattered light. So long as the 
slit width is small compared with the distance between orders, the true scattered 
light intensity falls off inversely as the square of the fractional order distance 
from the neighbouring order (see Bray (2)). ‘Thus the scattered light at the line 


being measured is } of the inter-order scattered light arising from the spectroscope 
imperfections. ‘The measured values were :— 


Inter-order intensity mm spacer = per cent 
2°5 mm spacer = 19°7 
Pure interferometric (calculated) 1:5 mm spacer = 0-9 
spacer= 
Inter-order true scattered light 1-5 mm spacer = 10°! 
2°5 mm spacer 
Scattered light at measured line 1°5 mm spacer= 2°25 
mm spacer= 


The profiles corrected for scattered light yielded the values for central 
intensities and equivalent widths shown in ‘Table I. 

The agreement of the equivalent widths for the two spacers is evidence that 
the scattered light correction has at least the correct ratio for the two spacers. 
‘The disagreement of the central intensities suggests that the difference in resolving 
powers with the two spacers is significant for the profiles concerned. 

2.3. Averaging the results from the two spacers.—I\t is desirable to average the 
results from the two spacers. ‘To achieve this, each point on each 1°5 mm 
profile was multiplied by a factor representing the difference in resolving power 
between the two spacers. ‘This factor was not difficult to compute, as all the 
profiles were found to be quite closely Gaussian. ‘The function turning a 1-5 mm 
profile into a 2-5 mm profile must, therefore, also be Gaussian. ‘The function 
was found from the Fe 6027 and Mn 6021 profiles, for which the best data were 
available, and applied to the others. The average profile for each line was then 


4 
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calculated. Fig. 1 shows the observed points obtained for Fe 6027, together 
with the mean profile. The final mean central intensities and equivalent widths 
of the profiles, with the root mean square errors of the mean, were as follows: 


Central intensities Equivalent widths (A) 


Fe 6027 0°468 + o-012 0°065 + 0°0013 
Mn 6021 0412+ 0°012 0089 + 00015 
Mn 6016 0°437 + 0°022 0°086 + 0°0031 
Mn 6013 + 0°012 + 00018 


Taare I 


Line Central intensities Equivalent widths 
mm 2°5 mm mm 2°5 mm 
Fe 6027 0°494 0'065 0°059 
0°408 0°430 0068 
0°§33 0°433 0064 
0°530 
Mean 0°493 0064 
Mn 6021 0°470 0°463 
0°430 0°424 0'092 
0°370 o'og! 
Mn 6016 0°447 
0°507 0°374 0°082 0°097 
Mean 0°490 0°374 0°083 
Mn 6013 0'085 
0°53! 
0'502 
Mean 0°524 ° 


| 


©.7 


Fic. 1.-—The observed profile of Fe 6027. Circles represent points from 2°5 mm spacer exposures 
and crosses points from 1°5 mm spacer exposures, converted as described in the text, Section 2.3 
Ordinates are intensities as a fraction of the continuum, and abscissae are wave-length differences in 
angstroms. 
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2.4. Systematic errors.—There may well be a number of systematic errors 
in the profiles. The methods of computing the scattered light and of averaging 
the results from the two spacers may not be thoroughly sound. However, these 
errors are likely to be the same for all four lines, and since the object is a comparison 
of the profiles and not an absolute determination, systematic errors are unlikely to 
be of much importance. 


3. Computation of manganese profiles from the iron profile 

3-1. Deduction of the solar absorption coefficient from the Fe profile.—The 
Milne—Eddington model with pure scattering was selected. For this model 
the residual intensity r for observations at the centre of the solar disk is given by 

_ + hv — + $3) + Ev 3)" 

- 
where n=1+5,/ky, €=1/n, ky=mass coefficient of continuous absorption, 
5, = mass scattering coefficient, and c = coefficient of limb darkening. From the 
data of Moll, Burger and van der Bildt, c=0-962. ‘To solve the equation for n, 
the method used by Plaskett (4a) for a more complicated problem was adapted. 
The equation may be written in the following form : 

n 3)! 
The right-hand side of this equation was evaluated for £ =0-0, 0°2, 0°4,...1°0 
for values of r in the observed profile in the iron line. For each value of r a graph 
of € against 1/n was plotted. ‘The value of & for which £ = 1/n was then selected, 
and this value of 1/n was the solution of the equation. n — 1, which is proportional 
to the solar scattering coefficient, was thus determined as a function of wave-length. 
The curve possessed a very slight asymmetry, which may have been accidental. 
The two halves were averaged to produce a symmetrical scattering coefficient. 

3.2. The hyperfine structure of the manganese lines.Vhe manganese lines 
AA 6013, 6016, 6021 have not been thoroughly investigated in the laboratory, 
but the hyperfine structures of the terms from which they arise have been deduced 
by White and Ritschl (6) from measurements on other lines. ‘The structure of 
the upper term has also been determined by Fisher and Peck (3). ‘The application 
of selection and intensity rules then permits the evaluation of the structures of 
the lines AA 6013, 6016, 6021. There is no evidence of perturbations in any of the 
laboratory investigations, so that there is no reason to doubt that the intensity 
rules are obeyed. ‘The hyperfine structures thus evaluated are given in ‘Table I] 
and shown diagrammatically in Fig. 2. 

3-3 Computation of the manganese profiles.—'Vhe absorption coefhcient of 
each hyperfine structure component was given the same shape as that of the iron 
line but a magnitude proportional to the relative intensity of the component. 
In this way an absorption coefficient for the whole line was evaluated. ‘lhe 
absolute magnitude of the coefficient was adjusted so that when the coefficient 
was substituted into equation (1) the equivalent widths of the computed and 
observed profiles agreed. After two or three attempts a computed profile of 
the correct equivalent width was produced. In Fig. 3 the observed and computed 
profiles are displayed. ‘The slight dip on the long-wave side of Mn 6021 is 
probably the predicted line Fe 6021-82 ; this dip was not included in the equivalent 
width for normalization purposes. It is seen that the agreement between the 
observed and computed profiles is well within the random experimental error. 
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Taste Il 


Mn 6021 Mn 6016 Mn 6013 
Wave-length Relative Wave-length Relative Wave-length Relative 
difference intensity difference intensity difference intensity 
(mA) (mA) (mA) 

o'o o'7 16°7 100°0 
29°0 25°5 61°4 
29°5 16 33°5 100°0 36°0 32°7 
§2°5 47°5 24'5 47°5 12°7 
52°5 54°5 52°5 20°5 
64°0 100'0 62°0 19°4 61°0 28°6 
18 67°5 25°5 64°0 27°3 
71'S 19°2 75°5 19°1 

1's 72°5 70°5 10°6 79°0 
82°0 16°5 83°0 8-2 
O's 88-0 16°7 89°0 5°5 
94°5 50°4 93'5 10°6 94°5 2°3 
99'5 33°0 94°5 24°2 
102'0 10°43 95°5 24°5 
103°0 19°4 
| 
Mn 6021 
| 
OO97 A 
Mn 6016 


| ll Lis. 


-—-O095A 
Mn 6013 


FiG. 2.——The hyperfine structure pattern of the manganese lines. 
the relative intensity of the component 


The length of each line represents 
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Mn 6021 


/ Mn 6016 


Mn 6013 


-O2 -O.1 


Fic. 3.—The profiles of the manganese lines The continuous lines are the observed profiles, 
and the broken lines the computed profiles Ordinates are intensities as a fraction of the 
continuum, and abscissae ave wave-length differences in angstroms 
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Since the three lines have significantly different central intensities but similar 
equivalent widths, such agreement is unlikely to be accidental. Further, Fig. 4 
shows the profile that one of the manganese lines would have if there were no 
hyperfine structure—this is the iron line profile with the same equivalent width 
as the manganese lines—and it is clear that the manganese line profiles are being 
substantially modified by hyperfine structure. 


T 


o8 


0.7 


-0.2 -O 0.0 0.1 0.2 


Fic. 4 The profile of Mn60143 computed with hyperfine structure (continuous line) and 
without hyperfine structure (broken line) Ordinates are intensities as a fraction of the continuum, 
and abscissae are wave-length differences in angstroms 

The observed profiles of the manganese lines have thus been satisfactorily 
accounted for by the inclusion of hyperfine structure as a cause of broadening. 
It is emphasized that the only major assumption lying behind the above procedure 
for comparing the iron and manganese line profiles is that the manganese lines 
and the iron line are formed in the same level of the solar atmosphere. ‘This is 
likely to be approximately true, as the intensities and excitation potentials of the 
lines are similar. ‘The solar model is no doubt deficient, and for this reason the 
solar scattering coefficient deduced for the iron line probably has little meaning. 
‘The comparison of the profiles must, however, be insensitive to the model chosen. 
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CONTINUUM INTENSITIES IN PLANETARY NEBULAE 
M. Seaton 


(Received 1955 January 22)* 


Summary 
Thirty-four planetaries are considered for which continuum intensities 
have been measured, 7,, may be calculated from the [O 111] intensities and 


N, may be calculated from forbidden line intensities or from the surface 
brightness. 

Values of B,/H4 measured by Aller (B, = Balmer continuum intensity) are 
in reasonable agreement with theory. About half of the nebulae considered 
appear to come closer to Case A than to Case B (A ~ optically thin, B = optically 
thick in Lyman lines). 

Total continuum intensities are calculated for H 11 regions allowing for : 
(1) all H and He 11 recombination and free—free continua and (2) the H 2s —. 18 
two-quantum continuum, collisional deactivation of H 2s being taken into 
account, ‘The calculated intensities, relative to Hd, are compared with the 
measurements of Page for both sides of the Balmer limit. ‘The mean results 
are in fair accord but there is a considerable scatter probably due to obser- 
vational errors. 

Recently measured Balmer discontinuities (A) are in good agreement with 
theory for 6 nebulae, but for NGC 6826 A (obs) is significantly smaller than 
A (calc). This object is unusual in that it is surrounded by a faint outer 
envelope. ‘The small value of A (obs) is probably due to the conversion of 
Ly « into H 2s .. 1s quanta in an extensive H 1 region. When such processes 
are not important A (obs) may be used to improve the estimates of 7, and N, 


1. Introduction.—-Spectra of planetary nebulae reveal the presence of bright 
emission lines and of background continua. ‘The present paper attempts a detailed 
comparison between theory and observation for the continuum intensities. 
This provides a useful check on our understanding of the physical processes 
and on our knowledge of the physical conditions in these objects. 

The most complete published observations are those of Page (1,2), who 
studied the region from 5000 to 3000A. ‘The continuum observed for 
A> 3647 A was termed by Page the “visual continuum”, V,. For A<3647A 
the intensity is much greater due to the Balmer continuum of hydrogen. From 
the discontinuity at 3647 A we may estimate the intensity per angstrom unit of 
the Balmer continuum at the Balmer threshold, which in the present paper will 
be denotedt by B,. Page reported large systematic differences between theory 
and observation for the ratio B,/H4. If confirmed, such differences would imply 
that the theory used to calculated the hydrogen line intensities was incorrect. 
The question is reconsidered in the present paper using measurements of Aller 
4, §) for Ha. 


* Received in original form 1954 April 27 


+ The intensity of the Balmer continuum is frequently denoted by Bag but Page (a) denotes by 
Ba, the total intensity for A<— 3647 A. 
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The best test for the theory of the continuum intensities is provided by 
recent measurements of the Balmer discontinuities made by Minkowski (6), 
by Barbier and Andrillat (7) and by Aller (8). An apparent discrepancy between 
the measured discontinuity for NGC 7027 and that to be expected using the 
electron temperature calculated from the relative intensities of [O 111] lines has 
been pointed out by Aller (8). ‘The question is discussed in the present paper 
using improved [O 111] intensities and allowing for deactivation effects. 

The processes most likely to be of importance for the continua are (i) recom- 
bination and free-free transitions and (ii) two-quantum emission from certain 
metastable states. 

The H-recombination continuum *, for which the emissivity may be calculated 
accurately, may be expected to be the most important process of type (i). After 
hydrogen the next most abundant element is He; according to Wurm (g) the 
He/H abundance ratio is about 3. Recombination to highly excited states of 
He will contribute to the continuum in the observable spectral region. Exact 
calculation of this contribution would be difficult, but since the states concerned 
will be nearly hydrogenic we may estimate that the total intensity of the 
Het-recombination continuum will not exceed 4 of the H-recombination 
continuum, ‘This contribution will be neglected. In certain high-excitation 
objects the He tl-recombination continuum may be more important due to the 
greater attractive power of the He** field. In the present paper we consider the 
H-recombination and the Hetl-recombination continua and neglect all other 
recombination continua, ‘The intensity of the H--attachment continuum has 
been estimated by Greenstein and Page (10), who consider it unlikely to be of 
importance, 

The continuum resulting from 2s->1s transitions in atomic hydrogen by two- 
quantum emission, which has been considered by Kipper (11) and by Spitzer 
and Greenstein (12), may be expected to be the most important process of type 
(ii). This will be referred to as the 2q-continuum. Possible contributions 
from similar transitions in He and in He* will be neglected. ‘The semi-classical 
calculations of Purcell (13) for collisional deactivation of the 2s metastable state 
of hydrogen, together with recent estimates of the densities in planetary nebulae, 
suggest that the 2q-emissivity may depend critically on the electron density. 
A detailed quantum-mechanical study of the deactivation process has therefore 
been made (14), and the results obtained, which should be of a sufficient accuracy, 
are employed in the present paper. 

2. Summary of theory formulae 


2.1. Hydrogen line intensities.—Bohm and Aller (1§) have shown that the 
electron velocity distribution in planetary nebulae will be very close to Maxwellian. 
For a given electron temperature 7’, the hydrogenic continuum emissivities are 
readily calculated, but certain additional assumptions have to be made in order 
to calculate the line emissivities. It is to be expected that conditions will always 
lie somewhere between the two extreme cases considered by Baker and Menzel 
(16); Case A refers to a nebula optically thin in the Lyman lines and Case B to 
a nebula optically thick in the Lyman lines. For both A and B it is assumed that 
the stellar Lyman line radiation is absorbed in the stellar atmosphere and lower 
nebular layers. 


* The total intensity resulting from recombination of electrons to X * ions, together with free-free 
transitions of electrons in the fields of X * ions, will be referred to as the X-recombination continuum, 


+ - 
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The number of bare nuclei with charge Z,(H_', He’ ‘, etc.) per cm* is denoted 
by N, and the number of electrons by V,. The energy emitted in transitions 
n--n' by the corresponding hydregenic ions (H, He ', etc.) is then given by 


= 10 NiN, erg sec ', (1) 


where t,=10°7, and X,,=15-78/n"t,. Baker and Menzel (16) give the coeffi- 
cients 5, as functions of 7, (Cases A and B) and also the Kramers—Gaunt factors 
Ban’ 

In the present paper all calculations will be made in the first place assuming 
Case B for both H and Het. Conversion factors for various combinations of 
Cases A and B will then be given in Section 2.8. 

2.2. The He** abundance.—Denoting the «-particle density by N, and the 
proton density by V,, Wurm (9) gives* 

N,/N, = (2) 
where 
Il = t1, A 4686)/7(HB, A 4861) (3) 
and 
&=e*+b (He 11) b,(H). (4) 
The quantity € has been calculated by Wurm with the following results (Case B): 


t.=10°*T, é 
0°37 
0°40 
0°44 


It will be shown in Section 3 that ¢, lies between 1-0 and 2°0 in most planetaries, 
We therefore take £ = 0°42, giving 


(5) 


The hydrogen Balmer lines n-+2 are blended with the He tt lines 2n-»4. It is 
readily shown that this does not have any marked effect on the Balmer decrements. 
2.3. Hydrogenic recombination and free-free continua. emission intensities 
are calculated in a form which is most readily compared with the observations of 
Page (1, 2); denoting the emission intensities per angstrom unit by /, we put 


J, = 108, /HB. (6) 


The contribution to J, arising from recombination to state a of a hydrogenic ion 
is given by 
where 
(5 
(v,t,) exp | ~(0°435 + 1°439v) 
and v is the wave number in reciprocal microns (i.c. v= 10*/A with Ain A). For 
the Kramers—Gaunt factors we take = 0°872, g,, and put all 
other g,,“" equal to unity. The error arising from this approximation will be 
of order 5 per cent or less (Menzel and Pekeris (17)). 


* The author is indebted to Dr K. Wurm for informing him of this work in advance of 
publication, 
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The contribution to J, arising from free—free transitions in the fields of ions 
of charge Z; is 
N 
= ¥(v, t,) (9) 
with (17) 
y\ 
0-078 (7) (1+ 1°39f, v). (10) 
i 
We consider the wave-length range corresponding to 2°0<»<3:0; the 
Balmer limit is at v= 2°74. Within this range the H-recombination contribution 


to J, is 


and the He tl-recombination contribution is 


J (He tt) = +d yy (12) 
The functions /,'(H) (J ,(H) for Case B) and J,'(He 1) (J (He 11) for Case B and 
N,/Ny*}) are given in Tables I and Il. In preparing these tables the 
summations over n for n 2 10 were replaced by integrals which are readily evaluated. 


Ne 


Tasie | 


J, 
o's 


nN 
| 6 
° 


5000 A ' O'S, 3'0 4°0 
4545 2'2 O'36 14 2°3 3°2 4°9 
4167 2°4 O'2, 3°3 
3346 it 2°2 3°3 
3047 2°74 O'ls 2‘! 3°3 
3047 2°74" 17°5 15°9 15°2 15°7 
357! 16°7 15°7 1§°2 15°9 
3333 30 13°8 14°5 14°9 16°6 


Il 


\ J,’ (He 11) 


2°2 2°7 3°6 56 
2'4 3°5 4°4 
2°6 2'1 3°4 4°5 6°5 
2°74 o'4 3°3 4°5 6°38 
2°74 10°6 10°5 10°38 12°0 
10°2 10°4 10°38 12°2 
7's 8-8 9°9 10°7 12°7 


2.4. The H2s two-quantum continuum.—Spitzer and Greenstein (12) give 
the emissivity of the 2q-continuum proportional to a quantity XY which is less 
than or equal to unity. In calculating X they considered the following two 
processes by which H atoms could enter the 2s state ; 


3 
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Process «.—-Atoms enter the 2s state by capture or by cascade from higher 
discrete states. X is then defined by 
_ Number of 28-+1s transitions by the 2q process (13) 
~~ Number of recombinations on all levels with n > 2° 3 


Neglecting deactivation, Spitzer and Greenstein estimated X to be between 


0°30 and 0-35 and adopted the value 0°32. Allowing for deactivation X is then 
given by 


0°32 

T+ WN 
where WN, is the probability per second of collisionally induced 28->2p transitions 
and A,, the total radiative transition probability for the 2q process. According 
to Spitzer and Greenstein A,,=8-227 sec*. Using the results of Breit and 
Teller (18) for electron impact producing 28-2p transitions, Spitzer and 
Greenstein concluded that WN,/A,,, was small compared to unity. However, 
Purcell (13) has recently shown, using semi-classical methods, that collisions 
with protons will be more important, a result which is confirmed by quantum- 


mechanical calculations (14). Assuming equal densities of electrons and protons 
we obtain* 


W=5:3x10'cm' sec! for T,=1-0x 104 deg. K 
and 


W=47x10'cm' sec! for 7,=2 x deg. K. 


Consideration of the probable sources of error in the quantal method used leads 
to the conclusion that the results for W should be correct to within 20 per cent. 
Similar calculations have been carried out for 38->3p transitions in Na produced 
by electron impact, for which good agreement with experiment was obtained. 
Adopting the mean value W= 5-0 x 10°* cm® sec”! we obtain finally 
0°32 
10 (15) 
Process B.—Atoms may enter the 2s state by collisionally induced 2p->2s 
transitions. For this process to be important a sufficiently large number of 


atoms must enter the 2p state by absorption of Ly « quanta from the ground state, 
The parameter X is now defined by 


x Total number of 28->1s transitions by the 2q process 


"Total number of Ly « quanta escaping , (16) 


The value of X depends on the total amount of neutral hydrogen, and in con- 
sequence a precise calculation becomes much more difficult; we may therefore 
consider some limiting conditions. If Case A is applicable, then process f 
is certainly unimportant. Spitzer and Greenstein concluded more generally 
that process 8 was negligible for processes occurring in regions where most of the 
hydrogen was ionized. Allowing for the larger value of W adopted in the present 
paper it would still appear that this conclusion is justified. On the other hand, 
if the ionized part of the nebula is surrounded by a sufficiently extensive neutral 
hydrogen region, X defined by (16) may approach the value unity. 

It is by no means certain that most planetaries are surrounded by H1 regions 
(Wurm and Singer (19)), and for a number of objects there are some grounds for 

* These results are about one-half of the published results of Purcell. The difference is largely 


due to the fact that Purcell did not allow for the reduced mass. 
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believing that Case A may be applicable (Wurm (9), and Section 4 below). For 
a first survey of the problem of the sources of the continuous emissions we will 
therefore neglect process £ ; its importance may be reconsidered if other processes 
prove inadequate to account for the observed continuum intensities. 

The 2q contribution to J,,J,'(2q), is given in Table III calculated for Case B 
with X = 0°32. 


Taste Ill 


J 


1'O, 1°2 1°3 
2°2 13 1's 1°6 1°7 1°8 
2°6 23 2°6 2°38 3°0 
2°74 34 3°5 3°7 
30 3°3 3°6 3°58 40 
30 3°7 471 


2.5. Corrections for space reddening.—The calculated relative intensities 
of the Balmer lines are comparatively insensitive to 7, and to the particular model 
assumed. ‘To a first approximation it may be assumed that the true relative 
Balmer emission intensities are the same for all planetary nebulae and that 
variations in the observed relative intensities are due to differences in the amount 
of space reddening. Berman (20) obtained a correction for space reddening by 
plotting the observed Balmer intensities against an effective path-length through 
the obscuring galactic stratum. In this way he obtained average decrements 
corrected for reddening which were found to be in good agreement with the 
theory. 

Since there is considerable patchiness in the distribution of the obscuring 
clouds, we shall adopt the procedure of determining the amount of reddening 
from a comparison of theory and observation for the Balmer decrements. 

In the wave-length range of interest the relation between reddening and 
absorption may be represented with sufficient accuracy by 


log (14(A)/To(A)) = — v9) (17) 
(Whitford (21), Divan (22)), where /,(A) is the true emission intensity, /,(A) 
the observed intensity and ¢ the reddening constant. We assume the optical 
properties of all obscuring clouds to be the same (Divan (22)) and take vy = 0°34 
(Whitford (21), Berger, Divan and Fringant (23)). Calculations will be made 
in the first place assuming Case B, for which the calculated value of ¢ will be 
denoted by ¢’; if Case A is applicable c=c’ + 0°21, the additive constant being 
practically independent of 7,. 

2.6. Determination of electron temperatures.—Electron temperatures may be 
calculated from the relative intensities of the [O 111] lines (Menzel, Aller and Hebb 
(24), Seaton (2§)). We obtain 

t, = 1°43/ {log (1(4959 + 5007)/1(4363)) 0°942 + log (1 + 0°0438 x 10° *N,t, **)}, 
(18) 
where t, = 10°*7, (7, in deg. K). In using this formula we allow for absorption 
and calibration errors by taking /(4959 + 5007) relative to Hf = 100 and / (4363) 
relative to Hy=50 (25). ‘The deactivation term involving N, is in all cases 
small but is none the less significant, 


+ 
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2.7. Determination of electron densities.—-Electron densities in planetary 
nebulae may be determined by two independent methods. ‘The first utilizes 
the relative intensities of forbidden lines (Seaton (2§)) and the second utilizes 
the absolute emission rate in the hydrogen spectrum as estimated from the 
absolute surface brightness and the nebular dimensions (Menzel and Aller (3)). 
The forbidden line method, which will be referred to as method (a), can at present 
only be applied to a limited number of objects, since in most cases insufficient 
observational material is available. For our present purpose the densities 
obtained by method (a) are to be preferred, since, if local condensations occur, 
this method gives an average local density effective for deactivation processes, 
It also has the additional advantage that it is not necessary to make any assumptions 
about the nebular distances or nebular geometry and that the results are practically 
independent of the theory assumed for the hydrogen lines. 

Considerable improvement in the measurement of the surface brightnesses 
of nebulae have recently been made, using photoelectric techniques, by Liller 
and Aller (26) and* by Liller (see Aller (8)). However, the published densities 
tend to be too small due to neglect of space absorption. Denoting the observed 
average surface brightness of Hf by Sg, and using equation (17), we obtain for 


the true surface brightness 10'7”S,. For a spherically symmetric nebula we 
may readily obtain from (1) 


log (| N, =log Sg+ 1°72¢ + 2 log A, + 11°47 +2D—logr, (19) 


where Sg is in erg cm * sec”', A, is the angular radius in seconds of are adopted 
for the measurement of Sg (the measured flux being proportional to A,*Sz), r is the 
distance of the nebula in parsecs, NV, and N, are the local electron and proton 
densities in cm™*, which may be functions of A, the angular distance from the 
centre measured in seconds of arc and finally 


D = log (20) 
Tasie IV 

o's 

o'176 

2°0 

4'0 0°467 


Table IV gives D as a function of 7, for Case B, denoted by D’; if Case A is 
applicable D= D’ +0121, the additive constant being practically independent 
of T,. 

In employing (19) we will put NV, = N,, an approximation which will certainly 
be justified if the level of ionization is low. ‘The most serious departures from 
this relation may be expected when nearly all the He is doubly ionized. We then 
have N.>N,+2N, with N,~4N,. However, in this case we should also allow 
for the Hew contribution to Hf; in (19) N,N, should be replaced by 
NAN, +0°6N,). It readily follows that the values of NV, calculated assuming 
N, = N, should be multiplied by 1-15 and hence that the error is small. 

In employing this expression we follow Liller and Aller in assuming the 
density to be uniform in the spherical shell defined by A, < A < A, and to be zero 


* The results obtained by Liller (8) have been adopted in preference to those of Liller and 


Aller (26) since they are apparently the more accurate. ‘The differences between the results of these 
two determinations are, however, in all cases small, 


20° 
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outside this domain (note that A, may equal zero, in which case we have a uniform 
sphere). With this assumption we obtain 
log N, = 5°97 + 4 log Sg— 4 log (A,r) + + —(A,/A,)*]. (21) 
The use of this relation, together with the values of Sg, A,, A, and r adopted by 
Liller and Aller or by Liller, will be referred to as method (b). It is readily shown 
that the densities calculated from (19) assuming Case A will be greater than those 
calculated assuming Case B by a factor of 2:0. 
For nebulae for which there is insufficient observational material to use either 
method (a) or method (b) we may obtain approximate density estimates using the 


photographic magnitudes m,. Denoting the photographic magnitude in 
monochromatic Hf light by mg we may put 


mg = m,, + Am( HB). (22) 
The quantity Am(H) may be expected to depend on the nebular spectrum and in 
particular to vary with the ratio /(4959 + 5007)/HB8; however, the values of 
Am(H 8) calculated by the method of Menzel and Aller (3) show this dependence 
to be comparatively small. Since mg= —2°5 log (A,*Sg)+const., we have 
log Sg +2 log A, +0°4m, = —0-4Am(HB) + constant. (23) 
Using the results of Liller and Aller or of Liller for Sg together with the values of 
m,, quoted by Vorontsov-Velyaminov (27) it was verified that the left-hand side 
of (23) showed no systematic dependence on log(1(4959+5007)/Hf). ‘The 
variation of Am(Hf) was therefore neglected. Adopting the mean for the 
15 nebulae studied by Liller and Aller and by Liller we obtain 
log Sg +2 log A, + 0-4m, = 4°35 + 0°31, (24) 
where +0°31 is the root mean square deviation. In view of the larger errors 
involved we adopt the simple approximation of considering the nebulae as uniform 
spheres with diameters in seconds d=2A,. ‘This gives 
log = 8-60 o-2m,, } logd— }logr + D + o'86c. (25) 
The use of (25), together with the values of m,, d and r listed by Vorontsov- 
Velyaminov (27) and assuming Case B, will be referred to as method (c). 
A similar equation quoted in a previous paper (28) was obtained using a con- 
version from m, to Sg essentially similar to that of Menzel and Aller (3). ‘This 
equation gives densities greater than those obtained from (25) by a factor of 2. 

It may be noted that the errors involved in using (25) may in certain instances 
be considerably greater than are indicated by the r.m.s. deviation in (24). The 
most extreme case is NGC 2440, for which method (a) gives NV, = 4x 10‘ cm * 
and method (c) gives V,=4x10?cm™*. ‘The relative exposures (RE) of Curtis 
(29) give a rough measure of the surface brightness. A plot of log(RE) against 
m,+5logd suggests that m,+5logd is greatly overestimated for this object. 
It also suggests that m, + 5 logd is somewhat overestimated for I 418, and this is 
confirmed by the photoelectric measures. 

2.8. Summary of theory relations.—Fortunately the ratios 6,(B)6,(A) are 
very insensitive to 7, and in consequence conversion factors from Case B to 
Case A may be obtained which are practically independent of 7. The final 
expression for the continuum intensities, J, = 10°/,/H5, may be written in the 
form 


J, = od m1) + (2q), (26) 
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where J,'(H), J,((Hem) and J,'(2q) are given in Tables I, II and III and 
where o and + are given in Table V for various combinations of Cases A and B. 
According to Wurm (9g) the combination of Case B for H and Case A for 
Heil is very improbable. 
Taste V 
Case, H B A A 
Case, He 11 B B A 
1°00 1°46 1°46 
T 1°00 1°46 
3. The observational material 


3.1. The data of Table VI.—We consider 34 planetaries for which continuum 
intensity measurements have been made and for which estimates of 7\, and 
AX. may be obtained. ‘Table VI gives the following data: 

(i) The catalogue numbers (right ascensions and declinations are given by 
Aller (§) and by Vorontsov-Velyaminov (27)). 

(ii) The space reddening constants c’ calculated from a comparison of the 
observed relative intensities of HB, Hy, H3, H, and H,, with the theoretical 
relative intensities for Case B with t,= 1-5. We do not use H, and Hy, since they 
are frequently blended with other features. 

(iii) The ratio Il = /(4686)/1(4861) (Il may be assumed to be small when no 
value is given). 

(iv) Values of t,=10°*7., calculated from (18). 

(v) Values of N, calculated by methods (a), (b) or (c). 

(vi) The Balmer continuum intensity relative to Hé according to Aller 
(3, 4,5). ‘The observed ratio, (10°B,/H8),,,,, must be corrected for reddening 
using (17) and for blending with Heit using (1), (5) and (7). ‘The quantity 
tabulated is 


log (10°B,, Hd) = log (10°B,/ Hd),,,, + 0-300" + log ( 


(27) 


‘The last term on the right of (27) is small (less than 0-10 in all cases). 

3.2. Accuracy of the calculated electron densities.-Comparison between the 
results obtained using methods (a) and (b) has been made elsewhere (30). It is 
seen that fair accord is obtained for NGC 6572 and for NGC 7027. For 1418 
method (b) gives a density greater than that obtained from method (a) by a factor 
of at least 3:1; however, if in place of equation (21) we use equation (19) together 
with the density distribution function of Wilson and Aller (31), this factor is 
reduced to 1°6 (30). For NGC7662 method (a) appears to give a density 
appreciably greater than that obtained from method (b), a result which might be 
explained in terms of local condensations. ‘This will be discussed further in 
Section 6.3. 

The r.m.s. deviation between the values of log N, calculated from methods 
(b) and (c) is o-2 but bigger errors may occur in certain cases. For one nebula 
of low density, NGC 40, methods (a) and (c) have been shown to be in fair accord 
(28). 

4. Balmer continuum intensities relative to Hd,—Equation (7) gives 

log (10*B,. Hd) = 0-880 — (0° 190 — log dy. (28) 


In Fig. 1 this function is plotted against 7, and compared with the observational 
results summarized in ‘lable VI. 
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Taste VI 
Nebuls, 
NGC, ete. 10~* N, 
An. 25™ 
I 351 
II 2003 
1535 


Notes 
(See p. 239) 


> 


0°29 (c) whe (1) 
0°35 (c) (1) 
(c) (1) 
0°26 (b) B . (1, 2) 
(b) A 

(a) (2, 3) 
(b) B 

5 (b) A 

0°28 (c) (1) 
017 (c) (1) 

4 (a) rq (1, 4) 
017 (c) (1) 
0°34 (b) B ' (1, 2) 
(b) A 

3:1 (b) B (2, 3) 
6 (b)A 

{c) (1) 
(b) B (2, 3) 

5 (a) (2, 3, 5) 
1°9 (b) B 

A 


1 418 


wn 


Il 2165 
J goo 
2440 
Il 4568 
Il 4593 


6210 


Il 4634 
6543 


6572 


Nw 
wn 


An. 47™ 
6741 

IT 4846 
6778 

6790 

6803 


vee (1) 
0°5 (c) ial (1) 
(c) ‘ (1, 6) 
0°06 (c) (1) 
10 (c) . (1, 6) 
o's (b) B , (1, 2) 
ro(b)A 
(c) (1) 
(b) B (1, 2, 7) 
os (b) A 
(b) B ° (2, 3) 
r2(b)A 
2°0 (c) (1, 6) 
1°4 (c) (1) 
(c) (1) 
io (b) B ° (2, 3, 8) 
20 (b)A 
(c) (1) 
4 (a) ' (2, 3, 9) 
4 (b) B 
A 
1°9 (c) (1, 6) 
(c) (1) 
(ce) (1) 

ree (1) 
5 (a) : (2, 3, 10) 
ri (b) B 
22(b)A 


6807 
6818 


6826 


63833 
6884 
6886 
7909 


O 


7026 
7027 


Il 5117 

An, 21% 
II 

An. 22" 
7662 


°° 
1° 
1° 
1 
I 
I 
1 
I 
I 
I 
I 
I 
I 
I 
2 
I 
I 
°° 
1° 
1° 
1° 
1° 
1° 
1° 


Notes 

(a) Densities obtained from relative intensities of forbidden lines (2g) 

(b) Densities obtained from photoelectric measurements of surface brightness (26, 8), 
Cases A and B being considered separately. 

(c) Densities obtained from the surface brightness as estimated from the photographic 
magnitudes. Case B is assumed. ‘The temperatures are not corrected for deactivation. 
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Fic. 1.—log (10"B,/H8) against T,. Theory curves are given for Case A and Case B using (28), 
The observed values are plotted from Table V1 


The main conclusion we draw from Fig. 1 is that there is no evidence for 
systematic differences between theory and observation. The present results 
differ from those of Page (2) because we have used: (i) higher electron tempera- 
tures (25); (ii) the Balmer intensities of Aller in place of those of Page (for a 
number of objects it appears that Page may have measured B, too small; see 
Section 5 below) ; (iii) a different method of correcting for reddening (the method 
of Page would be expected to give a bigger scatter); (iv) a correction for blending 
with Het. 

Of the 25 objects considered in Fig. 1, 11 come close to Case A, 9 come close 
to Case B and 5 are in poor agreement with theory. It is doubtful, however, if the 
observations are sufficiently accurate to decide which case is applicable for 
individual objects. ‘The mean value of the quantity c’ is 0-44 for the objects 
apparently coming close to Case A and is 0°24 for the objects apparently coming 
close to Case B. It would be somewhat surprising to find real differences in 
reddening between these two groups.* For Case A objects the correct space 
reddening constant is c=c’ + 0°21 (Section 2.5); we might therefore expect the 

* The possibility cannot, however, be ruled out. Johnson (32) considers that there is evidence 


for a real correlation between certain physical and spatio-kinematic characteristics of planetary 
nebulae 


(1) Line intensities, Aller (gs); (2) T, from [O 111) intensities adopted by Liller and 
Aller (26); (3) Line intensities, Aller (4); (4) 7, adopted by Seaton (ag); (5) N, calculated 
from the data of Liller and Aller (8, 26) differs from that given by Liller and Aller; (6) 
Stellar nebula with d uncertain (d= 2” adopted); (7) 7, for E side, continuum for W side (g) ; 
(8) 7, calculated from the data of Liller and Aller (26) differs slightly from that given by 
Liller and Aller; (9) The value of c’ is rather greater than that adopted by Aller (37), but 
Aller states that his value may be underestimated; (10) J (4959) 5007)/HB from MacRae and 
Stock (39). 
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quantity c’ to be smaller for the objects coming closer to Case A. It appears 
probable that some objects apparently come close to Case A due to c’ being 
overestimated, and that others apparently come close to Case B due to c’ being 
underestimated. However, if it be assumed that there are no systematic errors 
in the observations, it would still appear that Case A and Case B nebulae are about 
equally numerous. 

For Cases A and B the calculated values of log (10*B,, H4) differ by 0-17 and the 
calculated reddening corrections differ by 0-06. It is seen that it should, in 
principle, be possible to discriminate between Case A and Case B objects using this 
ratio. A better method might be to measure the ratio of B, to the intensities of 
the higher Balmer lines, as has been done by Greenstein (33) for the Orion nebula. 
The difference between log (10°B,/H,,) for Cases A and B is 0-10 and the reddening 
correction is very small. 

5. The total continuum intensities of Page.—\n the present section we compare 
the total continuum intensities* measured by Page (2) relative to Hd with the 
theory intensities obtained from (26). We consider Cases A and B for H but 
Case B only for Het. Page gives the maximum measured intensity shortward 
of the Balmer limit and the maximum measured intensity in the visual continuum. 
The maximum beyond the limit would be expected to occur at the limit and we 
therefore correct for reddening using the Balmer limit wave-length and compare 
with the theory results for the short-wave side of the limit, denoted by J 444,-. 
Page states that V,, (per A) is approximately constant in the range 3647 <A < 5000 A. 
‘Tables I, II, II] and VI show this to be in fair accord with theoretical prediction. 
We compare the maximum value of 10°V’,/H4 as given by Page with the values 
calculated for A= 4167 A (v= 24). 

We define the Balmer discontinuity as 


Vol. 115 


A = log (29) 


where /,,,. is the total continuum intensity on the long-wave side of the Balmer 
limit. Since V, is approximately constant we have 


Aslog [I s64;-/ Vc]. (39) 


The detailed results are given in ‘Table VII. ‘Tables | and V show that, 
considering H-recombination alone, J4,,,- should always be greater than 15 
for both Case A and Case B and for all values of 7,. Of the 21 measured values 
given in Table VII, 11 give Jyg4,-< 15. It is also seen that several of the nebulae 
with low measured values of J5.,,- have measured values of A close to zero, 
These results are difficult to reconcile with the theory. However, a comparison 
of Aller’s results for 10°B,/H8 with Page’s results for 10°(/5.,.-— V.) Hd suggests 
that the observational error in Page’s results is rather large and that this probably 
accounts for the apparent discrepancy with the theory. ‘The fact that the mean 
observed value of J5.,,- comes nearer to the theory value for Case B than to that 
for Case A is probably without significance. 

The comparison between theory and observation for V’, suggests that the 
processes considered are adequate to account for the main features of the observed 
visual continua but in view of the rather large scatter a more detailed discussion 
scarcely seems justified. 


* Aller (3, 4, §) does not give the total intensities 
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6. The Balmer discontinuity 

6.1. Theory relations.—The theory expression for the Balmer discontinuity 
may be written 

i, (H) + (Heit) + (X/0°32)ty 50 2q)) 
wheres denotes an intensity on the short-wave side andi, on the long-wave side 
of the Balmer limit. The value of X is discussed in Section 2.4 and the other 
functions occurring in (31) are given in Table VIII. The factor o-55II is obtained 
assuming that the same Case (A or B) applies for both H and Hew. If other 
combinations apply this factor should be multiplied by 7/o (Table V). 

6.2. Comparison of calculated and observed Balmer discontinuities.-Observed 
and calculated Balmer discontinuities are compared} in Table IX. ‘The results 
for NGC 6543, 6572, 6826, 700g and 7662 are those obtained by Barbier and 
Andrillat (7), who state that the internal precision for a single spectrum is of 
order 0-05 ; the number of spectra obtained are given in parentheses. ‘The result 
quoted for 1 418 is that obtained by Minkowski (6) and the result for NGC 7027 
is deduced from the statement of Aller (8) that the Balmer discontinuity is 
consistent with the theory of Spitzer and Greenstein for 7',= 10* deg. K. 


Taste VIII 


t. i_(He 11) i,(H) i,(He 11) 
o's 10°21 0°292 0°56 1°69 
7°51 4°43 0°65 1°38 
5°90 3°65 0°495 rig 
4°36 3°17 1°04 
1°6 4°14 o°631 0°95 
3°62 2°51 0°677 0°33 
2'0 3°22 2°30 0°707 0°78 


‘The units adopted are such that ** 


i NeNy is equal to the emission rate in erg cm~* 
sec” per reciprocal micron 


Taste IX 


NGC, Observed Calculated A 
etc A (a) (b) B (b) A N,-+0 
1 418 0°48 + 0°07 0°53 o'61 0°67 0'47 
6543 0°98 (4) 1°02 o'72 
6572 0°79 (1) 0°37 0°75 0°83 0°58 
6826 (4) 0°74 0°66 
7027 0°70 0°68 0°68 0°76 
7062 (2) 0°65 0°56 
Mean 0°76 0°575 


Note that NGC 6826 is excluded in calculating the means 


The Balmer discontinuities have been calculated assuming process 8 of Section 
2.4 not to be operative, i.e, using (15) for the calculation of X. Calculations have 
been made using the densities obtained by method (a) and by method (b) assuming 
both Cases Band A. ‘To illustrate the effect of the deactivation terms, calculations 
have also been carried out for the limiting case of N,->o. 


t A simular comparison has been made in a previous paper (34) using slightly less accurate 
values of 7, and N, 


{Barbier and Andrillat state that 


yoog was observd low on the horizon; Dr D 
Barbier informs the author that the observed value of A is of lower accuracy due to the large 
correction for differential atmospheric absorption. 


“ 
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It is seen that, with the exception of NGC 6826, the observations are in good 
agreement with the theory employed. It has previously been suggested (34) 
that the low value of A observed for NGC 6826 is due to the operation of process f. 
Long exposures obtained by Duncan (35, 36) show this nebula to be unusual 
in that it is surrounded by a faint and extensive outer envelope. It seems to be 
very probable that this outer envelope is an H1 region shining by the light of the 
H 2s->1s two-quantum radiation. An attempt to confirm this conclusion by 
obtaining the spectrum of the outer envelope would appear to be well worth while. 

6.3. Use of the Balmer discontinuity to obtain more precise temperatures and 
densities. —'The intensity ratios of forbidden lines and also the observed surface 
brightness give relations beiween 7, and N,. Provided that process 8 may be 
neglected, another such relation is provided by the measured Balmer discontinuity. 


su NZ 


SB. O/N 


4 


relations for NGC 7027 


Fic. 3 eo N, relations obtained from the (O U14\ ratio and from the Balmer discontimaty 
for NGC 6543. The observed Balmer discontinuity is 4 0°9% with a probable error of about 
40°03. To illustrate the sensitivity of the calculated density to errors in the measurement of 4 
curves are given for Ao-93 and 41-0 
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The relations between 7, and N, for NGC 7027 are shown in Fig. 2. The 
curves marked Su, Nu, Ot and O/N are those given in a previous paper (25), 
the curve marked O 111 is that obtained from (18) together with the intensities 
adopted in Table VI, and that marked A is obtained from (31) with A=o-7o. 
We also give the relation obtained from the surface brightness, assuming Case B, 
which is marked S.B. It is seen that all these relations are in fair accord, although 
the St curve indicates a density which is rather low. It is possible that a 
discrepancy for 5S 11 might arise from the fact that the atomic data are rather less 
accurate (25), but it seems more probable that the differences result from variations 
in the intensities across the disk, Aller (37) has recently given some intensities 
for a selected point in the brightest part of NGC 7027 which are somewhat different 
from those adopted for Fig. 2. Using these we obtain densities from S 11 which 
are appreciably greater than those obtained from O11 or N 11. 

The Balmer discontinuity is particularly sensitive to N, when 7\, is small 
and N, is not too large. This is illustrated in Fig. 3 for NGC 6543, for which 
apparently reliable values of 7, and of N, may be obtained from the |O111] intensity 
ratio and the Balmer discontinuity. It should be noted that the errors in 7, 
due to errors in the [O 111] ratio are much smaller when 7’, is small than when it is 
large. 

Improved estimates of the average local temperatures and densities, effective 
for excitation and deactivation processes, have been obtained from a combination 
of the data on forbidden line intensities and on the Balmer discontinuities. ‘These 
are given in ‘Table X. For 1418 and NGC 7662 the Balmer discontinuity is 
found to be in better agreement with the densities previously calculated from the 
forbidden lines than those calculated from the surface brightness, and in particular 
for NGC 7662 the measured discontinuity appears to confirm that the effective 
local density cannot be much less than that previously obtained. It appears that 
the average local density* exceeds, by a factor of at least 2, the root mean square 
density for the total volume obtained from the surface brightness. 


Taste X 
NGC, t,=10~T, 10 
et 
1 418 
6543 2°4 
6572 1°25 + 
2°5 
7027 4°5 
7662 5 


6.4. Variations of \ across the disk.—For NGC 7662 Page (1) obtained a 
distribution across the disk for V,, which was more uniform than that obtained for 
B.. This is in satisfactory agreement with what would be expected from 
deactivation effects. For 1418 Minkowski (6) states that A decreases inside the 
bright ring. It is evident from ‘Table IX that only a small decrease could result 
from a decrease in deactivation effects. A decrease in A might be explained in 
terms of an increase in 7\, but from a study of the factors determining 7\, Aller (38) 
has concluded that large variations in 1 418 areimprobable. It would be interesting 
to check this conclusion by measuring the variation of the [O 111] ratio across the 


* A typical average local density would be N, | be dv with S21. 
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disk and at the same time to obtain more definite numerical results for the 
variation of A. 
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Summary 


In Part I it is shown that the preference for near-commensurability of 
mean motions demonstrated in a previous paper extends to the small satellites 
of Jupiter, including the retrograde ones, and to the retrograde satellite of 
Saturn, implying that stability of near-commensurable configurations is the 
reason for such a preference. In Part II it is proved that if, at a certain epoch, 
a system of n gravitating point-masses has each radius vector from the 
(assumed stationary) centre of mass of the system perpendicular to every 
velocity vector (hereinafter called a “‘ mirror configuration’’), then the 
behaviour of each of the point-masses under the internal gravitational forces 
of the system after the epoch will be a mirror image of its behaviour prior to 
the epoch. It is further shown that, if a mirror configuration of the system 
exists at two separate epochs, then the orbit of each point-mass is periodic 
The authors argue that such periodic orbits are the more stable (under the 
action of external forces) the shorter the interval of time between mirror 
configurations, and they show that the frequent occurrence of mirror con- 
figurations between any two point-rmasses requires that the mean motions 
of the two point-masses be nearly commensurable. In Part III, various 
near-commensurable pairs of orbits in the solar system are shown to behave 
according to the above arguments. The relationship of the present work 
to the more restricted “‘ symmetry theorem ”’ of Griffin is discussed. 


Part | 


1. Introduction.—I\n a previous paper (1), hereinafter referred to as Paper I, 
the present authors have shown that, in the solar system, among the planetary 
and satellite systems, the occurrence of commensurability between pairs of mean 
motions is more frequent than in a chance distribution. It was suggested that 
the observed distribution was due either to a property of the mechanism of forma- 
tion of the solar system, or to an inherent stability of commensurable con- 
figurations. 


More recent work has led to the adoption of the latter hypothesis. ‘The 
authors were led to investigate this latter alternative by a consideration of bodies 
in the solar system not included in Paper I, namely Jupiter's satellites VI to XII 
and Saturn IX. 


2. The outer satellites of Jupiter and Saturn.—Using the nomenclature of 
Paper I, we have 


(2.01) 


where m,, m, are the mean motions of the two bodies, and A, and A, are integers 
(which, in Paper 1, were less than or equal to 7). 
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The Sun disturbs Jupiter’s satellites considerably, and the orbits of the 
retrograde satellites J VIII, JIX and J XI are not even approximately elliptic. 
Yet they have approximately equal mean motions, and if we consider the Sun 
as a satellite of Jupiter, the ratio of its mean motion to those of the retrograde 
ones is within the range of 1 to 7 previously considered. As is seen from ‘Table I, 
these satellites are nearly commensurable in mean motion with the Sun, the 
ratio being 1 to6. Indeed, the Sun’s mean motion about Jupiter is almost exactly 
one-sixth of the mean of the mean motions of J VIII, J 1X and J XI. 


Tasie I 
Ne Ne A, 
A A = 
Sun 
O'17055 I 6 0°003 BS 
Sun 
JIX 0°171 96 I 6 0°005 249 
Sun 
1X! 52 1 6 2006 
Sun and mean of 
J VIII, J IX, J XI 735 $9 
Sun 
0°057 84 17 0°000 
Sun 
JVI 0°060 02 I 17 20 
Sun 
1X 0°05867 I 17 0°000 15 
Sun and mean of 
JVI, JX 5883 I 17 0'000 O1 
Sun 
64 I 7 0'002 78 
lapetus 
Phoebe 1 7 0°00! 25 


Note : 0'003 68 


! I 
5 0°02381 


It is convenient to consider the other satellites of Jupiter here. The Sun's 
mean motion about Jupiter is nearly 17 of the mean motion of J XII, while the 
satellites J VI, J VII and JX form a group whose mean motions are nearly 
commensurable with the Sun’s at 1/17, and once again their average mean motion 
is almost exactly commensurable with that of the Sun. Finally, Phoebe (Saturn 
IX), disturbed by lapetus, pursues a retrograde orbit whose mean motion is 


commensurable with that satellite at 1/7. In Table | the mean motions are 


derived from the Handbook of the British Astronomical Association for 195%, 
except for those of J XII, lapetus, Phoebe and the Sun. ‘The last three were 
obtained from the Connaissance des Temps for 1954, while the mean motion of 
j XII was from Hilton and Brady’s orbit No. IV (2). 

These results provide additional evidence for the conclusions of Paper I, 
and suggest that the reason for the preference for near-commensurability is that 
such configurations are relatively stable. 


For, if the retrograde satellites of 
Jupiter and Saturn are captured bodies, it is difficult to understand why they 


should all have mean motions approximately commensurable with other bodies, 
unless they have reached positions of relative stability. 


| 
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Part Il 

3. The mirror theorem.—\n this section it is proved that, if n point-masses 
are acted upon by their mutual gravitational forces only, and at a certain epoch 
each radius vector from the (assumed stationary) centre of mass of the system 
is perpendicular to every velocity vector, then the orbit of each mass after that 
epoch is a mirror image of its orbit prior to that epoch. 

Consider an isolated system of n point-masses, centre of mass G, and let their 
masses be m,,i=1,2,3,...m. Suppose that at a certain time the radius vector 
from G to m, be R,, etc., and the velocity vector of m,be dR,/dt, etc. Then the 
equation of motion of m, relative to G is given, in gravitational units, by 

~4 R.-R 

The geometrical shapes and sizes of the orbits arising from this configuration, 
and their relative geometrical disposition, are completely specified by the masses 
m, together with magnitudes and relative configuration of the radii vectores R, 
and velocity vectors dR,/dt, i.e. by 

m, for all values of 7, 
R,.R, for all values of i and /, ae 
(dR,/dt). R, for all values of 7 and 
and (dR,/dt) .(dR,/dt) for all values of 7 and j. | 
Consider now the geometrical shapes, sizes and relative geometrical disposition 
of the orbits arising from the same configuration of radii vectores and velocity 
vectors, but with the new independent variable r= —?t. These orbits are 
completely specified by 
m, for all values of 7, 
R,. R, for all values of i and j, 
(dR, dr). R, for all values of ¢ and j, (3-03) 
and (dR,/dr) .(dR,/dr) for all values of i and j; 


the equation of motion (3.01) being invariant under the transformation of ¢ to r. 
The geometrical shapes, sizes and relative disposition of the orbits arising from 
the configuration in r are identical with those arising from the same configuration 
in ¢ if conditions (3.03) are identical with conditions (3.02). Unless certain 
symmetry relationships hold between the masses m,, conditions (3.03) are identical 
with conditions (3.02) if, and only if 


(dR, dr) .(dR, dr) =(dR, dt) .(dR, dt) for all values of 7 and j, (3.04) 
and (dR, dr). R,=(dR, dt). R, for all values of i and /. (3.05) 
But, from the definition of 7, (3.04) holds, and 
(dR, dr). R,;= —(dR, dt). R, for all values of i and j. (3.06) 
Condition (3.06) is equivalent to condition (3.05) if, and only if 
(dR, dr). R,=0=(dR, dt). R, for all values of ¢ and j. (3.07) 
Thus, if condition (3.07) holds at a certain epoch ¢ = r = 0, the motion of the system 
of n bodies after this epoch is geometrically similar to the motion of the system 
priortotheepoch. Ata mirror configuration the following conditions arise from 
the proot : 
(i) For all values of i, the geometrical shape and size of the orbit of the ith 
body after the mirror configuration is equivalent to the geometrical shape and size 


| 
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“Ss 
a 


No. 3, 1955 commensurable mean motions in the solar system 299 


of the orbit of the ith body prior to the mirror configuration, but described in the 
Opposite sense. 

(ii) For all values of i and j, the relative geometrical disposition of the ith 
and jth orbits after the mirror configuration is equivalent to the relative geometrical 
disposition of these orbits prior to the mirror configuration. 

(iii) For all values of 7, the ith orbit is continuous at the mirror configuration. 
‘These three conditions serve to define uniquely the nature of the “ reflection” 
of the orbits. 

4. The two possible cases of a mirror configuration.-Consider two of the n 
point-masses. Let NM,A, be the plane containing G and the instantaneous 
velocity vector of m,, and NM,A, be the plane containing G and the instantaneous 
velocity vector of m,, where M,, M, are the projections, on the celestial sphere 
centred on G, of the radii vectores Gm, and Gm, respectively, and A,, A, are 
the projections, on the same celestial sphere, of the instantaneous directions of 
motion of m, and m, respectively. ‘Then, if dR, dt. R, =o, 

M.A, = M,A, = M,A, = M,A, = 90 . 
‘Thus, either M, and M, are poles of the great circle A,A,, or A, and A, are poles 
of the great circle M,M,. ‘Thus, if longitudes are measured along the plane 
of the appropriate orbit, from the common node N, for a mirror configuration 
either 

(i) longitude of M,=o° or 180’, longitude of A,= +90 and longitude of 
M,=o0 or 180°, longitude of A,= +90, or 

(ii) longitude of M,= + 90", longitude of A,=o or 180° and longitude of 
M,= +90, longitude of A,=o or 180. 


Since for any pair of point-masses, cither (i) or (11) must hold, for a mirror 
configuration of the m point-masses, it is necessary either (i) that all the radu 
vectores R, are parallel or anti-parallel, and all the velocity vectors dR,/dt are 
perpendicular to all the R,, or (i1) that all the velocity vectors dR, dt are parallel 
or anti-parallel, and all the radii vectores R, are perpendicular to all the dR, dt. 

For any two of the n point-masses, the above cases correspond to a conjunction 
or opposition (with respect to G) with each body at either pericentron or apocentron. 
In case (i), the conjunction or opposition must occur at the common nodes of 
the two orbits, while in case (ii) the conjunction or opposition must occur go 
from the common node. 

‘To distinguish more clearly between these two cases, we discuss the nature 
of the “ reflection” of the orbits for each case separately. 

The equation of motion (3.01) may be written 

(x, 


where x, y, = are a rectangular set of axes with origin G, the centre of mass 
of the system of n point-masses. 
Case (i): — Suppose that at f= 7 =0, each m, is on the line y=0, and is 
moving in a plane parallel to yGz. Suppose also that 


Vir Zs ete. (4.02) 
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is a solution of (4.01). Then 
(4,=%,5 T= —t) ete. (4.03) 
is also a solution of (4.01), which is unchanged by the substitution. Hence 
(4.02) and (4.03) represent possible orbits. ‘The initial conditions of the set 
of orbits (4.02) are 
(x,, %,; dx,/dt, dy, dt,dz,jdt; t)=(a,0,0; 0,g,h,; 0) ete. (4.04) 
By the transformation (4.03), these initial conditions for the set of orbits 
represented by (4.03) become 
(x dydr,dzjdr; tr); =(ay0,0; h,; 0) etc. (4.05) 
(4.04) and (4.05) are identical, and therefore (4.02) and (4.03) represent the same 
orbit. Hence the orbit of m, for t -o is the mirror image of the orbit for to 
in the planes yGx and zGx. 
Case (ii): Suppose that at f= 7 =0, each m, is in the plane xGy and moving 
perpendicular to it. Suppose also that 
Vn ete. (4.06) 
is a solution of (4.01). ‘Then 
(x,=%5 T= ete. (4.07) 
is also a solution of (4.01), which is unchanged by the substitution. Hence both 
(4.06) and (4.07) represent possible orbits. ‘The initial conditions of the set 
of orbits (4.06) are 
(x, dxdt, dy, dt,dzjdt; t)=(a,,b,,0; 0,0,f/,; 0) ete. (4.08) 
By the transformation (4.07), these mitial conditions for the set of orbits 
represented by (4.07) become 


(X,, Vn dx,/dr, dy,/dr, dz,/dr; 7) =(a,,b,,0; 0,0,f,; 0) ete. (4.09) 


(4.08) and (4.09) are identical, and therefore represent the same orbit. Hence the 
orbit of m, for to is the mirror image of its orbit for t-—o in the plane xGy. 

5. The periodicity theorem.—-We define periodicity as follows:—In a system 
of  point-masses moving under their mutual gravitational forces only, their 
orbits are said to be periodic if, at periodic intervals of time, the same relative 
configuration of radii vectores and velocity vectors occurs with no change in scale. 

Hence the periodicity theorem :—If n point-masses are moving under their 
mutual gravitational forces only, their orbits are periodic if, at two separate 
epochs, a mirror configuration occurs, 

For, if the mirror condition (3.07) holds at t= —¢, and at =o, then it holds 
at t= +4), the orbits after ¢=o being the mirror images of the orbits prior to 
t=o. Hence the orbits are periodic, with period 2t,. 

6. Discussion of stability.—The above section gives no information on the 
stability of such orbits against perturbations by external forces. The manner 
in which stability is treated here makes no pretence to being exhaustive or rigorous, 
but the present authors believe that it throws some light on the possible importance 
of the periodic orbits discussed above, and leads to conclusions which are 
subsequently confirmed by appeal to observation. 

‘There are various accepted definitions of stability of an orbit. We will use 
Poincaré’s conditions for stability (4), namely : ‘ 

(i) the heliocentric distance of any planet cannot increase or decrease without 
limit, 
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(il) the system repeatedly passes through the configuration it had at time 
ty, Say a times /,, ty, etc., 

(ili) close encounters of any pairs of planets are ruled out. 

(The conditions for satellite orbits are analogous.) 

Condition (iii) implies that a body occupying a stable orbit will not stray fat 
radially from it. 

The discussion in Section 5 does not of itself satisfy Poincaré’s conditions. 
Although the orbits are periodic if the mirror condition is satisfied at two separate 
epochs, they might be unstable in the sense that a small change in their configura- 
tion would lead to non-periodic orbits and rapid departures from those orbits 
so far pursued. 

Let us consider the classical case of a system consisting of two massive bodies 
of unequal masses M, and M, pursuing circular orbits about their common 
centre of gravity, with a small mass m moving about the larger massive body M,, 
in the same plane as, and being perturbed by the other massive body, ™,,. 
Initially, let !,, M, and m satisfy the mirror condition, at time ft, say. ‘Chen m 
will move in a perturbed orbit and as a consequence its line of apses will rotate. 
In general, one synodic period of M, and m later, the velocity vector of m will 
be no longer perpendicular to the radius vector of m. But if the initial distances 
of mand M, have been chosen correctly, M,’s perturbation of m may be such that 
the time of revolution of m from one pericentron to some future pericentron is 
equal to one or more synodic periods of M, and m. ‘Then this new configuration, 
at f, say, satisfies the mirror condition again. ‘Therefore the behaviour of m 
after ¢, is a mirror image in the conjunction line M,mM, of its behaviour before ¢, 
‘Thus m will move in a disturbed periodic orbit, the perturbations having a period 
of 2(t;~t)). It seems likely that if mM, is such, initially, that the mirror 
configuration after t, does not arise to “reverse "’ the perturbations within a 
sufficiently short time, the disturbances might build up sufficiently to cause m 
to be captured by M,, or to collide with it or M,. Again, if the second mirror 
configuration does not arise sufficiently quickly after ¢,, the strongly perturbed 
periodic orbit of m might take it close to M,, when the perturbation of some 
fourth body might cause it to be captured by M, or sent to infinity. 

The periodic orbits studied by G. H. Darwin (3) are of interest, and may be 
said to be examples of the mirror theorem. Darwin considers the classical 
coplanar three-body problem, and confines his attention to simple periodic orbits 
which are periodic after one synodic period only, although the possibility of more 
complex periodic orbits with periods which are multiples of the synodic period 
is mentioned. In addition to isolating certain periodic orbits which correspond 
to the orbits discussed in this section, Darwin examines the stability of such 
periodic orbits. His results (in particular those relating to his “‘ Family A of 
planets ’’) to some extent corroborate the suggestions made in this section.* 

Hence the reasonable hypothesis is made that in this restricted three-body 
case, the periodic orbit of m will be the more stable against external perturbations 
the shorter the time-interval between the occurrence of mirror configurations. 

7. The effect of near-commensurability of mean motions of two bodies.-One 
major effect of near-commensurability is the occurrence of large inequalities, 
since the amplitudes of long-period terms involve the factor 

(Ayn, Ayn)". 
* We are indebted to the Referee for drawing our attention to this fact 
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Let us consider, after Newcomb (5), the mean motions n, and n, of two bodies 
orbiting about a primary in the same plane, and so related that 
An, — A\n,= 0, (7.01) 
A, and A, being integers, and n,<n,. Then A, revolutions of the first body 
will require the same time as A, revolutions of the second body, and at the end of 
this period they will have returned to their original positions, having been in 
conjunction (A, — A,) times at the same number of equidistant points. Calling 
these points conjunction points we put 
v= A, Ay 
But if the commensurability (7.01) is not exact, so that 
An, Ayn, 
then the mean conjunction points will revolve about the origin with velocity k 
given by 
Ayn, Ayn, 
(7.03) 
k approaching zero as « approaches zero. 
In the cases which we shall discuss in detail in this paper, 4, = A, +1, and 
therefore v= 1, so that we are concerned with only one mean conjunction point 
revolving uniformly with velocity 


k= Ayn, ~ Ayny. (7.04) 


By adjusting the mean motions, & can be given any value. 
It was seen in Section 5 that, if the mirror condition is satisfied twice, the 
orbits are periodic, and the reasonable hypothesis was made that the periodic 


orbit of m would be the more stable the shorter the time-interval between the 
two satisfactions of the mirror condition. ‘The shortest possible time between 
mirror configurations is obtained if one synodic period of m and M, equals the 
time taken by m’s pericentron to move from one conjunction to the next. But, 
by (7.04), if the mean motions are not exactly commensurable, the line of conjunc- 
tion revolves with a uniform velocity k, anda mirror configuration would occur 
at every conjunction if the mean motions of the bodies were such that 


hk=w 


where @ is the rate at which m’s pericentron revolves. 

We can now relax the conditions of the classical system considered hitherto, 
and consider two bodies moving in non-coplanar, non-circular orbits about a 
primary, and perturbing each other. We discuss two major cases : 

(a) If the inclinations are small, so that the orbits are approximately coplanar, 
and one eccentricity is much larger than the other, a first approximation to the 
frequent occurrence of mirror configurations is for the mean motions of the two 
bodies to be so nearly commensurable that the rate of revolution of the line of 
apses of the orbit of greater eccentricity equals the rate of rotation of the con- 
junction line of the two bodies. A perfect mirror configuration would occur 


only when We, OF 180, 


where w, and w, are the longitudes of the pericentrons. Hence a libration will 
ensue, whose period is the synodic period of the pericentrons. 

(+) If the inclinations are large, and the nodes of the orbits precess about the 
fundamental plane of the system at different rates, then by Section 4, a mirror 
configuration can occur only if conjunctions take place at either (i) the common 
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node of the two orbits, or (i1) go from the common node, with the bodies on the 
apse lines of their orbits in each case. In this case, frequent occurrences of 
mirror configurations can occur if the mean motions are such that conjunctions 
always take place according to (i) or (ii). 

It is of interest to make an estimate of the relative importance of (a) and (6) 
inan actual case. If the two orbits are coplanar, and one is circular while the other 
has a small eccentricity e, then it can be shown that, if a conjunction (or opposition) 
occurs at a small angie £ from the apse-line of the eccentric orbit, then the angle 
between the velocity vectors at conjunction (or opposition) is of the order of e€. 
On the other hand, if the two orbits are circular, and inclined at a small angle 7 to 
each other, then, if a conjunction (or opposition) (measured along one of the 
orbits) occurs at a small angle £ from the common node of the two orbits (or at an 
angle go + & from the common node), then the angle between the radii vectores (or 
between the velocity vectors) is of the order of £ sini. If the two orbits have small 
inclinations 7, and 2, (i, <i,) to the fundamental plane of the system about which 
their nodes precess, then the mutual inclination of the two orbits will vary between 
the limits (7, and (i, ‘Therefore, in a semi-quantitative way, we would 
expect that if 

sini, 
case (a) would be the more important, while if 


case (6) would be the more important, 
In Part II] we inspect certain systems in the solar system with outstanding 
commensurabilities, and see how far they support the above arguments. 


Parr III 

8. Saturn's system 

(a) Hyperion and Titan.-From ‘Table of Paper |, (ng/ny) | — 0-000 §7. 
The orbits are very nearly coplanar, and the eccentricities of ‘Titan and Hyperion 
are 0°02g and 0-104 respectively. Hyperion is very strongly disturbed by ‘Titan. 
On the basis of the arguments of Part II, we would expect that the mean motions 
of the two bodies would be such that the rate of revolution of the conjunction 
point would be the same as that of Hyperion’s perisaturnium. From (7.04), 
with n, = 22 °577009 per day (7), the point of conjunction retrogresses at a rate 
of 18°-66(9) per year. The rate of retrogression of the perisaturnium of 
Hyperion is 18°°66(3) per year (7). It 1s, indeed, well known (9) that Hyperion 
and ‘Titan come into conjunction with each other only at or near the aposaturnium 
of the former. More accurate elements derived since Newcomb’s time 
(6, 7,8) confirm this result. 

If ‘Titan had a circular orbit, a mirror configuration would occur at every 
conjunction and opposition. However, due to the eccentricity of Titan's orbit, 
a perfect mirror configuration would occur only when ‘Titan was also at its 
aposaturnium or perisaturnium. At other conjunctions, the satisfaction of the 
mirror condition is only approximate, and hence the reversal of perturbations 
at such conjunctions or oppositions is only approximate. A libration is therefore 
to be expected, whose period is the synodic period of the apse lines of Titan and 
Hyperion, i.e. 18°75 years (7), and a major libration of this period, of amplitude 
14° (7) is present in the motion of Hyperion. ‘The system of Hyperion and 
‘Titan is thus in excellent agreement with the arguments of Part II. 
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(6) Dione and Enceladus. —'Vhe orbits are equally inclined to the fundamental 
plane (the plane of Saturn’s equator), and retrogress at different rates. However, 
the inclinations are very small(~1'). The eccentricities of Dione and Enceladus 
are 0°0022 and 00044 respectively. Enceladus is the less massive body, and 
suffers the greater perturbations. 

For this system, a first approximation to the frequent occurrence of mirror 
configurations can be found by considering the orbits to be coplanar, and the 
orbit of Dione to be circular. ‘Then, for a mirror configuration to occur at every 
conjunction, the rate of motion of the perisaturnium of Enceladus should equal 
the rate of motion of the point of conjunction. Now, from ‘Table II of Paper I, 
(ny ny) = 4 +0°000643. Therefore, from (7.04), with ny = 262 -731 99 per day 
(7), the point of conjunction advances by 123°-41 per year. ‘The rate of advance 
of the perisaturnium of Enceladus, wy = 123°-43 per year (7). ‘Thus we have 
the relation noted by H. Struve (10) that 

2Np wy 
within the accuracy of the derived elements. Struve also noted that conjunctions 
of Enceladus and Dione occur always near the perisaturnium of Enceladus, 
about which point they oscillate. 

A second approximation to the actual problem is to consider the orbits to be 
still coplanar, but to consider the eccentricity of Dione’s orbit. ‘hen mirror 
configurations at conjunction will be only approximate unless the conjunction 
occurs when Dione is also at its perisaturnium or aposaturnium. A libration 
is to be expected whose period is the synodic period of the perisaturnia of Enceladus 
and Dione. Now w, = 30°°74 per year (6), so that 


within the accuracy of the derived elements (14). We would thus expect a 
libration in the longitudes of Enceladus and Dione of period 3°88 years. 
According to the Connaissance des Temps (7), a libration of amplitude 11°24 and 
period 3°89 years exists in the longitude of Enceladus, and a libration of amplitude 
~ 1’ of the same period in the longitude of Dione. 

A third approximation to the actual problem is to consider the inclinations 
of the orbits to the fundamental plane of the system. Exact mirror configurations 
can occur at conjunction only if the conjunction takes place at the common node 
or go from the common node of the two orbits. Owing to the fact that the two 
inclinations are sensibly the same, conjunctions which occur when the nodes of 
the two orbits on the fundamental plane are together will satisfy the mirror 
condition exactly, whatever the longitude of the conjunction point, since at this 
time the two orbital planes coincide. Now (7) 

= 328 — 152-71 
and = 276 — 31 
where @,, 0, are the longitudes of the nodes of Enceladus and Dione respectively, 
and ¢ is measured in Julian years from the epoch 1889 April o-s. 

.. the synodic period of the nodes = 2-96 years = B, say. 
But the synodic period of the apse lines = 3°88 years = A, say. 
or, within the accuracy of the derived elements (14), 
4B = 3A. 
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A libration in the longitudes of Enceladus and Dione of period 34 = 11-6 years 
is to be expected. According to the Connaissance des Temps, there is a libration 
in the longitude of Enceladus of amplitude 20’-0 and period 12-2 vears, while 
the longitude of Dione has a libration of the same period, and amplitude ~ 2’. 
According to the Berliner Astronomisches Jahrbuch (6), the amplitudes are 
14’ and 1° respectively, and the period is 11-1 years. It is seen that, within the 
accuracy of the elements, the libration period = 34, as predicted, 

As a further check, it was found that, according to the first set of clements 
(7), 4 conjunction of the nodes and a conjunction of the lines of apses occurred 
at t= —5°5 years, thus giving an exact mirror configuration at conjunction at 
this epoch and at all epochs separated from it by multiples of 11°6 years. The 
agreement is less satisfactory for the second set of elements (6), but it should be 
noted that, although the rates of retrogression of the nodes may be reasonably 
well determined, owing to the small inclinations the position of either node at a 
particular epoch is poorly determined. 

We may conclude that, within the limits of accuracy of the derived elements, 
the system of Dione and Enceladus conforms to the pattern of the arguments of 
Part IL. 

(c) Mimas and Tethys.—'The following elements have been given (7): 

Mean Motion Inchin. Eecen. Mass/primary 
Mimas 381-9045 36" (16 440 000)! 
Tethys 190 °697 95 0 ( 921 500)"' 
In this case, therefore, we have sini, -¢, and we would expect the mirror 
configurations to be related primarily to the position of the common node of the 
two orbits. 

For two orbits whose inclinations to the fundamental plane of the system are 
#, and i, (i, -1,), and whose nodes on the fundamental plane of the system have 
longitudes 4, and 4, respectively, it can be shown that the longitude of the common 
node 4, is given by 0, =O, +6, 
where tan ¢ (8.01) 

tani, .cott, — cos —4,) 
#, being measured along the fundamental plane of the system. 

As a first approximation to the actual situation of Mimas and ‘Tethys, we will 
consider that ‘Tethys moves undisturbed by Mimas, that the two orbits are 
equally inclined to the fundamental plane, and that both orbits have zero 
eccentricity. ‘Then, from (8.01), 


- 


If, therefore, conjunctions of Mimas and Tethys follow case (ii) of Section 4, 
the longitude of conjunction must always be (04 —4,) 2. Lf Oy, 0, are the rates 
of revolution of the nodes of Mimas and ‘Tethys on the fundamental plane, then, 
for this condition to hold, it is necessary that 


Ny = (Oy + 2 


But this is precisely the situation for these two bodies, as noted by H. Struve 
(11). Struve also noted that, if the mean longitudes (i.e. ignoring libration 


(6, — 
2 
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terms) of the satellites are ly, and /,, then 


4ly — 2ly — — Oy = 0. (3.03) 
Struve was thus led to enunciate the theorem that, ignoring libration, conjunctions 
of Mimas and Tethys occur always at a longitude which is the mean of the 
longitudes of the two bodies. Within the framework of the approximation 
considered here, this statement is equivalent to the statement that, ignoring 
libration, conjunctions occur always go from the common node of the two orbits. 
The mirror condition would, therefore, be satisfied at every conjunction. 
In passing, it should be noted that the improvement of the determination of the 
elements of the two orbits with time has rendered (8.02) more exact, viz. 


Tisserand des T. (7) B.A... (6) 
+3 +0254 

A second approximation to the actual problem is obtained by considering 
both orbits to be circular, but of different inclinations to the fundamental plane. 
Then the mean conjunction point is only exactly go’ from the common node if 
4,,~ 180 +4,. Allowing for the different inclinations, we would expect a 
libration in mean longitude of both bodies of a period equal to an integral number 
of synodic periods of the nodes. A libration of amplitude ~ 45° and period 
~ 70 years is present in the longitude of Mimas, and a small libration of the same 
period in the longitude of Tethys. For the satisfaction of the mirror condition 
at each phase of zero libration, on the present approximation, the period of 
libration, P, must be such that 

(Oy 67) P 2 = 360 m, 
or P=720 m, (Oy — Oy), (8.04) 
where m, is an integer. At the same time, the phase of the libration must be such 
that zero libration corresponds to an opposition of the nodes. 

With the elements of the Connaissance des Temps (7), an adjustment of 
0°37 years in the libration period, and of 0-14 years in the epoch of zero libration 
would suffice to bring agreement with (8°04) with m,=29. A comparison of the 
various sets of elements suggests that, owing to the long period of the libration, 
such adjustments are tolerable, but at the same time, the observational errors 
of both period and epoch being of the order of the synodic period of the nodes, 
little significance can be attached to such an agreement. 

It should be noted that, if the epochs of zero libration correspond to oppositions 
of the two nodes, then the epochs of maximum and minimum libration correspond 
to conjunctions of the nodes. At such epochs, the mutual inclination of the orbits 
is at its smallest value of (7, —7,), and whatever the longitude of conjunction of the 
two bodies, a reasonable satisfaction of the mirror condition would occur. 

A further approximation to the actual case of Mimas and ‘Tethys is obtained 
by considering the eccentricity of Mimas’ orbit. In this case, it is necessary 
that the epochs of zero libration correspond to conjunctions of the bodies when 
Mimas is at either its perisaturnium or its aposaturnium. With the mean motions 
used in Paper I, Table II, the rate of retrogression of the conjunction point of 
the two bodies is 218°-77 per year. ‘Therefore, if wy is the rate of advance of 
the apse line of the orbit of Mimas, it is necessary that 


(w+ = 180 my, (8.05) 


where my is an integer. At the same time it is necessary that the position of 


x 
t * 4 
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perisaturmium be such that wy, =4 4, +90 at the epochs of zero libration, where 
wy 18 the longitude of the perisaturnium of Mimas. 

While it is possible to derive a set of values for #,,, 6,, P and wy which satisfy 
(8.02), (8.04) and (8.05), and which are probably consistent with the observationa! 
accuracy of the determined values of these elements, it does not seem possible 
for wy to have the required value, unless there is an error ~30 in the positior 
of perisaturnium at the epoch of the published elements. 

It may be concluded that the system of Mimas and ‘Tethys agrees with the 
predictions of the arguments of Part II of the present paper to the first, and 
major, approximation, but it is uncertain whether the agreement is exact. 

g. Enforced eccentricities.—As is well known, the occurrence of near- 


commensurability of mean motions of two bodies produces large inequalities, 
which may take the form of enforced eccentricities. ‘The usual viewpoint ts 
that the nature of the enforced eccentricities, in particular the rates of revolution 
of the apse lines, are consequences of the magnitudes of the mean motions. ‘The 
viewpoint of the present paper is in no way in conflict with this traditional view, 
but it suggests that the mean motions have their observed values because, in such 


orbits, the enforced eccentricities will be such as to give the frequent satisfactions 
of the mirror condition. As is seen from the theorem of Newcomb (15, 16): 

‘In the case of two orbits, originally circular, if the movement of the point 
of conjunction, supposed unique, is small compared with the movements of the 
two satellites, the perturbations caused by one of the bodies engenders in the 
other an eccentricity, and the apse line of each orbit passes always through the 


point of conjunction.” 


The theorem thus states that, for near-commensurable mean motions, the 


enforced eccentricities will provide frequent satisfactions of the mirror condition. 


10. Jupiter's system.—\n this system, there are two near-commensurabilities 
where the ratio of integers is of the form(A, +1) A,, namely JIL Jland JUL JE 
From Table II] of Paper I, we have 


Ny Nn, = 0-001 82 


and Ny Ny = | 0°003 65. 


In addition, there are the Laplace relations 


nN, — + 2N, =O, (10.01) 


1, + 2l, = 180 (10.02) 


where / denotes the mean longitude (12). ‘The orbits have small inclinations 


and small but variable eccentricities. 
Ihe relations (10.01) and (10.02) mean that when JII and JIII are in 
conjunction, J | is in opposition to them, giving a close approximation to the 


mirror condition which is repeated every 7-050g days, which is less than a single 
period of J 

11. Uranus’ system.—The satellite system of Uranus is unique among the 
satellite systems of the planets in that the orbits of the satellites are sensibly 
circular and coplanar, as is shown by the extensive investigation of Newcomb 
(13), who, in giving the elements, remarked that there is but slight evidence 


* According to Griffin (17), whenever J IT is in conjunction with J I or J LLL, the inner of the 
pair is Near perijove and the outer of the pair near apojove, which would give an even clove 
approximation to mirror conditions 


| | 
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of any real eccentricity of the orbits, and no evidence of any mutual inclinations. 
Thus, within the framework of the arguments of this paper, for pairs of satellites 
the system of Uranus presents a completely degenerate case, in that any conjunc- 
tion of a pair of bodies therein will satisfy the mirror condition to a high degree 
of approximation. We would thus expect to find no preference for near- 
commensurability of mean motion among pairs of satellites of Uranus. Reference 
to Table IV of Paper I shows that this expectation is, indeed, confirmed. This 
may be considered evidence in support of the given interpretation of the preference 
for near-commensurability of mean motion found in the other satellite systems 
and the planetary system. 


12. Griffin's analysis.-After the above investigation had been completed, 
the attention of the authors was drawn to an investigation of periodic orbits by 
Griffin (17). For the restricted case of n finite satellites moving in almost-circular, 
coplanar orbits, Griffin deduced his “ symmetry theorem "’: 

“If a symmetrical conjunction occurs at any instant ¢ = f,, then the orbit 
of each satellite before and after the conjunction issymmetrical, both with 


regard to geometric equality of figures and with regard to intervals of 
time 


where a ‘symmetrical conjunction” is the coplanar case of a “ mirror 
configuration”, Griffin's “symmetry theorem” is thus a special case of the 
general mirror theorem” deduced above. 

Gjriffin further deduced that commensurabilities of the type 


A, (12.01) 


are excluded tor periodic orbits. ‘This is in conflict, not only with the detailed 
conclusions of the present paper, but also with the actual occurrence of com- 
mensurabilities in the planetary and satellite systems, as is seen at once on 
reference to Tables I-IV of Paper I. 

The contradiction is resolved by an examination of the basis of Griffin's 
analysis. He begins with a system of infinitesimal satellites moving in circular 
orbits, the periodicity of configurations being assured by the introduction of 
commensurabilities into these original orbits. He then considers the orbits that 
are obtained when the masses of the satellites are increased from zero, and obtains 
a series of periodic, nearly-circular orbits which “ grow out of’’* the original 
circular orbits. He shows that in the special case (12.01), periodic orbits do not 
grow out of the original circular ones. 

His conclusion is thus that there 1s no continuity, with varying satellite masses, 
between the finite and the infinitesimal cases when (12.01) holds. His analysis 
does not preclude the situation whereby the perturbations of finite satellites 
obeying (12.01) build up into a pair of periodic elliptical orbits, beyond the limits 
of his approximation. ‘Thetheorem of Newcomb cited in Section g isan example 
of this process. Griffin specifically points out that any commensurabilities of 
the form (12.01) can only yield periodic orbits if the orbits are elliptical, and the 
orbits would suffer large perturbations. 

As evidente of the non-existence of periodic orbits obeying (12.01), Griffin 
cites the asteroid gaps. ‘This interpretation of the asteroid gaps has been 


criticized, after Brown and Shook’s analysis, in Section 7 of Paper I. 


* "This term is Griffin's. 
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13. Conclusions.—The mirror theorem in its two aspects, together with the 
periodicity theorem, provides reasons why near-commensurabilities of the form 
(A, +1) A, should be preferred among the mean motions of the planets and 
satellites, except in the case of coplanar, circular orbits, when they should be 
absent, in agreement with observation, It appears that such orbits satisfy 
Poincaré’s conditions of stability if only that for a considerable part of the system's 
lifetime the bodies will remain near these orbits and suffer only librations. 

The frequent occurrence of mirror configurations is most likely to occur in a 
system of nearly-circular, nearly-coplanar orbits, when it is satisfied at every 
conjunction and opposition (for a pair of bodies only). An extension of the 
arguments of this paper suggests, therefore, that the most stable configuration of 
orbits for a gravitating system of point-masses (i.c. neglecting tidal friction) ts 
a system of nearly-circular, nearly-coplanar orbits, and that any gravitating 
system of bodies will spend most of its lifetime in such a configuration, whatever 
the original distribution of orbital inclinations and eccentricities may have been. 
‘The preference in the solar system for orbits of small inclinations and small 
eccentricities may be further evidence for the arguments of this paper 

‘The authors realize that these arguments are by no means rigorous (excepting 
the proofs of the mirror theorem and the periodicity theorem given in Sections 3, 
4 and 5, which are rigorous) but they hope that the support given to these argu- 
ments by the observational evidence will encourage more interest in this problem 
of celestial mechanics. It is hoped that future work will show whether or not the 
other cases of near-commensurability of mean motions in the solar system fall 
within the pattern of the arguments of the present paper. 
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Summary 


The stability of the Einstein universe is investigated using the new 
boundary conditions at surfaces of discontinuity, discovered by O’Brien and 
Synge. On certain assumptions it is shown that the Einstein model may 
represent a universe with inhomogeneities, but that, in general, these must be 
static otherwise the model is unstable 


In the case of an instantaneous disturbance in the Einstein universe, the 
processes which cause instability must involve discontinuities in the pressure 
or its time derivatives. If the matter present satisfies an equation of state, no 
physically plausible disturbance of the static condition is possible, and the 
model is stable. If there is no equation of state, small disturbances involving 
changes in the pressure or its derivatives can occur, and will, in general, lead to 
instability, though in certain special cases condensations may begin to form 
without initial instability 


1. Introduction,—Between the years 1930 and 1935 there appeared several 
papers on the stability of the Einstein model of the universe. ‘The interest in 
this problem arose mainly because it was thought that the present expansion of the 
universe might have started from a perturbation of the Einstein equilibrium 
state. Although the writers all agreed that the Einstein model is unstable in 
some way or other, there was no agreement on the nature of the process in the 
model which would cause instability, or on the conditions in which an expansion, 
rather than a contraction, would begin to take place. ‘These questions were not 
settled, and for various reasons cosmologists turned their attention to the point- 
source models of general relativity. However, these models are not without 
their difficulties—such, for example, as that of the formation of condensations 
in them—which make it premature to rule out the disturbed Einstein model. 
It therefore seems worth while to reconsider the stability of the Einstein universe 
in the light of certain recent developments in the theory of relativity. 

A study of the previous work on the stability of the Einstein model clearly 
shows that much of the discordance in the results is due to the use of different 


boundary conditions at surfaces of discontinuity. The usual procedure was to 
suppose some disturbance to occur in the model at a definite time, say f= 0, and 
to try to find, under certain assumed conditions of discontinuity of the metric 


and of the pressure and density, what happens for fo. As twenty years ago 
the question of boundary conditions to be used for such problems in general 
relativity had received no systematic attention, it is not surprising that different 
sets of conditions were used, and different conclusions reached. 

Quite recently this subject of boundary conditions at surfaces of discontinuity 
in general relativity has been studied by O’Brien and Synge (1), and by 
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Lichnerowicz (2). By the use of boundary conditions obtained by these authors 
it is possible to give a rather more definite answer than formerly to the question 
of the stability of the Einstein universe. Of course, in a subject like general 
relativity, where the fundamental equations are so far removed from physical 
intuition, sets of boundary conditions must depend on certain mathematical 
assumptions, so that in the absence of direct experimental confirmation they 
are arbitrary to a certain extent, and must be used tentatively. However, since 
essentially the same sets of conditions have been reached by O’Brien and Synge 
and by Lichnerowicz, starting from plausible but different assumptions, and 
since, as | have pointed out elsewhere (3, 4), the results of applying them to certain 
problems seem highly realistic, it is reasonable to use them in studying the 
stability of the Einstein universe. 

In this paper | shall use these conditions in two problems connected with the 
stability of the Einstein model. ‘The first of these, treated in Section 3, is the one 
stated above, of finding what instantaneous disturbances are possible, and what 
effects they would have if they occurred. ‘The second problem, dealt with in 
Section 4, is that of the existence of models similar to Einstein's in being on 
a large scale homogeneous and static, but containing local inhomogeneities and 
motions. Before discussing these problems | give a short historical survey of 
the subject in Section 2. The conclusions which | reach, given in Section 5, 
are different from widely expressed views on the subject. ‘They are, roughly, 
that the Einstein model is not necessarily unstable, nor essentially homogeneous, 
but that instability would result from certain disturbances in it which involve 
changes in the pressure. 

2. Historical survey.--Using co-moving coordinates, the line-element of the 
homogeneous cosmological models of general relativity may be written 


1 
ds* dr* + 4 sin® Ad¢*) + dt?, 


where R& is a constant which may be real, imaginary or infinite, and g(f) is an 
undetermined function of the time in terms of which the proper-density p and the 
proper-pressure p of the matter in the model are given by 

"4 


Here’ denotes differentiation with respect to ¢, and A is the cosmological constant, 
‘The special case, corresponding to a static model, in which g(t) 0 and RK is real, 
is called the Einstein universe, and for this the line-clement may be put in either 
of the forms: 


2\~1 
ds* = - (: dr*® — r*(d@* + sin* + dt*, 


ds* (: (dr* + + r* sin® Odd*) + dt*. 
4k 
The statement that the Einstein universe is unstable was first made in 1930 
by Eddington (5), whose argument was as follows. If # is real then Re! *’( =a) 
is the radius of the model (2.1), and (2.2) and (2.3) give 


3d*a 


>). 2.6 
A —4n(p + 3p) (2.6) 


| | 
(2.3) 
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For an Einstein universe with p=o we have 4mp=A; but if in such a model 
there is a slight disturbance causing a decrease in density, it is clear from (2.6) 
that d*a/dt* is positive, so that the radius starts to increase. Moreover, since this 
expansion will cause a further decrease in density, it will continue, so that the 
model is unstable. A similar argument may be used, with the same result, if 
the pressure is not zero; and evidently increases in p or p will, on this reasoning, 
cause instability due to contraction. Eddington did not give any clear indication 
how the initial change in the density would occur, but considered that the formation 
of condensations might be a cause of the instability. 

Shortly afterwards the instability of the Einstein universe was studied in 
some detail by Lemaitre (6, 7), who tried to find, in particular, the nature of the 
process which would cause instability. He concluded that the process required 
to start an expansion from the Einstein equilibrium state was a reduction of the 
pressure in the model, which he called a “ stagnation’’. The formation of 
condensations, except in so far as it is accompanied by a change in pressure, has, 
according to Lemaitre, no effect on stability. 

At about the same time McCrea and McVittie (8, g, 10) attacked the problem 
by considering the differences between the volume of the completely homogeneous. 
Einstein universe and of various static models similar to the Einstein universe 
but containing condensations represented as point singularities. “They eventually 
concluded that, to the first order of approximation, the presence of condensations 
does not affect the total volume of the model, and that their method was not suitable 
for judging the effect of condensations on stability. 

Subsequent studies of the instability were made by Dingle (11), ‘Tolman (12) 
and Sen (13, 14), with conflicting results. The method of these authors was 
of the type referred to in Section 1, and it is at this stage that it becomes clear that 
the results obtained must depend to a large extent on the boundary conditions 
used. It should, however, be remarked that the importance of such conditions 
had evidently been realized in the earlier work of Lemaitre. 

3. The conditions for instability.-In this section we shall consider the stability 
of the Einstein universe to a disturbance which takes place at a definite instant 
to, and which preserves isotropy of pressure and also spherical symmetry 
about the origin of coordinates. By a disturbance is meant an event which 
corresponds to discontinuities in the functions which characterize the model ; 
that is to say, in p, p and the metric tensor g,, (but not in A). | suppose that 
the discontinuities must conform with the boundary conditions of O’Brien and 
Synge (1), which are that at the boundary f(~x,)=0, the following quantities 
are continuous: 

where i, k= 1, 2,3, 4.and v= 1, 2, 3, and is the usual energy tensor of general 
relativity. 

The most general expression for a non-static universe, using co-moving 
coordinates and supposing spherical symmetry about the origin, but not isotropy 
for all observers, is 


ds* e‘dr* — + sin® 0dd*) + e” dt*, (3-2) 
where A, w and v are functions of r and ¢. In accordance with the assumptions 


of the previous paragraph, we suppose that the line-element has form (2.4) for 
t<o, and form (3.2) for t>o, and that the conditions (3.1) apply att=o. The 


3 
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problem is then to find the restrictions imposed on A, w and v by the boundary 
conditions and by the field equations 


Since co-moving coordinates are being used, and the pressure is supposed 
isotropic, the only non-zero components of 7',! are 

=T{=T; p, TA=p, 
where p and p are the proper-pressure and proper-density. 

In a previous paper (3) (hereafter referred to as I) I considered this problem 
for similar disturbances in a model which for ¢-—o is the non-static homogeneous 
universe (2.1), and I proved two theorems which, adapted to the special case of 
an Einstein universe, are as follows. 

Theorem 1. If for t>-o the matter present satisfies an equation of state 
p=S(p), where f is a single-valued function, then (3.2) reduces to (2.1). 

Theorem 2. If the pressure is continuous at f= o and if all its time derivatives 
at f=0+ are zero, then for ¢ 0 (3.2) reduces to (2.4). 

I shall refer to (2.1) as a Friedmann universe. Certain differentiability 
conditions are required in the proof of these theorems. Sufficient conditions 
are that p, p and g,, possess time derivatives of all orders at f= 0 + , and that these 
derivatives be at least twice-differentiable functions of r. We shall assume 
throughout this section that these conditions are satisfied. 

From Theorem 2 it follows that disturbances of the type envisaged here can 
arise only if there is some change in the Einstein model which affects the pressure or 
its time derivatives: otherwise, the model is completely stable. ‘Theorem 1 shows 
that if such a change takes place, and if there is an equation of state for t-o, the 
model will start expanding or contracting, and will, in fact, become a Friedmann 
universe. 

If there are discontinuities in the pressure and its derivatives, the stability 
may be studied, as was done in | for Friedmann models, by finding the initial 
rate of change of an element of proper volume dv of (3.2), given by 
dv sin @ From it follows that 


 A(log dz 
onl =O, 
as 


log 


where s is the proper time of an observer at rest in the coordinate system of (3.2), 
and + and — refertot=o+ andf=o-—. From (3.4) it follows that an instan- 
taneous decrease of pressure at t=o causes an expansion. ‘This was the result 
found by Lemaitre, who called the decrease of pressure a “ stagnation’’. Wa 
see also from (3.5) that the expansion will be reinforced if 0p/dt is initially negative. 
If we suppose that the matter in the model is such that the pressure will be further 
diminished by an expansion, (3.5) indicates the instability of the model to 
reductions in the pressure. Evidently, a similar argument indicates that an 
increase of pressure will cause the Einstein model to contract. 
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In reaching his result Lemaitre used arguments of questionable physical 
significance, and his methods were strongly criticized by McCrea and Me Vittie 
(10). Itis therefore interesting that it is possible to arrive at his result by a method 
which uses as its only essential assumption (apart from the spherical symmetry 
and isotropy of pressure) the boundary conditions (3.1). 

It was shown in I that p, is not a function of r, but p. may be. ‘This implies 
that if there is a discontinuity in the pressure, it must be universal, but that the 
pressure could begin to change in some isolated region. If there is an equation 
f state for ¢_-0, the universal discontinuity in pressure is necessary for instability 
to develop, because of the following argument. If there is instability, the result, 
by ‘Theorem 1, is the Friedmann model (2.1); now the boundary conditions 
require 


so that from (2.3) 


Sn(p.—p_) 


If the pressure is continuous, p, =p sothatg,=o. But owing to the equation 
f state, we then have 


Op df Op df 


Further, differentiating (2.3) we obtain 


Equations (3.6) and (3.7) give 
=0. 


In this way it may be shown that (d"g/dt"), =o for all n, so that (2.1) reduces to 
the Einstein model and there is no instability. 

Unless, therefore, we are prepared to suppose that there is a universal discontinuity 
m the pressure, we have to rule out instability when there is an equation of state for 
t -o. ‘The most realistic possibility for the development of instability seems to 
be the case where there is no such equation, and where the pressure is continuous, 
but dp dt (or a higher derivative) is discontinuous in some small region of space. 
In this case expansion or contraction starts according to (3.5). ‘This case will be 
studied in the next section. 

4. Approximate Einstein models,—I\n the previous section | have given a solu- 
tion of the problem of the stability to certain types of disturbance of the Einstein 
model considered as completely homogeneous and static. One would naturally 
wish to interpret the homogeneity as the result of averaging the density and 
pressure over regions of a certain size, because inhomogeneities might be expected 
ina physical realization of the model. But for this interpretation to be permissible 
one must be sure that inhomogeneities could in fact exist in the Einstein universe. 
Stated rather more precisely, the question is whether there are solutions of the 
field equations representing static, inhomogeneous models which are similar 
to the Einstein model on a sufficiently large scale. If such solutions do not 
exist, then we are not entitled to regard the Einstein universe as a convenient 
model for something more complicated than it actually represents: in this case 
we must appreciate that any physical realization of the Einstein model must be 
strictly homogeneous, at least on a macroscopic scale. 
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Such static, inhomogeneous solutions do indeed exist; some have been 
published already, and I shall prove the existence of others later. Weare therefore 
allowed to think of condensations in the Einstein universe without being in 
danger of any logical contradiction. 

Ihave dealt with this point at some length because some writers on the 
Einstein universe seem to think that it does not allow the presence of condensations 
at all. ‘There has been a tendency to confuse the existence of condensations 
with the formation of condensations. ‘The latter is a different question, involving 
non-static solutions of the field equations. 

The formation of condensations may be regarded from two points of view, 
First, we may consider processes in the strict Einstein model involving changes 
of the pressure and its derivatives, discussed in Section 3, which lead to changes in 
density and metric. Secondly, we may suppose that we start with a non-static 
model which resembles the Einstein model in the large, but which includes 
regions where density changes are taking place. In both cases we want to know 
whether the process of condensation can continue in bounded regions without 
causing the model to expand or contract as a whole. But in the second case we 
do not have to satisfy initial conditions at (= 0, because the model is not supposed 
to start from the strict Einstein state. 

To study these problems | consider what | call ‘approximate Einstein 
models”. In these the variations from the Einstein metric are either confined 
to small regions of space, or (if not) are small except over limited regions. 

Several approximate Einstein models were obtained by McCrea and Mc Vittie 
(8,9, 10), who found line-clements corresponding to Einstein universes contain- 
ing particles represented as point-singularities. As an example, a special case 
if a solution of McVittie (15) may be quoted: 


1 + 2 


(m, 2r)[1 + 
1 +(m 2r)[1 +(r* {1 + (r*/4R*)?* 
(dr* + r*d0*-4-r* sin*® Odd*). (4.1) 


ds? 


dt® 


lf m=o this reduces to the Einstein model in the form (2.5), and if R= & it 
gives the Schwarzschild metric in isotropic form. ‘The pressure and density, 
obtainable from expressions given by McVittie, are functions of A, R, m and r, 
ind we may not put po without abolishing m. 

Models such as (4.1) suggest that inhomogeneities could exist in an Einstein 
universe. ‘This is the conclusion also from the work of Raychaudhuri (16), 
who found that there are exact static solutions of the field equations corresponding 
to an Einstein model containing spherical pockets of radiation of variable density 
Indeed, as we shall see later, models can be constructed in which static, spherically 
symmetric regions of arbitrary density are embedded in an Einstein model. 

In the approximate Einstein models which we shall consider, we suppose 
that the inhomogeneity, in which there may be motion, has spherical symmetry 
about the origin of coordinates, and lies within the sphere r= a (constant). We 
assume that co-moving coordinates may be used so that the line-clement has the 
form (3.2) inside r=a, and also that the pressure is isotropic. Outside r=a 
the line-element is supposed to be that of the Einstein model (2.4), and across the 
boundary r = a, we suppose that the conditions (3.1) are satisfied with the obvious 
change that since the boundary is now x, = constant, the quantities which must 
be continuous are 
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where i, k= 1, 2, 3, 4 and yw, v= 2, 3, 4. Evidently the model may contain any 
number of non-overlapping spheres of this type, each corresponding to an inhomo- 
geneity. ‘This model is somewhat similar to one used by Lemaitre to study 
condensations in the expanding universe. 

The conditions (4.2) give the following restrictions on A(r, t), w(r, t) and 


u(r, t) in (3.2): 


2 


A(a, t) = log t)=2loga, t)=o, 


(4-3) 
(a t)=o. 


The continuity of p, required by (4.2), is secured by (4.3). 
Let us consider first the case of a purely static inhomogeneity. ‘Then we may 
take (3.2) as 
ds* = — e*dr* — r*(d0* + sin*® 6dd*) + 
where A and v are functions of r only. The field equations (3.3) give (Tolman (17)) 
I 
8r7',' =- Srp = —€ ‘( + x) + (4-4) 


where ‘ denotes d/dr. Equation (4.6) may be integrated to give 
Sar 
where 


| prtdr+ K, 


the density p being an arbitrary function of r, and the constant of integration K 
being chosen so that g satisfies the following, which is required by (4.3) : 


(1 I 
(7-34): 


Equating the right-hand sides of (4.4) and (4.5), and substituting for A from (4.7), 
we find a second-order differential equation for v. It is clear that this admits 
a solution v(r) containing two arbitrary constants which may be chosen to satisfy 
(a) = v'(a)=0 

without further restriction on g(r). Hence for an inhomogeneity with an arbitrary 
density distribution there exists a solution satisfying the boundary conditions at 
y=a. This leaves no doubt that static inhomogeneities can exist in the Einstein 
universe, though they may not all be non-singular. 

One further point about the static solutions is worth noticing. If p is constant 
for r<a, it follows from the identity 

dp|dr = —\v'(p +p) 

that v=o and the solution reduces to the homogeneous Einstein metric. In 
particular, in a model with p=o, static inhomogeneities cannot exist. 

‘Turning to non-static inhomogeneities, we face a much more difficult problem 
because it is necessary to integrate the field equations obtained from (3.2). 
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I shall not attempt the complete solution of this, but confine attention to the 
approximate solution which arises by supposing that the metric (3.2) inside the 
inhomogeneity differs only slightly from the Einstein metric (2.4). 

First, let us suppose that v, 0A/0t, dw/dt and all their derivatives are small 

quantities of the first order. Then, neglecting small quantities of order higher 
than the first, the field equations for (3.2) (which are given in I) reduce to the 
following : 
~8nT A= —82p = + — + A, (4.10) 
=0= —w' —4w'(w—A), (4.11) 
where * denotes 0/dt, and ’ denotes 0/dr. 

Integrating (4.11), we have 

= e’w'*h(r), (4.12) 


where A(r) is an arbitrary function, which for convenience we give the value 


2 


where k(r) is another arbitrary function. 
We now suppose further that 
=r*{1 +a(r, (4-14) 
where « and its derivatives are small quantities of the first order. ‘Then (4.12), 
(4.13) and (4.14) give, to first order, 


A 
€ = (1~ ja) (1+a+ra’ +k). (4.15) 
With (4.14) and (4.15), equations (4.8), (4.9) and (4.10) give after some calculation 


— 8nT;,' = =A ja 


Equating the right-hand sides of (4.16) and (4.17), we have 


a’ ak 


~8nTA= —82p=A- 


which gives, on integration, 


= (: pa) je + ~ = p(t). (4.19) 
‘The boundary conditions (4.3) require that a, v and v’ vanish at r=a for all t, 
and as the continuity of A requires k(a)=0, we deduce that p(t) in (4.19) must 
be identically zero. If we put 


k 2 2 
Ar, t)=a— (: (4.20) 


22* 
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(4.19) may then be written 
B- ja) (4.21) 
and we find from (4.16) and (4.18) 
3p 
p Po) = Re (4.23) 


where Py, po are the values of p and p in the homogeneous part of the model. 

Equations (4.22) and (4.23) give the pressure and density inside the inhomo- 
geneity in terms of two functions v and 8, which have to satisfy (4.21) and also the 
following conditions at r=a: 


Wa, th=v'(a, th=o (4.24) 
B'(a, t ( k’( 
t)=o, (a, t)= ay a) = constant. (4.25) 


In deriving (4.25) we remember that k(a)=0, but k’(a) is an arbitrary constant. 
The conditions (4.24) and (4.25) require 


pla, t)= Po, p'(a, th=o, (4.26) 
pla, t)— py = constant. 


‘To investigate the solution further, additional assumptions are necessary. 
Let us suppose that inside the inhomogeneity matter obeys a single-valued 
equation of state p=/(p) (excluding, for the moment, p= const. and p=const.), 
that p and p have derivatives with respect to r of all orders for r<a and 
left derivatives with respect to r of all orders at r=a, and that p, p and these 
derivatives are at least twice-differentiable functions of ¢. Using (4.22) and 
(4.23), we find from (4.21) 


and from the equation of state we may write this as a partial differential equation 
for p in the form 


where F is a known function. If in this we put r= a and use (4.26) we find that 
p” is constant atr=a; and differentiating (4.27) with respect to r, it then follows 
easily that p’” also is constant atr=a. In this way it may be shown that all the 
derivatives of p with respect to r are constant at r = a, so that p is not a function of 
the time immediately within r=a. Using the differentiability condition on p, 
we can use (4.27) to extend this result to any point inside r=a. ‘Thus p does 
not depend on ¢ for r<a, and it is easy to see that v, p and 8 do not depend on ¢ 
either. It can also be shown that if p= const. the solution is not a function of f. 
Thus given an equation of state (other than p = constant) and appropriate differentia- 
bility conditions of p and p, the inhomogeneity must be static. ‘This result would 
be expected from Section 3, in which it was shown that in an Einstein model 
with an equation of state any disturbance must be universal. 

If there is no equation of state, we may suppose that we are given v(r, ¢) 
(i.e. the pressure) for r<a. In this case the right-hand side of (4.21) is known, 
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and the equation may be integrated to give 

B= A(rje® + B(rje + H(r, t). (4.28) 
where A(r) and B(r) are arbitrary subject to conditions at r=a, and H(r, t) is 
a particular integral, known if vis given. (‘This form of 8 applies also in the special 
case where p=const.=p,, which gives H(r,t)=0.) ‘The functions A and B 
would be determined if we were given appropriate initial conditions (e.g. p{r, 0) 
and p(r, 0)). 

From (4.28) it is clear that 8 will increase or decrease exponentially with time 
except in special cases, such as, for example, that in which A and H vanish. This 
means that the metric (3.2) for <a will depart more and more from the Einstein 
state, and the assumptions of our method of approximation will eventually cease 
to apply. Before, however, we conclude that this means instability, we must 
examine a little more closely what is implied if 8 is an increasing or decreasing 
function of the time. 

From the assumption that for r>a the metric is that of the Einstein universe, 
it follows that the inhomogeneity is necessarily confined within rea, The 
effect on stability of the motion inside r =a is to be judged by the changes in the 
proper radius and proper volume of the inhomogeneity, which will affect the 
proper radius and proper volume of the whole model. If the proper radius and 
volume of the inhomogeneity tend to infinity with time, then the whole model is 
unstable because of expansion ; if these quantities tend to zero, the whole universe 
contracts, though in this case we cannot be certain that this contraction will 
continue after the volume of the inhomogeneity becomes zero. If, however, 
the change in proper radius and proper volume of the inhomogeneity is zero, or 
remains small, then there is no reason why the motion inside r= a should cause 
instability. 

The proper radius / and the proper volume v of the inhomogeneity are, to 
the first order of approximation, 


“a a 


12 
| (: 7a) (1 +48 + 4rB’)dr, (4.29) 


2\-12 
sind dd db dr=4n\ r* (: (1+ + dr. 
(4-39) 


In these expressions A(r) does not appear because the terms which involve it may 
be integrated, and then vanish on the assumption that k(a)=o (already made), 
and that k(r) tends to zero with r more rapidly thanr. It may be stated here that 
the physical significance of k(r) is connected with the spatial geometry of the 
model for r<a. When the other unknown functions, corresponding to the 
density and pressure, are determined, this geometry is not completely fixed, 
being unspecified to the extent of the arbitrary function k(r). ‘The arbitrariness 
disappears if appropriate initial conditions are given, as we shall see later in an 
example. 

The contributions to the integrals (4.29) and (4.30) which could make / and v 
change rapidly with time are those involving, through f and fp’, A(r)e and 
H(r, t). In general, whatever the form of H, the dependence on Ae will make 
/ and w change rapidly and cause instability. However, in certain special cases 
A and H may be such that / and v remain nearly constant in time, even though 
at certain points in ra, 8 may become large. For example, this will be so if 


a 
0 
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H ~o (which implies p= p,) and A satisfies 


with Asso. In other words, the inhomogeneity may remain roughly the same 
size, though the density (which depends on 8) changes considerably at some 
points in it. Thus condensation, or rarefaction, may begin without causing the 
whole model to contract or expand, though this will be so only for these special 
cases, and even in these we cannot be sure (since our method of approximation 
breaks down when £ becomes large) that the integrals will not eventually begin 
to change rapidly, thus making the model unstable. 

The solution, which is given essentially by (4.20) and (4.28), depends on a 
function -(r, t), whose prescription amounts to the specification of the pressure 
p(r, t); on the functions A(r), B(r) which may be considered as determined by 
the specification of p and dp/dt as functions of r att =o; and also on A(r), determined 
by the metric att=o. We have not so far supposed the model to have originated 
in any particular way: we may indeed regard it simply as a refinement of the 
Einstein model which represents, in so far as it is stable, a universe static and 
homogeneous on a large scale, but containing local irregularities and motions. 
On the other hand, we may ask whether a model of this type could have originated 
from the strict Einstein universe by a perturbation of the kind described at the 
end of Section 3. The result of this enquiry, to which we now turn, may be 
expected to restrict the arbitrariness of the functions A, B and k. 

In addition to the conditions (4.3) at r= a, we must now impose the following 
at t=o (which result from (3.1)): 


Ar, ©) = (1 - Mr, 0) =0, 
0) = 2logr, Ar, 0) =0, 
ur, 0) =o. 
Equations (4.3), (4.14), (4.15), (4.20) and (4.31) together give 
alr, 0)=a(r, 0) =0, 
= 0, 
a(a,t)=«a'(a, t)=o, 
(a, t)=v'(a, t)=0, 
alr, t). 
Using (4.28) for 8, and putting t= 0, we see that from (4.32) 
A(r, 0) = A(r) + B(r) + H(r, 0) = 0, 
A(r, 0) = R-{A(r) — B(r)] + H(r, 0) =0, 
which determine A and B since H(r, t) may be supposed known as (7, ¢) 1s given. 
In addition, (4.32) requires that 


A(a)e!® + Blaje + H(a, t)=0, 
A'(a)je!® + B(aje"*® + H'(a, t) =o. 
It is easy to see that a function »(7, t) may be chosen so that the problem has 
a solution. For example, if we take 


v= “t, 


(4.33) 


(4.34) 


+ 
3 
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where « is a positive constant and g(r) is a polynomial, we get a non-singular 
solution satisfying all the conditions (4.33) and (4.34). ‘The solution corresponds 
to a disturbance in the Einstein model which takes place inside r = a, and in which 
the pressure begins to change from the equilibrium value at =o. Although the 
pressure difference dies away eventually and may be arbitrarily small at all times, 
it causes instability in the general case because (4.33) requires A(r) a2 © so that the 
integrals (4.29) and (4.30) will eventually change rapidly, causing the whole model 
to contract or expand. However, in special cases this model can correspond to 
condensation or rarefaction inside ra without initial instability because the 
constants in g(r) can evidently be chosen so that / and © do not change rapidly 
with time, even though A(r) <0. 

I shall now summarize the results of this section. We have seen that exact 
static solutions of the field equations exist corresponding to an Einstein universe 
with inhomogeneities. If there is an equation of state (other than p= const.) 
no non-static approximate models with inhomogeneities are possible, but if there 
is no such equation, non-static models exist and are in general unstable, though 
there are special cases in which motion inside r = a can occur without instability, 
including some in which condensations may begin to form. If a non-static 
approximate model results at ¢ =o from some change in the time derivatives of p 
inside a regionr = aof the strict Einstein model, its future is completely determined 
given p(r, t) for to, r <a, and in general it will be unstable no matter how small 
is the pressure difference in the homogeneous and inhomogeneous parts of the 
model; but in special cases, condensations can begin to form without instability, 

5. Conclusion.-From Section 4 it seems that we are justified in regarding 
the homogeneity of the strict Einstein model (2.4) as a smoothed out representation 
of a static universe with inhomogeneities (except in the case of a completely 
pressure-free model). Its static character must, however, be interpreted more 
strictly because it appears that, except in particular cases, systematic motions 
are accompanied by instability. Speaking roughly, we may say that the Einstein 
universe is not essentially homogeneous, but it is static*, except for the special 
motions mentioned in Section 4, and also, possibly, for rotational motions, which 
have been excluded from this work by the assumption of spherical symmetry, 
and which may not be accompanied by instability. 

If we accept this view of the Einstein universe, it is clear that instability will 
result from any process which introduces motion, subject to the exceptions referred 
to above. [Except in the case p ~o, such a process will not follow from the mere 
presence of condensations of matter, because static solutions corresponding to 
them exist. From Section 3 it appears that the processes concerned are those 
which involve discontinuities in the pressure or its derivatives, since in the absence 
of these no instantaneous disturbance can occur. However, if there is an equation 
of state after the disturbance, the only permitted discontinuity of the type we 
have considered is a universal discontinuity in the pressure itself, and this is not 
likely to occur. (If it did, a Friedmann model would result.) ‘Thus a model 
with an equation of state would probably be stable because no change whatever 
could take place in it. 


* I think it would be wrong to conclude from this that molecular motions contributing to the 
pressure are thernselves sufficient to cause instability. In all ordinary cosmological investigations, 
including this one, p and p are macroscypic quantities, and the theory is not intended to apply 
at the molecular level. 
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In a model without an equation of state, instability could arise through a 
gradual pressure change (i.e. a discontinuity in the time derivatives of p at t =o) 
inside a limited region; such a disturbance,e ven if small, will in general lead 
to condensation or rarefaction, and will cause instability through an expansion 
or contraction of the region in which it occurs. However, in special cases, the 
disturbance could die out or lead to condensation or rarefaction without immediate 
instability. 

A universal decrease in the pressure in an Einstein model (if such could occur) 
would initiate an expansion, and an increase would start a contraction. For a 
gradual pressure change in a limited region, whether there is an expansion or 
contraction depends on the sign of di dt and dv dt, which in turn depends on the 
pressure distribution, and the initial distributions of p and dp dt 

If the Einstein universe contained evolving stars, changes of the pressure 
of radiation would presumably be sufficient to cause instability. If, however, 
one pictures it as filled with uniform gas at low temperatures (Bondi, 18), or as 
a static mixture of gas and condensed, non-radiating matter, this work suggests 
that it would be stable. 


As the boundary conditions (3.1) are still to some extent tentative, it is worth 
while emphasizing how convincingly they seem to work in the problems dealt 
with here. All the functions which they leave arbitrary have a reasonable 
physical interpretation, and the data required to determine the solutions are very 


much what one would expect on classical grounds. This has certainly not been 
true of sets of boundary conditions used in some of the previous work on this 
subject, 
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Summar\ 


The intrinsic colours of type I galactic cepheids are determined from the 
photoelectric colours of Eggen, using a method in which the observed colours 
are corrected for reddening by means of a model of the reddening medium 
which gives the colour excesses of B stars as a function of distance and the 
galactic coordinates. Estimates of the distances and reddening of the cepheids 
are found from the model by converting colour excess into photographi 
absorption and solving the distance modulus equation using Baade’s zero 
point of the period-luminosity curve. Solutions by least squares are then 
made for the coefficients in linear relations for the intrinsic colours at 
maximum and minimum light for cepheids belonging to Eggen’s types A, B 
and C. The following solution for the intrinsic colours on the P, V System 
is adopted 

‘Types A, B and (P-—V) max log P 
Types A and B (P-V) min of log P 
Type C (P-V) min log P 


It is shown that these intrinsic colours give colour excesses and distance 
moduli which support the hypothesis that cepheid distances and periods 
are not correlated. ‘The intrinsic colours are much bluer than hitherto 
supposed, and at maximum light there is some similarity with the colours 
of cepheid variables in the Small Magellanic Cloud 

The intrinsic colours are applied to an examination of the connection 
between differential galactic rotation and the zero-point of the period 
luminosity curve. It is found that Baade’s zero-point gives an Oort constant 
A= 20 km/sec/kpe for a group of nearby cepheids for which the average 
absorption is about 2 magnitudes per kpe. The recovery of Joy's original 
value of A for the cepheids is shown to be due to the fact that the greater 
absorption required by the intrinsic colours compensates for the increas 
in Joy’s distance moduli made necessary by Baade’s revision of Shapley’s 
zero-point of the period—luminosity curve 


Introduction.—'Vhe determination of the interstellar absorption of light by 


means of the observed colours of stars requires a knowledge of their unreddened 


(or intrinsic) colours, and a value for the ratio of the total photographic absorption 
to the selective absorption or colour excess. In the case of the B stars, there is an 
adequate sample of nearby objects in high galactic latitudes from which the first 
requirement, the intrinsic colours as a function of spectral type, can be obtained 
by supposing that the amount of reddening which their light has undergone is 
negligibly small (1, 2), or can be estimated (3, 4, §). On the other hand, cepheid 
variables belonging to Baade’s population type | are strongly concentrated in the 
galactic plane. They are, moreover, intrinsically more luminous than the B stars 
and thus can be seen at greater distances. Consequently it is unlikely that there 
are many cepheids whose observed colours can be considered as unreddened 
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Attempts have been made to avoid the difficulty presented by the spatial 
distribution of the type I cepheids by an appeal to the colours of non-variable 
supergiant stars. ‘The method involves two steps (i) the assignment of spectral 
types to the cepheids, and (ii) the determination of the intrinsic colours of super- 
giants as a function of spectral type. Originally (i) was carried out by means of 
a relation between period and median-spectrum, but this has been modified by 
more precise spectral classifications on the Yerkes System which lead to well- 
defined period-spectrum relations at rnaximum and minimum light (6). This is 
probably the most certain part of the method. In step (ii), in addition to having 
to decide upon how much reddening should be allowed for in assessing the 
intrinsic colours of supergiant stars in a variety of galactic latitudes, it is felt 
that some consideration should be given to the basic assumption that variable and 
non-variable supergiants of the same spectral type have the same intrinsic colour. 

The method which refers to the colours of non-variable supergiants is 
discarded in this paper in favour of a more direct procedure, similar to that used 
in (4) for the determination of the intrinsic colours of the B stars, in which indivi- 
dual corrections for reddening are applied to the observed colours, the amount 
depending upon distance and galactic latitude. Earlier attempts (7, 8) to 
determine the intrinsic colours of cepheids in this way were based upon photo- 
graphic observations and upon a cepheid distance-scale which has been 
considerably modified by Baade’s revision of the zero-point of the period- 
luminosity curve. In the present work a determination of the intrinsic colours 
is made from photoelectric colours using the new cepheid distance-scale. 
Eggen’s observations (9) of cepheid colours reduced to the P, V System (10) 
are used, thereby making it possible to compare the colours of galactic cepheid 
variables with the P, V colours of cepheids in the Small Magellanic Cloud (11). 

1. The method 

Let P and V be the photographic and photovisual magnitudes on the P, V 
System, P-V being the colour, and let us suppose that the cepheid colours at 
maximum and minimum light can be written in the form 


(P-V)max = a+ blog P+E, 
(P V)min C+ d log P+ E, 
where P is the period and E, the colour excess, is some complicated function of 
distance and the galactic coordinates. ‘The first two terms in these expressions 
are the intrinsic colours at maximum and minimum light. ‘The intrinsic colour 
is assumed to be linear in log P since the colour amplitudes 
A\(P V)= (P V) min V) max, (2) 
which are probably unaffected by space reddening, were found by Eggen (9) 
to exhibit a strong linear correlation with log P. 
Since the total photographic absorption 4,, due to selective absorption 1s 


proportional to the colour excess, i.e. A xE, we have, considering only one 
of the equations in (1), 


(1) 


ve 


x{(P-V) —(a+ blog P)} = A,,. (3) 


If A,, for a given cepheid can be inferred from the colour excesses of nearby 
B stars on the system of Stebbins, Huffer and Whitford, we have A,, = x,£, and 
equation (3) becomes 


(x xi){(P-V) — (a+ blog P)} = £,. (4) 
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Now the photoelectric colours are probably sufficiently accurate to enable £, 

to be treated as the observational quantity in the equation of condition (4). In 

this case, it is legitimate to apply the ordinary method of least squares to solve 

for a, b and y/y,. But the determination of y y, is properly a photometric 

problem, and it would be more appropriate to adopt a value for the ratio y y, 

thereby throwing all the weight of the solution into the determination of a and /. 
In (8), equation (3) was used in the form 


a+blog P+cA,,=(P-V), 


and a least squares solution was made for a, 6 and c, using photographic observa- 
tions of cepheid colours. ‘This procedure is, however, only correct if A ,, is known 
accurately or with an error considerably smaller than the error in the colours, 
otherwise a multi-error least-squares solution would be more appropriate, if 
manageable. In the present work, the observed colours are corrected for 
reddening before applying the least squares method, the correction being derived 
from the colours of B stars. 


9 


If the colour excesses of the B stars can be approximated by a function involving 
the distance 7, the galactic latitude 4, and if possible the galactic longitude /, the 
formula for colour excess may be used to estimate the reddening experienced by 
a given cepheid. ‘The correction for reddening is then (y,x)E,(r, 4, 4), or 
(x,/x)E\(r, 6) if the average over the observable sector of longitudes has to be 
taken. In order to establish a suitable interpolation formula for £,, it 1s necessary 
to adopt a model which is capable of representing the averaged structure of the 
reddening medium. 


2. The reddening medium 

It is well known that the variation of stellar density with distance = from the 
galactic plane can be adequately represented by the relation D(z) = Dye *. 
For the A stars alone, 2 is approximately 100 pc, but considerably greater values 
of 6 are found for later spectral types. ‘Taking all the stars together, Oort's 
data (12) can be well represented for z < 500 parsecs by the simple exponential 
formula with B=230 pe. But from the appearance of extra-galactic nebulae 
seen edge-on, the obscuring matter might be expected to be more concentrated 
towards the galactic plane than the stars in general. ‘The exponential distribution 
for the obscuring matter with 8 = 105 pe was adopted by Williams (13). 

The averaged density distribution of the reddening material could be 
determined from the colour excesses of the B stars if a curve showing the 
dependence of FE sin 6 upon 2(=rsinb) could be established from the data. 
‘Thus if e(z) is the coefficient of selective absorption at height 2 above or below 
the galactic plane, we have 

(z)=dE sinh dz, 


and the density of the reddening material, which is proportional to ¢(z), follows 
by differentiation of the curve. With c(z)= ce *”, the colour excess for an object 
at distance r in galactic latitude 4 is given by 


sind 
E(r, b) | exp ( )ar 


( ( ysind 


This form for the reddening has been used extensively (14, 1§, 16, 17). 


(6 
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In the present application, use is made of three independent determinations 
of the parameters ¢, and 8, two of which give the colour excess distribution 
averaged over the sector of galactic longitude covered by the data on B star colours, 
the third giving the distribution within small zones bounded by circles of galactic 
latitude and longitude in each of which average values of the parameters are used. 

2.1. Reddening according to Oort.—-In Oort’s work (3), a correction of 

+ o”-04g was applied to the zero-point of the colour excesses of the B stars on the 
system of Stebbins, Huffer and Whitford (SHW). Using distances r calculated 
from secular parallaxes and the rotational effect in the radial velocities with 
A= 1g km sec kpe, Oort obtained a curve giving E, sin basa function of z. This 
curve was differentiated by Parenago (15) who obtained 


log (dE, sinh dz) = — 3°42 —0°00356 2. 


This corresponds to the exponential distribution of c(z) in equation (5) with 

(y= 038 per kpc, and B=122 pe. ‘These values of the parameters give the 

general distribution of the reddening medium averaged over the longitude sector 
: from /=o° to 200. ‘The value of c, determined in this way is rather sensitive to 
the correction applied to the zero-point of the colour excesses. 


2.2. Reddening according to Parenago.-Equation (6) written in the form 


( ( rsinb (> 
sinh )). 7) 


has been developed and applied extensively by Parenago in two papers (15, 16). 
In the first paper, the following mean values of the parameters were adopted for 
the whole sky: ad) = per kpc, B= 100 pe. The second paper takes account 
of irregularities in the absorption, both in latitude and longitude. ‘The deter- 
mination of the parameters a, and £ has been described by Parenago (16, p. 131) 
as follows : 


We determine for each star whose photoelectric colour index is known, for every place 
with counts of extra-walactix nebulae, etc., the values of Ang and 7 Assuming some 
initial value of / (that taken was [f= 100 pe), we calculate for each determination of Ags 
and r the value of ay from formula (7) and place it on a diagram as a function of / and + 
Further, having taken account of the unavoidable errors in the determination of A,, 
and r, we group our determinations of ay in such a way that without exceeding the limits 


of these errors, we determine regions on the diagram such that a, may be regarded as 
constant in each, In this way we get a whole series of such regions for each of which 


we can determine by the method of least squares the most probable a, and the correction 
to the assumed // 


Parenago found that within the limits of error of the determinations, the value 
of 6 was independent of the longitude, the mean value being 8 = 100 + 4 pe. 
Values of a, were assigned to 122 regions of which 81 lie within +10 of the 
galactic equator. for regions close to the equator it is reasonable to suppose that 
most, if not all of the weight in the determination of the values of a, comes from the 
colours of B stars. Parenago used Florya’s intrinsic colours (4, p. 12) which 
correct the zero-point of the SHW colour excesses by + 0-04, and obtained the 
absorption from A,,=10F,. The distances of the B stars were calculated for 
adopted absolute magnitudes. 


. 
7 
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2.3. Reddening according to van Rhijn.—I\n an investigation on the differential 
absorption averaged over the longitude inters al from /= 320 to 210 , van Rhiyn 
(17) obtained 
o”-031 


‘This corresponds to = 0™-25 per kpe and = 124 pe in equation (6). Van Rhyn 
used the SHW colour excesses with a zero-point correction of + 0-03, as calculated 
in(5). ‘The distances were derived from interstellar line intensities, the rotational 
effect in the radial velocities with A = 18 km, sec/kpe, secular parallaxes, and from 
the density law and luminosity function. It was shown that the difference 
between the observed mean colour excess E, and that found from the formula 


with distances estimated in the above mentioned ways was quite generally small. 


‘The formula was also shown to hold for the longitude sectors 320° to 110°, and 
110° to 210° separately, although derived from stars at all longitudes combined. 


3. Solutions for the intrinsic colours 

According to Eggen (9) cepheids may be divided into types A, B or C on the 
basis of relations between period, light-amplitude and colour-amplitude, ‘The 
light-amplitudes of types A and B are strongly correlated with their periods, 
whereas the type C amplitudes are nearly uncorrelated. ‘The colour-amplitudes, 
on the other hand, are strongly correlated with the periods for all three types, 
but the A and B types average about o™-2 greater in colour-amplitude than the 
C types. Following Eggen, we shall suppose that cepheids of types A, B and C 
at maximum light satisfy the same intrinsic colour relation. We then have to 
find two relations for the intrinsic colours at minimum light, one for types A and B, 
and another for type C. The solutions given in this Section are made from 
,7 cepheids of which 28 belong to types A and B,andg totype C. The intrinse 
colours of type C cepheids at minimum light cannot be determined from th 
limited amount of data. ‘The relation may, however, be inferred by adding the 
type C colour-amplitudes to the intrinsic colours of types A, B and C at maximum 
light. 

The equation of condition for the intrinsic colours of th cepheids 1s 


a+blog 
where (P-V),,,,, the observed colour corrected for space reddening, is given by 


(P-V) corr = — 4). 


lhe reddening has been estimated by means of the distributions of colour excess 
in Section 2, E,(r,6) being calculated for each cepheid by successive 
approximations in the equation for the distance modulus 


b)—5, 


where M is the absolute magnitude from Shapley’s period-luminosity curve (18) 
with a zero-point correction of —1°4 magnitudes. Values for x, and y,/x are 
also required. We have adopted y,=9 as confirmed by Whitford (19), and 
x,/x=1°9 from the work of Seares and Joyner (20). 

The solutions for the intrinsic colours are given in ‘Table I, ‘The root mean 
square errors of the coefficients in the linear relations are also given in the ‘Table. 
‘The root mean square error of an equation of condition is +0™14 at maximum 
and minimum light for the solutions designated Oort and van Rhijn, in which the 


\ 
| | 
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Tasie I 


Intrinsic colours of Cepheids at maximum and minimum light 


Reddening Maximum Minimum 


adopted 
pest ‘Types A, B and C Types A and B Type C 


> 


‘08 log P -68 log P o™- 10 o-67 log P 
Oort 
+ + +:°10 +-o7 


-18 log P -78 log P | o!-77 log P 
Parenago 
+ + -o6 + 'O7 + + 


0-12 log P o™-12 0-74 log P -o4 0-71 log P 
van Rhiyn 
+ ‘07 + + + + -o8 


mean colour excess distributions were used. ‘This error is reduced by 36 per cent 
to +0"-og in the Parenago solution which takes account of irregularities in the 
longitude distribution. Since the difference between the intrinsic colours in 
the solutions Parenago minus Oort, and Parenago minus van Rhijn are not statis- 
tically significant, we shall adopt the intrinsic colours according to the Parenago 
solution on account of its greater weight. ‘Thus we adopt the following relations 
for the intrinsic colours of type I galactic cepheids 


Types A, Band C: (P-V)max = +017+0™18log P, | 


Types A and B: (P V)min + + log P, (8) 
Type C: (P-V)min= —0™-09 + 0™-77 log P. | 
lhe colours which lead to this solution are shown in Fig. 1. 


4. Some properties of the adopted solution 
Eggen (9) has given a set of relations for the intrinsic colours which were 
derived from eight cepheids. ‘These are the bluest, for their periods, of the type I 
cepheids for which photoelectric observations of colour are available. Eggen’s 
intrinsic colours on the P, V System are as follows (21): 
Types A, Band C; (P-V)max o™-04 + o™-60 log P, | 
Types A and B: (P-V) min o™-13 + 1™-22 log P, (9) 
Type (P V)min o™-32 + 1™-20 log P. J 
The relations (8) and (g) are compared in Fig. 2. The most striking feature of 
the comparison is the marked difference between the slopes of the lines in the two 
sets of intrinsic colours. ‘The difference amounts to about seven or eight times 
the root mean square error in the present determination of the coefficients of log P. 
The significance of this difference is now examined. 


| | | 
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Fic, 1.—The intrinsic colours of cepheid variables at maximum and minimum light according 
to the relations (8). The observed colours corrected for space reddening are shown by open and 
filled circles. 
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Fic, 2.—Comparison between the adopted intrinsic colours (full lines) and Eggen's intrinsic 
colours (broken lines). Open circles give the observed colours at maximum light for cepheid variables 
in the Small Magellanic Cloud. 
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4-1. Correlation between distance and period.-\f the distances of the cepherds 
are calculated from the relation 

logr=1+h(m—M— xB), (10) 

where E is the difference between the observed colour and the intrinsic colour 

for a given period, and y= 4°75, it is found that Eggen’s intrinsic colours (g) 


give distances which show a strong correlation with the periods, as shown in 
Fig. 3. We wish to enquire if such a correlation is physically or statistically 


os 
LOG. PERIOD 


Fic. 3 Correlation bet on cephiid distances and period The linear regression of lo pr 
m log P is plotted in the figure @5 repe on log P The upper and lower panels give the correlatior 
for distances based upon Eqgen's intrinsic colours (9), and the adopted solution 0S) respective 


(pen circles denote group means 


possible. On the physical side, the only known correlations involving the periods 
is a correlation between the z-distribution and log P, due to Parenago and 
Grigorieva (22) who found that the concentration towards the galactic plane is 
greater for the longer periods than for the shorter periods, and the related 
conclusion of Shapley and McKibben (23, 24) that the longer-period cepheids 


favour regions of high stellar density. On the statistical side, we have to decide 


if the large correlation which Eggen’s intrinsic colours give could be ascribed to 
chance if the distances and periods of cepheids in the parent population were 
uncorrelated. 

The statistical significance of the correlation may be tested by means of 
Fisher's x—{ distribution for the correlation coefficient which is approximately 


normal with mean 


and standard deviation 
12 


2(n—1)? | 


n being the number in the sample and p the correlation coefficient in the parent 
population, Consider first the case p =o, that is when the distances and periods 
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of cepheids in general are uncorrelated. ‘Ihe correlation coefficient for the 
linear correlation of Eggen’s values of logr with log P is + 0°64 from a sample 
of 38 stars.* The deviation in the z—¢ distribution is 4-470, and the probability 
of finding logr and log P positively correlated with a correlation coefficient 
0°64 or greater is less than 1 in 2 x 10°. Now the distances calculated by means of 
the intrinsic colours (8) give a correlation coefficient +036, corresponding to 
a deviation in z—{ of 2-240. In this case the probability is less than 1 in 80. 
Consider now the possibility that the cepheids have an artificial correlation between 
log r and log P due to observational selection. We can take account of this 
by supposing that the parent population itself is similarly correlated with a 
correlation coefficient p. If we suppose that p is as large as 0°25, we get a deviation 
in s~¢ of 2-960 in Eggen’s case, and the probability is less than 1 in 650 that by 
chance alone we should find logr and log P positively correlated in the sample 
with a correlation coefficient 0°64. On the other hand, for distances computed 
with the help of equations (8) the deviation in z—{ is o-6ge, and the correlation 
found in the sample is entirely due to chance. In Eggen’s case, however, the 
correlation is still not due to sampling, and the conclusion must be drawn that the 
intrinsic colours (g) which lead to such a large correlation are questionable. 

Evidently the eight cepheids which Eggen originally adopted as unreddened 
are, in fact, appreciably reddened with increasing values of their periods. ‘The 
reason for this can be traced to the coincidence that if the eight cepheids are 
arranged in order of increasing log P, the last six are in order of increasing distance, 
ranging from the nearest at 0-3 kpc to the farthest at 4:0 kpe, and the last five are 
in order of decreasing galactic latitude. ‘The spatial distribution is clearly 
favourable towards appreciable reddening for the longer periods. 

4.2. The colour-spectral type relation.—The spectral classification of the 
cepheids on the Yerkes System has been shown by Code (6) to give period- 
spectrum relations in which the spectral types at maximum light are nearly the 
same for all periods, most, if not all of the well known progressive lateness with 
period shown by the spectral type at median phase being due to the variation at 
minimum light. Code’s assignments of spectral type can be represented 
accurately by the relations 

Sp. (max) = 4°8 + 1°5 log P, 


(11) 
Sp. (min) = 3°0 + 12°3 log P, 


in which the spectral type is in sub-class units measured from type Fo. The 
relation for Sp. (min) has only been calculated for cepheids belonging to Eggen’s 
types A and B since the type C spectra probably obey a slightly different relation 
at minimum light. From equations (8) and (11) we obtain the following colour- 
spectral type relations : 


Types A, Band C: (P-V)max 0°40 +0'12 Sp. (max), 


Types A and B: (P-V) min +.0°0635p. (min). 


These relations are shown in Fig. 4 with the intrinsic colours of selected cepheids 
for which P, V colours and Yerkes spectral classifications are available. The 
colours plotted in the diagram are the observed colours corrected for space 
reddening as in Fig. 1. 

* One of this number, with log P«0o'84, was omitted in the solutions in Section 2 as it wave 
a colour residual greater than three times the root mean square error of an equation of condition 
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- Fic. 4.-—-Relations between intrinsic colour and spectral type for cepheid variables 

ied at maximum and minimum light, 

ap An interesting feature of the relations (12) is that the small range in spectral type 
? which occurs at maximum light for a considerable range in period (for example, 
s the range from log P= 0-2 to 2°0 only involves a change in type from F5 to F8) 


is accompanied by a small range in colour. In Eggen’s work (9, p. 403) a large 
range in colour at maximum light was obtained, but this was due, as we have shown, 
to the neglect of space reddening in establishing the intrinsic colours. 

4.3. Comparison with colours of SMC cepheids.—The P, V colours for 
Small Magellanic Cloud (SMC) cepheids at maximum light are shown in Fig. 2. 
Gascoigne and Kron (11) have drawn attention to the fact that Eggen’s colours 
of type II cepheids in high galactic latitudes appear to be strongly related to the 
colours of the SMC cepheids. According to Eggen (21), the relation for type I 
a cepheids is 


(P-V) max = + + 0™-21 log P. (13) 
from equations (8) and (13) we have 
(P-V) max: Type | minus Type II + o™-og — 0-03 log P, 


the root mean square error in each coefficient being + 0-05 on the assumption 
of error-free values in equation (13). ‘This suggests that there is also a strong 
relationship between the Type I and Type II colours at maximum light, 
sg particularly in the slope of the intrinsic colour relation. The present work 
therefore lends some support to a recent suggestion, due to Eggen (21), that all 
the cepheids may have the same colour at maximum light-—both Types I and II, 
and that the high latitude galactic cepheids of ‘T'ype II and the SMC cepheids 
ave the only unreddened cepheids known. It may, however, be necessary to 
modify this suggestion if the obscuration of the extra-galactic nebulae in directions 
surrounding the Small Cloud in galactic coordinates /= 268°, b= —45° is 


| 
2 | 
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approximately equal to the average obscuration. According to Oort (3, p. 237) 
the average photographic absorption for extra-galactic objects in galactic latitude 
b is given by A,,=o™ 31cosech. If, as seems likely, most of the obscuration 
is due to selective absorption by matter near the galactic plane, the average colour 
excess for objects in the Small Cloud is 0-314/2/4°75 or o-og magnitudes. A 
colour absorption of this order would probably place the intrinsic colours of the 
SMC cepheids beyond the statistical limits set by the uncertainty in the solution 
adopted for the colours of Type I galactic cepheids. Coincidence between the 
two sets of colours would then require more space reddening than is at present 
indicated by the B stars. However, the interesting fact would remain that the 
intrinsic colour relations at maximum light for cepheids belonging to different 
stellar populations have approximately the same coefficient of log P. 


5. Application to galactic rotation 

An examination is now made of the connection between differential galactic 
rotation of the nearby cepheids and the zero-point of the period-luminosity curve, 
using cepheid distances corrected for absorption by means of colour excesses 
based upon the intrinsic colours (8). This can be carried out by means of a 
distance scale-factor similar to that introduced by Camm (25) and Mineur (26). 

5.1. The distance scale-factor.—Let M be the absolute magnitude of a cepheid 
as given by the adopted period-luminosity curve, and suppose that the true 
zero-point is AM magnitudes fainter than adopted, that is, the true absolute 
magnitude is M+AM. If the distance modulus is corrected for absorption 
by means of a colour excess which does not depend upon AM, we have the following 
relations for the distance modulus : 


m,~M=s5logr—5=AM +5 logr,— 5, (14) 
where r and r, are the distances calculated for absolute magnitudes M and M + AM 
respectively. In this case we obtain 

(15) 
where A, the distance scale-factor, is related to the zero-point correction by 
AM = —5§ loga. (16) 


But the colour excesses based upon the intrinsic colours (8) are dependent upon the 
zero-point of the period-luminosity curve. Equation (14) should then be 
5 logr,—5, (17) 


where AE is the increase in colour excesses which would arise if the solution for 
the intrinsic colours (8) had been made with absolute magnitudes M+AM. 
The zero-point correction and the distance scale-factor are then related by 


5 loga 
AM’ (18) 


\M 


It can be shown that the coefficients of log P in the relations (8) are approximately 
independent of AM and that the intrinsic colours are increased by o”-o07 per 
magnitude decrease in absolute brightness. If we adopt \E/ AM~=o-07, equation 
(18) becomes AM=7°5 log A. 

5.2. Oort’s constant and the zero-point correction.—The amplitude of the 
double-wave term in Oort’s approximation to the effect of galactic rotation on the 


wy 
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radial velocities is proportional to the true distance 7,. Since r,~Ar, Oort’s 
equation may be written 
V’ = AAr sin 2(1—1,) cos* b, (19) 

where V’ is the radial velocity corrected for solar motion. The solution of equation 
(19) by least squares gives 
XV'rsin 2(1—/,) cos*b 


(20) 


In this way we determine not A, but AA, where A is related to the zero-point 
correction AM by equations (16) or (18). 

The solution (20) assumes that the longitude of the galactic centre and the 
corrections for solar motion are known. We shall adopt /, = 325°, and corrections 
for the standard value of the local solar motion (27). ‘The data on cepheid 
colours used in Section 3 are here augmented by provisional colour excesses for 
eighteen cepheids kindly made available by Eggen. Of the fifty-six stars for 
which colour excesses and radial velocities are available, the distances of all but 
four are less than 2 kpc. It was considered necessary to impose the condition 
r <2 kpe since Oort’s formula is not accurate for large distances. ‘This reduces 
the sample to fifty-two stars which, although small, is quite representative of 
galactic cepheids in so far as the period distribution and velocity dispersion is 
concerned, ‘The following solution was obtained for AA and the root mean square 
error: 

AA = 19°8 + 2°3, (21) 


the units being km per second per kiloparsec of true distance r,, ‘The root 
mean square error of an observational equation was found to be + 10°8 km/sec, 
and the average velocity residual without regard to sign 8-8 km/sec. The mean 
distance of the sample r=o-87 kpe gives an amplitude AAr = 17-2 km/sec for 
Oort’s double wave (19). 


Taser Il 


Oort's constant and the zero-point correction to the P-L curve 


A= (AA) 
§ 


— 


km, sec/kpe | km/sec/kpe 
14°6 | 
17°0 
19°58 19°58 
24°9 
26°9 
39°6 


The connection between A and AM may be written 
A = (AA) x 10°" (22) 


where a= 5 or 7°5 according as equation (16) or (18) for AM is used. From 
equations (21) and (22) we obtain the values of A given in Table I]. These 
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are the values which the sample would give for Oort’s constant if Baade’s revised 
zero-point of the P-L curve were corrected by \.V magnitudes. ‘The marked 
dependence upon AM demands that some explanation be offered for the fact 
that when AM =o (Baade’s zero-point) we recover the value A = 20 km/sec) kpe 
originally found by Joy (28) using Shapley’s P-L curve (29) for which AM = + 1°5 
magnitudes relative to Baade’s zero-point. If the dependence of the intrinsic 
colours on AM is neglected, the cepheid distances are halved and Oort’s constant 
doubled for AM= +1°5. Thus we should obtain A= 40 km/sec kpe on the 
basis of Shapley’s zero-point, whereas Joy found A=20km/sec/kpe. The 
difference is due to the amount of interstellar absorption allowed for in the two 
investigations. In Joy’s work, the model of uniform absorption was adopted, the 
coefficient of absorption being 0-85 per kpc for all galactic latitudes and longitudes. 
In the present work, the absorption depends upon both latitude and longitude, 
and it is necessary to define average values before a comparison can be made with 
Joy’s absorption. 

Following Stebbins and Huffer (1, p. 247) we shall consider two averages for 
the coefficient of selective absorption given by the colour excesses: c,= XE/Zr, 
and ¢c,=(1/n)LEr. Changing to total photographic absorption, the averages 
for the fifty-two stars used in the above solutions are 1-83 and 2™-30 respectively. 
We can obtain a third estimate of the average absorption by means of the B star 
colour excesses used in the determination of the intrinsic colours of the cepheids. 
Referring to the exponential model of the absorbing medium, the average 
absorption for objects within the distance z from the galactic plane is given by 
(see 4, p. 43) 


n+l Ja\ 


3/6 ann} 


Now the average value of a, for the 37 stars used in the Parenago solution in 
Section 3 is 2-74 per kpc, and of the 52 stars used in the solution (20) there are 
21 with z<50 pe and 40 with z<100 pe, the average value of 2 being 69 pe. 
From equation (23), with 8 = 100 for the absorbing medium, the average absorp- 
tion for objects within 50 pe or 100 pe from the galactic plane is 2"-43 per kpe or 
2™-19 per kpe respectively. ‘These three estimates place the average absorption 
near 2 mag. per kpc compared with Joy’s value o™ 85 per kpe. 

The recovery of Joy’s original value of Oort’s constant in the present work is 
due to the fact that the greater absorption compensates for the increase in the 
distance moduli introduced by Baade’s correction of Shapley's zero-point. ‘The 
effectiveness of this compensation within the context of uniform absorption 
can be seen as follows: Let r be the mean distance of the sample for Baade’s 
zero-point and the coefficient of absorption a, 7, beingthe mean distance for 
Shapley’s zero-point AM magnitudes fainter and for the coefficient of absorption 
a,. Averaging over the stars in the sample we obtain 


M=m+5—s5logr—ar—5C 
where C and C, are corrections for dispersion in distance given by 


('=logr—logr, logr, log 


ay I 
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The equation for the mean distance is then 
slog7, + ar, =5(C —C,)+5logr+ar—AM. 

In this equation we take r=o-87 kpc, AM=+1°5, a,=0™85 per kpc, and 
a=2™ per kpc. If we neglect the term 5(C —C,) on the ground that the effect 
of the above mentioned compensation is to equalize the dispersions, we obtain 
1,073 kpe. Considering the limitations imposed by the model of uniform 
absorption, and the difficulty in obtaining a suitably defined average coefficient 
of absorption, we may conclude that the mean distance of the sample when 
referred to Shapley’s zero-point and Joy’s coefficient of absorption does not differ 
significantly from the mean distance referred to Baade’s zero-point and the higher 
value of the coefficient of absorption. 

Conclusion.—The intrinsic colours of Type I cepheids adopted in the paper 
give an average absorption of about 2 magnitudes per kpc for objects near the 
galactic plane, the result being dependent to some.extent upon how the absorption 
is averaged. Evidence for the high absorption was obtained in 1938 by Oort, 
who found an average coefficient of selective absorption in the galactic plane of 
o™-32 per kpe from the colour excesses of B stars on the system of Stebbins, 
Huffer and Whitford. In Oort’s work, a correction +0™-049 was applied to 
the zero-point on the SHW colour excesses and the minimum value 4,,/E, = 5-5 
was adopted. It now appears that the correction should not be greater than 

+ o™-o4, and the ratio A,,/E, =g seems well-established. Oort’s data then give 
an average coefficient of selective absorption o™26 per kpe for objects within 
100 pe of the galactic plane, and an average photographic absorption 2™-3 per kpc. 
A similarly high value of the average absorption near the galactic plane has been 
advocated by Parenago and van Rhijn. 

Some of the direct consequences of the high absorption are (1) the similarity 
between the colours of cepheid variables in the Small Magellanic Cloud and the 
intrinsic colours of ‘Type I galactic cepheids, the latter being hitherto supposed to 
be redder, and (2) the recovery of Joy’s original value of Oort’s constant when 
Baade’s correction is applied to Shapley’s zero-point of the period—luminosity 
curve for galactic cepheids. 
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SPECTROPHOTOMETRIC MEASUREMENTS OF EARLY-TYPE 
STARS 
I]. ResULTS POR STARS OF TYPE BI* 


W. M. H. Greaves, E. A. Baker and R. Wilson 


(Communicated by the Astronomer Royal for Scotland) 


Summary 


The Edinburgh spectrophotometric programme aims at the measurement 
of equivalent widths of stellar spectral lines with the highest accuracy reason- 
ably attainable using the available equipment. In this paper a general 
description is given of the technique which has been developed ; it supple- 
ments the more detailed description previously published. All the lines 
measured lie in the spectral range 4100 A to 6800 A over which tolerable focus 
can be secured with the spectrograph used. For each star the equivalent 
widths given depend on the mean of five or six negatives and each negative has 
been measured at about 1500 fixed wave-lengths. Experience has confirmed 
that the methods used achieve a substantial gain in accuracy. Systematic 
errors in calibration have been reduced to the order of one per cent. Plate 
grain fluctuations have been reduced by combining the measures for the 
several negatives of each star. A further reduction in accidental fluctuations 
can be secured by averaging for several stars of the same type, and the 
formation of a mean plot for five B1 stars (27 negatives), when combined with 
the information given in Moore’s Tables, has made it easier to choose con- 
tinuous spectrum wave-lengths clear of weak lines with central depressions 
of the order of one per cent. In this way the falsification of the interpolated 
continuum which can result trom the presence of weak lines is reduced 
On the mean Br plot the standard errors of measurement at one wave-length 
are + 0°008 and + 0°006 magnitudes in the blue and red respectively. ‘This 
mean Br plot also reveals the presence of faint blending lines and a reasonably 
satisfactory treatment of blends becomes possible 

Equivalent widths of 28 spectral lines are given for each of seven stars. 
For each line the limits of integration used and the standard error of the 
equivalent width are also given. The standard error depends, inter alia, on 
the range of integration ; for the weaker lines the s.e. is about + 0°07 equivalent 
angstroms 


For lines other than hydrogen or helium there is a systematic dependence 
of equivalent width on stellar luminosity, the equivalent widths being larger 
for the higher luminosities. ‘This progression covers a range of more than two 
to one in equivalent width. For the Balmer lines Ha, Hf and Hy the 


progression is reversed, the equivalent widths being smaller for the higher 
luminosities 


The results for the He 1 lines are of special interest. When the singlet 
and triplet lines 6678 and 5876 are included in a discussion of the material, it 
becomes impossible to reconcile the equivalent widths with the theoretical 
curve of growth. ‘There is a deficiency in the measured equivalent widths of 
6678 and 5876 which seems to increase with decreasing luminosity. It is 
suggested that cyclical transition processes are operative, and are more 
pronounced for the lower luminosities 

The Balmer lines of 8 Cephei have been found to be variable, the variation 
being especially pronounced for Ha. From the present material no statement 
can be made as to whether the variation is periodic or irregular. 


* The full text of this paper appears in Publications of the Royal Observatory, Edinburgh, 
1, No. 6, 1955 
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